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CHAPTER 1 
GEMPAK OVERVIEW 


This document is the user manual for GEMPAK. It is intended to 
help users running GEMPAK programs and to serve as a reference 
guide. 

Information on how to run GEMPAK is contained in Chapter 2. 
Chapter 3 contains descriptions of the variables used by the 
GEMPAK programs. Chapter 4 contains descriptions of the GEMPAK 
programs . 
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CHAPTER 2 
USER INTERFACE 


The Transportable Applications Executive, TAE, is the user inter- 
face for the GEMPAK programs. A simple, optional interface, called 
the no-TAE, is also available. This chapter gives a brief descrip- 
tion of these interfaces. 


2.1 TAE INTERFACE 


The TAE provides menus and tutors which are helpful in running the 
GEMPAK programs. The logical names GEMUSR, GPUSR, and TAE must be 
assigned to point to the GEMPAK, GEMPLT , and TAE software. These 
names may be system logical names or may be assigned in your 
LOGIN.COM file, which resides in your home directory. 

In order to run GEMPAK programs you must enter the TAE. This is 
done by executing the command: 

aGEMUSR : GEMSTRT 

Define the following symbol as a system symbol in your LOGIN. COM : 
GEMPAK :=*- aGEMUSR : GEMSTRT 

This will allow you to get to the main menu by simply typing: 

GEMPAK 

Once you have logged into the TAE, you will be presented with the 
GEMPAK root menu displaying several options. Select the option 
that you want by entering the appropriate number. When you have 
selected a program you will enter tutor mode. After selecting 
par ame ters in the tutor, enter RUN to execute the progr am . T o 
exit the program, type EXIT. If you need help at any time, type 
HELP . 
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USER INTERFACE 


Most of the GEMPAK programs will enter an interactive tutor mode 
after executing the program or encountering an error. This 
dynamic tutor simply provides a way to execute the program again 
with new values for the variables. The dynamic tutor is a 
NOSCREEN tutor. It is functionally the same as the usual SCREEN 
tutor. It is used in GEMPAK so that the terminal screen will not 
be cleared within a program. Instead, the user will be prompted 
wi th 


GEMPAK - PROGRAM> 


Any of the input parameters may be changed. The following commands 
can be entered: 


p a rm - value 

RUN 

LIST 

LIST parm 

DISPLAY 

DISPLAY parm 

SCREEN 

NOSCREEN 

EXIT 

SAVE n ame 
RESTORE name 
HELP 

HELP parm 


change the value of an input parameter 
run the program 

list the parameters with their current values 

list the value of the par ame ter 

list par ame ters, level 1 helps, and values 

list parameter, level 1 help, and value 

enter the SCREEN tutor mode 

enter the NOSCREEN tutor mode 

exit and return to the TAE 

save the program variables in a file, name 
restore the program variables from name 
list help for the progr am 
list help for the par ame ter 


This interactive mode will not be entered if the value of the global 
parameter, SRESPOND, is NO. 


2.2 NO -TAE INTERFACE 


The no-TAE mode is also available to run programs on systems where 
the TAE is not installed. It can also be used to run programs 
directly, without entering the TAE. Any program can be executed 
by typing: 

RUN GEMEXE : PROGRAM 

The program will display a tutor similar to the NOSCREEN dynamic 
tutor in the TAE. All the options mentioned above are available, 
except the SCREEN mode. Parameter values will be retained between 
programs. In the no-TAE, ^RESPOND and SMAPFIL can be set in the 
tutor for any program. This allows a map file to be changed for 
later runs of the program. 
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CHAPTER 3 
GEMPAK VARIABLES 


This chapter describes variables used in the GEMPAK5 programs. The 
values of these variables are retained between programs. 


3 . 1 ADDSTN 


ADDSTN is a logical variable which indicates whether stations 
which are in STNFIL, but not already in the station file, 
should be added to the file. 

Enter YES to add stations. Enter NO to update current stations 
but not add new ones. 


3.2 ANLYSS 


ANLYSS is the information to be stored in the analysis 
block of a grid file. This information is used in the 
objective analysis programs. 

The average station spacing and the grid extend region are 
separated by a slash. The station spacing is in degrees 
of latitude and is used to compute the weighting functions 
in the Barnes analysis programs. The grid extend area is 
specified by four integers, separated by semicolons, which 
are the numbers of grid points to extend the grid left, 
down, right, and up. Only data within the extended grid 
area can be used after the first Barnes pass. 


3-1 



GEMPAK VARIABLES 


The default for the average station spacing is twice the 
grid spacing. The default for the grid extend region is 
2 ; 2 ; 2 ; 2 . 6 


3 . 3 AREA 


AREA is the data area. Only data within the area specified 
will be processed. 


Areas may be defined containing subareas. Subareas must be 
separated by slashes. Each subarea is additive, + , or 
subtractive, - , depending on the first character following 
the slash, with + being the default. Additive subareas add 
stations to the list of valid stations; subtractive subareas 
eliminate stations which were previously valid. 


Subareas may be specified in the following ways 


1 • 1 a 1 1 ; 1 onl ; lat2; lon2 

This defines a 1 a t i tude/1 ong i tude range where 

Clatl, lonl) is the lower left corner and (lat2, lon2) 

is the upper right corner. West longitude is negative. 

#clat;clon;dlat;dlon 

This defines a 1 a t i tude / 1 ong i tude range by the center 
latitude and longitude. The lower left corner is 
iclat-dlat; clon-dlon); the upper right corner is 
Ulat+dlat ; clon+dlon). No corrections are made for 
the poles or the International Date Line. 


2. GEOG 

abbrevi ^ ion for a geographic area defined in 
the GEMPAK geographic table which includes abbreviations 
for states, provinces and countries as well as other 

GEOT S TRl If # ??°? iS entered ’ the user’s geographic table, 
GEOG.TBL, will be searched. * or - after the name may 
be used to reduce/expand the area. 


the 


3 . STN 

lnCS an area centered a station found in 
GEMPAK station table, which currently contains US 
Canadian and Mexican surface stations. * or - after the 
name may be used to r educ e / expand the area. 


4. DSET 

This includes all the stations in the current data set. 
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5 . @ST 


This area includes those stations located in the state, 
province or country defined by ST. Only some countries 
are recognized (US ,CN ,MX,CI , BW.AU) ; other countries may 
be specified using method 6. 


This area includes those stations located in the country 
defined by CN . 


7 . 


aSTNl ;STN2; . . . ; STNn j 

This area includes the stations listed, where 

be a station identifier or a station number. 


STNi may 


8. SHDR: i loval : ihival 

This area defines a range of integer values for the 
station header, SHDR. Valid keywords for SHDR are: 


COUN 

SELV 

SLAT 

SLON 

STAT 

STID 

STNM 


country 

elevation (in meters) 
latitude (in degrees x 100) 

longitude (in degrees x 100, West is negative) 
state 

character identifier 

station number (WO 5-digit number for upper 
air, 6 digits for surf ace- -usua 1 ly the VMO 
5-digit number followed by a 0) 


where COUN, STAT and STID are not very useful, since the 
integer representation of characters is system dependent. 


For example, SELV:0:2000 specifies stations whose 
elevations are less than 2000 meters. 


3 . 4 BORDER 


BORDER is the color, line type and line width of the plot 
background. The three numbers must be s epa r a t ed wi t h 
slashes . 

If the color is 0, no background will be drawn. The defaults 
for color, line type and line width are each 1. 
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3 . 5 CCCMP 


CCCMP specifies the color components. 


When CTYPE - NAME, CCCMP should contain a color 
such as red, blue, maroon, etc. Only the first 
letters are used to search for the name. 


name 

three 


When CTYPE - RGB or CTYPE - HLS , 
values between 0 and 1. should be 
slashes . 


three real-number 
entered, separated 


by 


When CTYPE - INIT, all the device colors are 
to a standard set of colors. 


initialized 


Co lor table 


NAME 

RED 

GREEN 

BLUE 

RED 

1 .000 

.000 

.000 

GREen 

.000 

1.000 

.000 

BLUe 

.000 

.000 

1.000 

YEL 1 ow 

1 .000 

1.000 

.000 

CYAn 

.000 

1 .000 

1 .000 

MAGen t a 

1.000 

.000 

1.000 

WHI t e 

1 .000 

1.000 

1 .000 

BLAck 

.000 

.000 

.000 

GRAy 

.500 

.500 

.500 

VANi 11a 

1 .000 

. 900 

. 750 

ORAng e 

1 .000 

.500 

. 000 

DKOr ang e 

. 863 

.471 

. 000 

GBR own 

. 500 

.281 

.281 

BROwn 

.500 

.294 

.196 

BE I ge 

1.000 

.672 

. 500 

SANd 

. 851 

. 600 

.298 

MUD 

. 395 

.100 

. 160 

AQUa 

.000 

. 718 

. 718 

SKY 

.000 

.640 

1 .000 

LTB 1 u e 

.000 

. 500 

1 .000 

NAVy 

.000 

.000 

. 500 

DKBlue 

.000 

.000 

. 750 

INDi g o 

.578 

.000 

1.000 

LTG r e e n 

.600 

1 .000 

.000 

AVOcado 

.250 

. 500 

. 125 

DKGreen 

.000 

.562 

.000 

TAN 

.750 

.500 

.000 

GOLd 

1 .000 

. 750 

.000 

BLOnd 

.900 

. 820 

. 390 

LEMon 

. 900 

1.000 

.250 


HUE 

.000 
.333 
.667 
. 167 
.500 
.833 
.000 
.000 
.000 
. 102 
.083 
.092 
.000 
.051 
.055 
.092 
.970 
.500 
.558 
.583 
.667 
.667 
.764 
.232 
.280 
. 333 
.113 
. 129 
. 143 
.186 


LIGHT 

.333 
.333 
.333 
.667 
.667 
.667 

1 .000 
.000 
.500 
. 883 
.500 
.445 
. 354 
. 330 
. 724 
.583 
.218 
.479 
.547 
.500 
. 167 
.250 
.526 
.533 
.292 
.187 
.416 
.583 
. 703 
.717 


SATUR 

1.000 
1 .000 

1 .000 

1.000 

1 .000 

1 .000 
.000 
.000 
.000 
. 320 
.816 
.710 
.216 
.269 
.464 
.456 
. 308 
.678 
. 831 
.816 
.632 
.816 
. 821 
. 824 
. 343 
. 679 
.632 
. 862 
.490 
. 738 
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PINk 

1.000 

.600 

DKP i nk 

1 .000 

. 500 

MARoon 

.578 

.000 

APRi cot 

1 .000 

. 700 

VIOlet 

.625 

. 250 

PURp 1 e 

.500 

.000 


.600 

.000 

.733 

.426 

.500 

.000 

.667 

.500 

.000 

.000 

.193 

.691 

.375 

.000 

.692 

.553 

.625 

.833 

.500 

.354 

.500 

. 833 

.333 

.500 


3.6 CINT 


CINT is the contour interval, minimum and maximum values 
separated by slashes. The scale factor is now entered in 
the variable SCALE. 

The contour interval may be any real number. If it is not 
specified or if the value is 0, the program will select a 
contour interval which will generate 5 to 10 levels. 

The minimum and maximum values specify the range of data to 
use in selecting contour levels. If either value is not 
specified, the value will be obtained from the range of 
values in the dataset. If the minimum and maximum are 
equal, that value will be used and only one contour level 
will be dr awn . 

A list of two or more contour levels may be entered using 
semicolons to separate the individual values. In this case, 
the minimum and maximum are ignored. 


3 . 7 CLEAR 


CLEAR determines if the screen is to be cleared before 
plotting. 

Enter YES to clear screen or NO to leave current 
graphics . 
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3 . 8 COLORS 


COLORS specifies a list of color numbers which must be 
separated using semicolons. 

The color list will be repeated, if necessary, to determine 
colors for all the parameters or contours to be drawn. The 
colors corresponding to the color numbers used in each entry 
are device dependent. 

If no color is input, the default is color 1. If a color 
number of 0 is given, no plotting will be done. 

A specific color may be assigned to a color number by 
entering -NAME after the color number. The first three 
letters in NAME will be used to identify the color in 
the GEMPAK color table. For example, COLORS-l-red ; 2-blue , 
will set color number 1 to red and color number 2 to blue 
If the last character in COLORS is a ?, then the current 
color names will be listed. 


3.9 CPYFIL 


CPYFIL identifies the location of the grid navigation and 
analysis information to be stored in a grid file. Three 
options are available: 

1. If CPYFIL is blank, the information is taken from the 
inputs for PROJ , GRDAREA, KXKY, and ANLYSS. 

2. If CPYFIL begins with a #, the information is read 

from the line in the grid navigation table which 
has the grid name or number corresponding to the 
rest of CPYFIL. * 

3. If CPYFIL is the name of a current grid file, the 
grid and navigation information will be copied from 
that file. 


Note that if CPYFIL is not blank, 
GRDAREA , and KXKY will be ignored. 


the values of PROJ, 
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3.10 CTLFLG 


CTLFLG is a logical flag which indicates 
characters are included in a raw surface 
decoded . 


wh ether control 
data set to be 


The raw surface data set must contain either bulletins 
from the Domestic Data Service or single 8 1 a 1 ‘ ° n . I e \?l l [ { e 
If the data set contains bulletins, the records 
must be 80 bytes each, and each bui 1 e t in mus t beg in with a 
<CNTL>A in column 1. If the input cent a > ns r epo: r s , e ach 
report must be written on a single line. CTLFLG IS J 

variable indicating whether control characters aje present. 
If CTLFLG - YES, the input file must contain bulletin . 
CTLFLG - NO, the file must have reports. 


3.11 CTYPE 


CTYPE specifies the manner in which the colors are to be 
changed : 

INIT -- colors are initialized to their default values. 
NAME -- color name is given in CCCMP. 

RGB - - RGB c omp onents are given in CCCMP . 

HLS -- HLS components are given in CCCMP. 


3.12 CURVE 


CURVE is a number corresponding to the method to be used 
to fit the curve. Generally, CURVE should be set to 2. 

The curve types currently available are: 

1 - piecewise straight line 

2 - cubi c spl ine 

21 - cubic spline with linear ends 

22 - cubic spline with parabolic ends 

23 - cubic spline with extrapolated ends 

Note that types 2 and 21 are the same. Types 22 and 23 produce 
slightly different results at the ends of the splines. 
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3.13 CXSTNS 


In SNCROSS , CXSTNS is the list of stations, separated with 
semicolons, to be used for the cross section. Stations may 
be entered using either character or numeric identifiers. 

In GDCROSS, CXSTNS contains the endpoints of the cross- 
section line separated using a >. Each endpoint may be 
defined as foil ows : J 

1. a station character or numeric identifier; 

2. a latitude and longitude pair separated by a semicolon; 

3. an e followed by a grid x and y coordinate pair separated 
using a semicolon, e.g., @1 . 5 ; 2 . 3>@25 . 9 ; 30 . 


3.14 DATTIM 


DATTIM is the date and time to 

The standard format for DATTIM 
YYMMDD/HHMM, where: 


be used by GEMPAK programs, 
is an 11-character string 


YY is the last two digits of the 
MM is the month 
DD is the day of the mo nth 
/ is the date and time separator 
HH is the hou r 

MM is the minutes past the hour 


year (ex: 84 for 1984) 
(ex : 03 for March) 

(ex: 14 for 2:00 pm) 


The part before the / is the DATE • 
the TIME. 


the part after the / is 


DATTIM may be abbreviated. If the input has no / it is 

e i thTr^TIMF be *£ IrV™ ** abbreviated version of 

either TIME or DATE is assumed to be the part closest to the 

Um, U ISI'fU, F S,i,ndard , time iS ° b,ained fr °“ the last 

» th, tV,, ' Iample ’ if ">• >*** time in the file 

84U515/1200 the following translations will be done: 


13/11 

13 

13/ 

0412/1300 


> 840513/1100 

> 840515/1300 

> 840513/1200 

> 840412/1300 
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A list of times may be entered for DATTIM. Times in the 
list must be separated with semicolons. For example: 

13/11 ;0412/13 ;0515/6 

In the above example, the times are: 

840513/1100 
840412/1300 and 
840515/0600 . 

DATTIM may also be entered as a range. The first and last 
times must be separated by a minus sign, for example: 

13/11-15/14 

A range with an increment may also be entered using minus 
signs as separators. The format of the increment is hhhmm. 
If the increment is one or two digits, it is assumed to be 
in hours. 

The following options are also valid for DATTIM: 

LAST -- the last time in the file 

LIST -- lists all times and waits for input 

ALL -- all the times in the file 

/ALL -- all the times for a single DATE. 

If /ALL is entered, all the times for the most recent date 
are processed. If /ALL is appended to a specific date, 
then all the times for that date will be processed. 


3.15 DELTAN 


DELTAN is the average station spacing in degrees of latitude. 
The Barnes objective analysis programs use this number to 
compute weights for data interpolation. 


3.16 DELTAX 


DELTAX is the spacing be t we en grid points in the x direction 
on CED grids. This value is in degrees longitude. 
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3.17 DELTAY 


DELTAY is the spacing between grid points in the y direction 
on CED grids. This value is in degrees latitude. 


3.18 DEVICE 


DEVICE is the graphics device. 

The following valid graphics devices may be available: 

VT : VTIOO Re t r o g r aph i c s terminal 

DI : Ditto graphics terminal 

QM : QMS laser printer 

CT : Tektronix 4107 color terminal 

II : IIS Image Terminal 

RE : VT240 REGIS graphics device 

LA : LA100 printer plotter 

LN : LN 03 laser printer 

HP : Hewlett Packard Pen Plotter 

MC : Macintosh with Versaterm 

SG : Silicon Graphics Personal Iris 

AP : Ap olio 

GKS : GKS standard 

BGL : BGL printer 


3.19 DTAAREA 


DTAAREA defines the area over which station data will be 
input to the Barnes objective analysis. Since data must be 
interpolated from the first pass grid back to stations, only 
data wi thin the EXTEND area wi 11 be used after the first 
pass. If the DTAAREA is not specified by the user, it will 
default to the EXTEND area in OAGRID or to the data area 
stored in the grid file analysis block in OABSFC and OABSND . 


DTAAREA can be specified in three ways: 

1. 1 a 1 1 ; 1 on 1 ; 1 a 1 2 ; 1 on2 

Th is defines a latitude/longitude range wh ere 
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( 1 a 1 1 , 1 onl ) is the lower 

is the upper right corner. 


left corner and (lat2, lon2) 
West longitude is negative. 


lat;clon;dlat;dlon 

This defines a 1 a t i tude / 1 ong i tude range by the center 
latitude and longitude. The lower left corner is 
(clat-dlat; clon-dlon); the upper right corner is 
(clat+dlat ; clon+dlon). No corrections are made for 
the poles or the International Date Line. 


This is an abbreviation for a geographic area define 
in the GEMPAK geographic table which includes 
abbreviations for states, provinces and countries, as 
well as other names. If #GEOG is entered, the user s 
geographic table, GEOG.TBL, will be searched. *, - 

after the name may be used to reduce/expand the area. 


3 . STN 

This defines an area centered on a station found in the 
GEMPAK station table, which currently contains US and 
Canadian stations. *, - after the name may be used 

to reduce/expand the area. 


Note that the other ways 
DTAAREA in OAGRID but are 


of specifying AREA are not valid for 
valid in both OABSFC and OABSND. 


For all projections, the lat/lon corners defined will be 
used exactly. For p r o j e c t i ons wh i ch include the pole, if 
latl - 1 a 1 2 and lonl - lon2, then latl will specify the 
range of data from the pole and lonl will specify the 
central longitude. West longitude is negative. 

Note that the stations actually used for the data area may be 
changed in the objective analysis programs using DTAAREA. 


3.20 EXTEND 


EXTEND specifies the numbers of grid points beyond the 
GRDAREA which define the grid extend area in the Barnes 
objective analysis. The first pass is computed on the 
extend area to reduce edge effects on the GRDAREA. 

EXTEND is specified as four integers, which are the number 
of grid points to extend the grid left, down, right, and 
up. The values are separated by semicolons. The default 
for EXTEND is 2;2;2;2. 
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3.21 FILTER 


FILTER is a logical flag used to filter data to eliminate 
plotting of overlapping stations. 

If FILTER is set to YES, the data will be filtered. 
Otherwise, all data will be plotted. 


3.22 GAMMA 


GAMMA, the convergence parameter, is a multiplier for the 
weight and search radius for passes after the first pass of 
the Barnes analysis programs. GAMMA must be within the range 
0 ~ , value outside this range will default to a value 

of .3. If GAMMA is 0, the number of passes will be set to 1 
The recommended value for GAMMA is .3. 


3.23 GAREA 


GAREA is the graphics area. This is the area which will be 

displayed on a graphics device. 

GAREA can be specified in three ways: 

1 . lat 1 ; 1 onl ; 1 a t2 ; 1 on2 

This defines a 1 a t i t ude / 1 ong i tude range where 

(latl, 1 onl ) is the lower left corner and (lat2, lon2) 

is the upper right corner. West longitude is negative 

#clat;clon;dlat;dlon 

This defines a 1 a t i t ud e / 1 ong i t ude range by the center 
latitude and longitude. The lower left corner is 
(clat-dlat; clon-dlon); the upper right corner is 
(clat+dlat; clon+dlon). No corrections are made for 
the poles or the International Date Line. 

2. GEOG 

This is an abbreviation for a geographic area defined 
in the GEMPAK geographic table which includes 
abbreviations for states, provinces and countries, as 
well as other names. If #GEOG is entered, the user’s 
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geographic table, GEOG.TBL, will be searched. *» 
after the name may be used to reduce/expand the area. 

3 . STN 

This defines an area centered on a station found in the 
GEMPAK station table, which currently contains US and 
Canadian stations. *, - after the name may be used 

to r educ e / expand the area. 

Note that the other ways of specifying AREA are not valid for 
GAREA. 

For all projections, the lat/lon corners defined will be 
used exactly. For projections which include the pole, if 
1 a 1 1 - 1 a t2 and lonl - lon2, then latl will specify the 
range of data from the pole and lonl will specify the 
central longitude. West longitude is negative. 

If satellite projections are used, PROJ must be set to AOI 
for AOIPS/2 navigation or to NPG for Naval Postgraduate 
School navigation and GAREA must be the file name of an AOIPS 
image for AOI or the navigation base file name for NPG. 


3.24 GDATTIM 


GDATTIM is the date/ time for the grid. 

Grids may contain two date/ time fields, in which case the 
fields mu st be separated using a colon. If the grid to 
be selected contains only one time, the colon and second 
time may be omitted. 

The standard format for grid times is a character string 


YYMMDD / HHMM t h h hmm 


wh ere: 

YYMMDD is the year, month, day 
HHMM is the hour, minute 
/ is the date and time separator 

t is the type ( F-forecast A-analysis 

h h hmm is the forecast hour, minute 


G-guess ) 


This format allows the two time fields to be used, for 
example, for the difference of two times. If only a single 
time is needed, TIME (2) is blank. If t is blank, an analysis 
grid is assumed. If hhhmm is blank, 00000 is assumed. If 
hh hmm has one or two digits, they represent hours. With 
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three or more digits, zeros will be added at the beginning 
of the f i e Id . 

The field YYhMDD is called the date; HHMM is the time; 
t is the type; hhhnm is the forecast time. The fields may 
be abbreviated. If the input has no /, it is assumed to be 
the time part. Any abbreviated version of either date or 
time is assumed to be the part closest to the /. The rest 
of the date and time is obtained from the last time in the 
file. If the type is missing, then A is assumed, and a 
missing forecast time is replaced by 00000. 

The type and forecast time fields are used with forecast 
model data. If GDATTIM is 8 8083 1 /0000F24 , the grid to be 
found is the 24-hour forecast from the model run at 00Z on 
Aug 31. 

Examples (with the last time - 851205/1600 ): 

GDATTIM TIME1 TIME2 


851205/1200 851205/1200 

12:15 851205/1200 851205/1500 

LAST 851205/1600 

/ 00F24 85 1 205/0000F24 

The value in GDATTIM may be overridden by specifying "GDATTIM 
with the grids to be found. For example, 

GFUNC — SUB (TMPF"2 8/12 ,TMPF"27/12) will compute the same 
layer difference as GFUNC - LDF (1MPF) with GDATTIM-28/12 : 27/12 . 


Note that a grid from an objective analysis will have forecast 
type and time A00000 , which is the default for a blank, forecast 
type and time. 


3.25 GDEFIL 


GDEFIL is the name of the grid edit file which will be used 
to update a grid file. The edit file must contain only 
c omp lete grids. 

GDEFIL is a text file which may be created using the program 
GDLIST with F as an output device. A text editor may be 
used to create or change the grid edit file. 

GDEFIL must contain the following information before the 
c omp lete grid data. 
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The time, level, vertical coordinate and parameter name MUST 
appear on the same line in the order given. The time must 
be a fully qualified GEMPAK time. The parameter name is the 

name of the grid. 


The first and last row and column information must appear 
on one line. Column i nf o rma t i on mus t be preceded by the 
word COLUMNS:. Row information must be preceded by the 
word RCW: . The first row and column must be 1. The 
last row and column must be the same as the grid size in 
the file. For example: 


COLUMNS: 1 8 RCWS : 1 6 

If the string FACTOR:, appears and is followed by an integer, 
the data will be divided by 10**SCALE before being stored in 
the grid file. In the following example, the data will be 
divided by 10**5. 


SCALE FACTOR: 10**5 

The string COLUMN: signals the beginning of the grid data. 

If the required grid information has not already been found, 
an error will result. 

The data for the grid follows the grid keywords. The data 
for a grid is listed from the top (last) row to the bottom 
row in the grid. The beginning of each row of data must be 
preceded by the string RCW. Data consists of numeric 
values for each data point; missing data should be entered 
as -9999.0. 


3.26 GDFILE 


GDFILE is the name of the file which contains gridded data 


3.27 GDNUM 


GDNUM allows the user to select grids by number. 
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GDNUM ma y be specified in the foil owi n g wa y s ; 

GDNUM — range 

The first and last grid numbers in the range are 
separated by a If an increment is included, 

it will be ignored . 

GDNUM — list 

The grids numbers in the list are be separated 
using semicolons. 

GDNUM - ALL 

All the grids in the file will be included. 

GDNUM - LIST 

All the grids in the file will be listed. The user 
will be prompted to select the grids to be included. 


3.28 GDOUTF 


GDOUTF is the output grid data file name. 

GDOUTF is used in programs which create a new grid data 
file. It is also used in programs that move data from an 
input file to an output file. In that case, GDFILE is the 
original file. 


3.29 GFUNC 


See the GPARM documentation. 


3.30 GLEVEL 


GLEVEL is the vertical level for the grid. 

Grids may contain two levels separated by a colon. If 
the grid to be selected contains only one level, the colon 
and second level may be omitted. In this case, the second 
level is stored in the grid file as -1. 
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Note that the vertical coordinate system for GLEVEL is 
specified by GVCORD. 

The value in GLEVEL may be overridden by specifying ©GLEVEL 
with the grids to be found. For example, the following two 
computations are identical: 

GFUNC * SUB (TMPF@850 ,TMPF@500) 

GFUNC - LDF (TMPF ) and GLEVEL - 850:500. 


3.31 GNAME 


GNAME is the geographic name abbreviation and full name 
to be saved in the geographic table. These names must be 
separated with a slash. 

The abbreviation may later be input as an area name in AREA 
or GAREA. The full name will be saved in the table but will 
not be used. 


3.32 GPACK 


GPACK is the number of bits and packing type to be used to 
pack the grid data in the output file. The two values 
mu st be separated using a slash. 

The valid packing types are: 

NONE No packing 

GRIB Data is packed in GEMPAK GRIB format 

DIF Data is packed in GEMPAK DIF format 

DEC Data is packed in GEMPAK GRIB format. 


If the number of bits is 32, the data will be stored as 
real numbers, using packing type NONE. 

If the number of bits is less than 32, the packing type 
specified will be used. If no packing type is given, the 
data will be packed using the GEMPAK GRIB scheme. 


If the packing type is DEC, the first value is the precision 
used to pack the data instead of the number of bits. 
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In general, data should be stored using packing type NONE 
or in GRIB format using 16 bits. Care should be taken using 
the other packing schemes. 


3.33 GPARM (GFUNC , GVECT ) 


GFUNC and GVECT are the scalar and vector grid functions. 

They are input as nested strings of operators and 
operands. The operand list for an operator is enclosed in 
parentheses or square brackets, with operands separated by 
semicolons or commas. 

The following grids will be computed automatically from grids 
in the grid file: 

TEMP DWPC MIXR* THTA* DRCT 

TMPC EWPF SMXR STHA SPED 

TMPF DWPT RELH* THTE* STHE 

where * indicates names which may be used as operators also. 
Other special scalar parameter names denote constant value 
grids : 

DTR Conversion factor for degrees to radians - 

PI / 180 

E Base of natural logarithms - 2.71828182 

GRAVTY Gravitational constant - 9.80616 (note spelling) 
KAPPA Gas cons tant/ spec i f i c heat - 2/7 
PI 3.14159265 

RTD Conversion factor for radians to degrees = 

180 / PI 

nnn Any number (e.g., 2, -10.2) 

Another class of special parameter names depends on the grid 
navi ga t i on : 

CORL Coriolis force = 2. * OMEGA * SIN ( LATR ) 

LATR Latitude in radians 

LONR Longitude in radians 

XVAL Value of the x coordinate in graph coordinates 

YVAL Value of the y coordinate in graph coordinates 

A grid may be identified by its number in the grid file by 
prefixing the number with the symbol #, e.g., #5. 

Standard vector grids are: 
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OBS Observed wind (with conversion, if needed) 

GEO* Geostrophic wind 

AGE* Ageostrophic wind 

ISAL* Isallobaric wind 

THRM* Thermal wind 

where * indicates names that also may be used as operators. 

Time, level, and vertical coordinate as specified through 
the user interface may be overridden by in-line parameters: 

"time ©level %ivcord 

appended to an operand in any combination. 

Grid operators may be nested. Note that layer and time range 
operators expect operands read directly from the grid file. 

In the following list of diagnostic operators, scalar 
operands are named Si and vector operands are Vi. Vector 
components are denoted by u and v. All meteorological gri s 
are in MKS units, except as noted. POL following the 
description indicates that the computation currently can 
only be performed on polar (R, THETA) grids. In the 
trigonometric functions, the angles are expressed in radians. 

SCALAR OUTPUT GRID 


ABS 

(S) 

Absolute value 

ACOS 

(S) 

Arc cosine function 

AS IN 

(S) 

Arc sine function 

ATAN 

(s) 

Arc tangent function 

ATN2 

(SI ,S2) 

Arc tangent function 

COS 

(S) 

Co sine f unc t i on 

EXP 

(SI ,S2) 

Exponential to real 

EXP I 

(SI ,S2) 

Exponential to integer [uses NINT(S2)J 

LN 

(s) 

Natural logarithm 

LOG 

(s) 

Base 10 logarithm 

SIN 

(SI ,S2) 

Sine function 

SQRT 

(S) 

Squa r e root 

TAN 

(S) 

Tangent function 

ADD 

(SI ,S2) 

Addition 

MUL 

(SI ,S2) 

Mu ltiplication 

QUO 

(SI ,S2) 

Division 

SUB 

(SI ,S2) 

Sub t racti on 

ADV 

(s,v) 

Adve c t i on 

AVG 

(SI ,S2) 

Av e r a g e 

AVOR 

(V) 

Absolute vorticity 

DDEN 

(PRES, TMPC) 

Density of dry air 
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DDR 

(S) 

DDT 

(S) 

DDX 

(S) 

DDY 

(s) 

DEF 

(V) 

DIRN 

(V) 

DIRR 

(V) 

DIV 

(V) 

DOT 

(VI ,V2) 

DTH 

(S) 

FCNT 

(s) 

FRNT 

(THTA.V) 

JCBN 

(SI ,S2) 

KNTS 

(S) 

LAP 

(S) 

LAV 

(S) 

LDF 

(S) 

MAG 

(V) 

MASS 


MDIV 

(V) 

MIXR 

(EWPC , PRES ) 

MSDV 

(s,v) 

NORM 

(V) 

PLAT 

(S) 

PLON 

(s) 

POLF 

(S) 

PVOR 

(V) 

RELH 

(TEMP.EWPT) 

RICH 

(V) 

ROSS 

(VI ,V2) 

SAVG 

(S) 

SAVS 

(S) 

SDIV 

(s,v) 

SHR 

(V) 

SM5S 

(s) 

STAB 

(TMPC) 

STR 

(V) 

TANG 

(V) 

TAV 

(s) 

TDF 

(s) 


Partial derivative with respect to R (POL) 
Time derivative 

Partial derivative with respect to X 

Partial derivative with respect to Y 

Total deformation 

Direction relative to north 

Direction relative to grid 

Divergence 

Vector dot product 

Partial deriv. with respect to THETA (POL) 

Coriolis force at grid center (POL) 

Frontogenes i s 

Jacobian determinant 

Convert meters / second to knots 

Laplacian operator 

Layer average (2 levels) 

Layer difference (2 levels) 

Magnitude of a vector 

Mass per unit volume in a layer from PRES 
Lay e r - a ve r ag e mass divergence 
Mixing ratio g/g internally, g/kg on output 
Layer-avg. mass-scalar flux divergence 
Normal component (for cross sections) 
Latitude at each point (POL) 

Longitude at each point (POL) 

Coriolis force at each point (POL) 

Potential vorticity in a layer 
Relative humidity 

Richardson stability number in a layer 

Rossby number 

Average over whole grid 

A v erage over subset grid 

Flux divergence of a scalar 

Shearing deformation 

5 -point smoother 

Lapse rate over a layer in K/km 
Stretching deformation 

Tangential component (for cross sections) 
Time average 
Time difference 


(TMPC, PRES) Potential temperature 

( PRES , TMPC , EWPC) Equivalent potential temperature 
(V) North relative u component 

(V) Grid relative u c omponen t 

(V) North relative v c omp o n e n t 

(V) Vor t i c i ty 

(V) Grid relative v component 

( S ) Av erage along a grid r ow 

( S ) S um along a grid r ow 

(S) Average along a grid column 

(S) Sum along a grid column 
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VECTOR OUTPUT GRID 

AGE 

(S) 

CIRC 

(V,S) I 

GEO 

(S) I 

GRAD 

(S) 1 

I NAD 

(VI ,V2) 

I SAL 

(s) 

KNTV 

(V) < 

LTRN 

(s,v) 

OBS 

I 

NRMV 

(V) 

ROT 

( ang 1 e , V) 1 

SMUL 

(S,V) ] 

SM5V 

(V) 

TNGV 

(V) 

THRM 

(S) 

VADD 

(VI ,V2) 

VAVE 

(V) 

VAVS 

(V) 

VECN 

(SI ,S2) 

VECR 

(SI ,S2) 

VLAV 

(V) 

VLDF 

(V) 

VMUL 

(V1.V2) 

VQUO 

(VI ,V2) 

VSUB 

(VI ,V2) 


kno t s 

of a scalar 
cross sections) 


Circulation (for cross section) 

Geostrophic wind 
Gradient of a scalar 
Inertial advective wind 
Isallobaric wind 
Convert me t e r s / s e cond to 
Lay e r - ave r ag ed transport 
Observed wind 
Vector normal wind (for 
Coordinate rotation 

Multiply a vector’s components by a scalar 
5 -point smoother 

Vector tangential wind (for cross 
Thermal wind over a layer 
Ad d the c omp onents of two vectors 
Average vector over whole grid 
Average vector over subset grid 
Create vector from north relative 
Create vector fr om grid relative 
Layer average for a vector 
Layer difference for a vector 
Multiply the components of two vectors 
Divide the components of two vectors 
Subtract the components of two vectors 


sections) 


c omp onents 
c omp onents 


3.34 GPOINT 

GPOINT is the grid location to be used for the plot. 

GPOINT can be entered in the following ways: 

1. station numeric of character identifier 

2. a latitude and longitude pair separated with a semicolon 

3. @ followed by a grid x and y coordinate pair separated 
with a semicolon. 

If necessary, the grid data will be interpolated using a 
bilinear interpolation to the point specified. 

Ex amp 1 e s : 

GPOINT = al;l.S Grid point (1,1.5) 


3-21 



GEMPAK VARIABLES 


GPOINT — 30;-120 Grid point at 30 latitude and 

-120 longitude 

GPOINT - 5;5 Grid point at 5 latitude and 

5 longitude. Note that this is 
probably an error and the input 
should have been a5 ; 5 

GPOINT - BWI Point located at Ba 1 1 imo r e -Wa sh i ng t on 

International Airport. 


3.35 GRDAREA 


GRDAREA specifies the area to be covered by the grid. This 
area is combined with PROJ to define the region over which 
the grid is evenly spaced. 

GRDAREA can be specified in three ways: 

1. 1 a 1 1 ; 1 onl ; 1 a 1 2 ; 1 on2 

This defines a 1 a t i tude / 1 ong i tude range where 

Clatl, lonl) is the lower left corner and (lat2, lon2) 

is the upper right corner. West longitude is negative. 

#clat ;clon;dlat ;dlon 

This defines a 1 a t i t ud e / 1 ong i tude range by the center 
latitude and longitude. The lower left corner is 
(clat-dlat; clon-dlon); the upper right corner is 
(clat+dlat; clon+dlon). No corrections are made for 
the poles or the International Date Line. 

2. GEOG 

This is an abbreviation for a geographic area defined 
in the GEMPAK geographic table which includes 
abbreviations for states, provinces and countries, as 
well as other names. If #GEOG is entered, the user’s 
geographic table, GEOG.TBL, will be searched. *, - 
after the name may be used to reduce/expand the area. 

3. STN 

This defines an area centered on a station found in the 
GEMPAK station table, which currently contains US and 
Canadian stations. *, - after the n ame ma y be used 

to r educ e / expand the area. 

Note that the other ways of specifying AREA are not valid for 
GRDAREA. 
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For all projections, the lat/lon corners defined will be 
used exactly. For projections wh ich include the pole, if 
1 a 1 1 = 1 a 1 2 and lonl = lon2, then latl will specify the 
range of data from the pole and lonl will specify the 
central longitude. West longitude is negative. 


3.36 GRDNAM 


GRDNAM is the parameter name for the output grid. If 
this name is blank, the default name generated by the grid 
diagnostics package will be used. 


3.3 7 GUESS 


GUESS contains the information to use as a first guess for 
objective analysis programs. The name of the grid file 
containing the guess field and the time to be used to 
extract the field mu st be entered using * as a separator. 


3.38 GVCORD 


GVCORD is the vertical coordinate of the grid to be selected. 
The valid values are: 


NONE 

PRES 

THTA 

IIGHT 


for 

data 

data 

data 


surface data 

in pressure coordinates (millibars) 
in isentropic coordinates (Kelvin) 
in height coordinates (meters) 


The value in GVCORD may be overridden by specifying 9bGVCORD 
with the grids to be found. For example. 


GFUNC = SUB ( TMPC® 8 50%PRES , TMPC® 1 5009bHGHT ) 

will compute the difference between temperatures on the 8 50- 
mb level and the 1500- meter level. 
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3 . 39 GVECT 


See the GPARM documentation. 


3.40 IDNTYP 


IDNTYP is the type of station identifier to use for input or 
output. The valid values are STID and STNM. STID is used 
to specify station character identifiers; STNM specifies 
station numbe r s . 

If the value in IDNTYP is not STNM, the default of STID 
will be used. 


3.41 KXKY 


KXKY specifies the size of a grid. Two numbers separated 
by a semicolon must be input. These numbers are KX, the 
number of grid points in the x direction, and KY, the number 
of grid points in the y direction. 

If the projection is CED, these numbers may be DELTAX and 
DELTAY, the grid spacing in degrees in the x and y 
directions. If the spacing is input, the first character 
i n KXKY mu s t be # . 


3.42 LATLON 


LATLON specifies the latitude and longitude grid lines 
to be drawn. The line color, line type, line width, label 
frequency and increment are separated by slashes. The 
increment specification contains the latitude and longitude 
increments in degrees separated by a semicolon. 

If the color is 0 or LATLON is blank, grid lines are not 
drawn. The latitude and longitude increments will default 
to 10 degrees. The label frequency defaults to 1. 
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Ex amp 1 e s : 

LATLON - 1 

Lat/lon grid lines are drawn every 10 degrees 
in color 1 using line type 1 and line width 1. 

LATLON - 4/8/3/2/5;5 

Lat/lon grid lines are drawn every 5 degrees 
in color 4 using line type 8 and line width 3. 
Every other line is labelled. 


3.43 LEVELS 


LEVELS is the variable specifying vertical levels to be 
extracted from the data set. The coordinate system for the 
levels is specified in the variable VCOORD as PRES, HGHT or 
THTA. 


LEVELS may be a list separated by semicolons. The following 
items may be included in the list: 

a single level; 

MAN for the mandatory levels below 100 mb; 

VAS for the standard VAS levels ; 

a range of levels with an increment separated by - . 

The following items are also valid, provided they are not 
part of a list: 

ALL for all levels; 

a range of levels without an increment. 

SFC or 0 may be entered for surface data. TOP or -1 is the 
top level at the station. These values may not be entered 
for a range with an increment. 


3.44 LINE 


LINE is the color, line type, line width, and line label 
freqency separated by slashes. The individual values in 
each group are separated by semicolons. For example, 
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LINE - 1 ;2/7;8/4;5/2 

assigns colors 1 and 2 to alternate lines, dashing patterns 
7 and 8 to alternate lines, line widths 4 and 5 to alternate 
lines, and labels every other line. 

If no color is input, the default is color 1. If a color 

number of 0 is given, no plotting will be done. 

A specific color may be assigned to a color number by 
entering =NAME after the color number. The first three 
letters in NAME will be used to identify the color in 
the GEMPAK color table. For example, COLORS=*l = r ed ; 2 = blue , 
will set color number 1 to red and color number 2 to blue. 

If the last character in COLORS is a ?, then the current 
color names will be listed. 

If line type or width is zero, the current default will be 
used. There are ten distinct line types: 

1 - solid 

2 - shor t dashed 

3 - medium dashed 

4 - long dash short dash 

5 - long dash 

6 - long dash three short dashes 

7 - long dash dot 

8 - long dash three dots 

9 - me d i um dash dot 
10 - dotted 

These patterns can be expanded or compressed by prefixing the 
single digit with a number from 1 to 9. Aprefix of 1 
compresses the pattern, 2 is the default and prefixes 3 — 9 
expand the basic pattern. For example, 32 expands line type 2 
while 12 compresses the same pattern. 

If the line type and line width are entered as more than 
one number separated by semicolons, then successive lines will 
have different characteristics corresponding to the different 
numbers. The pattern established will be repeated to 
accommodate all lines. 

If the label info rma tion is a single numb er, n, then every 
nth line will be labeled. If a label is a sequence of n umb e r s 
separated by semicolons, lines corresponding to positive 
integers will be labeled. Th e pattern established will be 
repeated to acc ommo date all lines. 
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3.45 LSTALL 


LSTALL is a flag indicating whether a full listing of the 
contents of a file is to be given. For ex amp 1 e , in 
GDINFO , if LSTALL=YES , the grid identifiers will be listed 
if LSTALL=NO , grid identifiers will not be listed. 


3.46 MAP 


MAP is the map color, line type and line width separated 
by slashes. 

If the color is 0, the map is not drawn. If the map 
color, line type, or line width is blank, a default of 1 
is used. 


3.47 MAPFIL 


SMAPFIL is the name of the map file to be used for maps 
drawn by GEMPAK programs. If no directory is specified, 
the GEMPLT map files in GEMMAPS will be used. 

The map files in GEMMAPS are named by concatenating the 
resolution, map boundaries, and area with the three-letter 
source file type. For example, the medium-resolution 
political world map from GSFC is called MEPC7AD.GSF. 


RESOLUTION 

BOUNDARIES 

AREA 

SOURCE 

High 
MEd i um 
LOw 

Political 
COa s 1 1 i n e 
REg i ona 1 
CouN t y 

TOr 1 d 

NW quadrant 
NE quadrant 
SE quadrant 
SW quadrant 
WEst hemisphere 
North Hemisphere 
South Hemisphere 
US 

Ma r y 1 anD 

GSFc 

WI S c on s i n 

Th e default 

map is GEMMAPS 

:MEPOWO.GSF. 
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The following map files are currently available: 


HICNMD.GSF 
HICNUS.GSF 
HI PONE . GSF 
HIPONW.GSF 
HIPOSE.GSF 
HIPOSW. GSF 
HIPCMO.GSF 


MECCMD .WI S 
MEPOWE .WI S 
MEPOW.GSF 


LOCONH.GSF 
LOCOSH.GSF 
LOCONW.WI S 
LOREWE .WI S 


3.48 MARKER 


MARKER is the marker color, type, size, line width, and 
hardware/sof tware flag separated by slashes. 

If the marker color is 0, no markers will be drawn. If the 
marker color is not specified, a default of 1 will be used. 

The marker type specifies the shape of the marker to be 
drawn. If the type is unspecified or zero, the current 
marker type will be used. The software marker types are: 


1 

plus sign 

9 

Z wi th bar 

2 

hexagon 

10 

Y 

3 

t r i ang 1 e 

11 

box with diagonals 

4 

box 

12 

asterisk 

5 

sma 1 1 X 

13 

Xwith top, bott om bars 

6 

d i amond 

14 

star 

7 

up arrow 

15 

dot 

8 

X wi th top bar 

16 

large X 


Hardware marker types may differ. 

The marker size is a real number multiplier for the default 
marker size. If the size is zero or unspecified, the 
current size will be used. 

The hardware / software marker flag must be HW or SW. 
Otherwise, the current value is assumed. 
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3.49 MAXGRD 


MAXGRD is the maximum number of grids that can be stored in 
the grid file being created. 


3 . 50 MIXRLN 


MIXRLN specifies the color, line type, line width, minimum, 
maximum, and increment for the background water vapor mixing 
ratio lines. The values should be separated by slashes. 

If the color is 0, or MIXRLN is blank, no lines will be drawn. 


3.51 MRGDAT 


MRGDAT contains a flag indicating whether data is to be 
merged or unmerged. The type of merged data follows a 
slash with the default being 3. The valid types are: 

1 - mandatory data below 100 mb ( TTAA ) 

2 — ma ndatory and significant level data bel ow 100 mb 

( TTAA, TTBB , PPBB ) 

3 — mandatory and significant level data below and 

above 100 mb ( TTAA, TTBB, PPBB, TTCC, TTDD, PPDD ) 

If MRGDAT is used to create a sounding file, it indicates 
whether the file is to contain merged or unmerged data. 
Unme rged data sets will store the data as separate parts. 
Merged data sets contain values for each parameter at each 
level . 

If MRGDAT is used in a program to list data, the data will 
be merged if MRGDAT = YES. If MRGDAT - NO, the separate 
parts requested by the type will be listed. 
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3.52 NPASS 

NPASS controls the number of passes. Valid values are in the 
range 1-5. Note that two passes are STRONGLY RECOMMENDED. 

3.53 NTRACE 

NTRACE is the number of traces to be drawn in the meteogram. 

If NTRACE is less than 5, the user input in the traces after 
ntrace will be ignored. 

3.54 OLDGRD 

OLDGRD is the n ame of a GEMPAK3 grid file. 

3.55 OLDSFC 

OLDSFC is the name of a GEMPAK 3 surface file. 

3.56 OLDSND 

OLDSND is the name of a GEMPAK3 sounding file. 

3.57 OUTPUT 

OUTPUT determines the output devices. The valid devices are 
the t e rm inal, T, the default printer, P, a file, F or 
none , N . 
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Specify the output devices desired by entering the appropriate 
letters with no separators. 

If output is sent to a file, the file name may be specified 
after a slash. If no file name is entered, the output will 
be written to PROGRAM. FIL, where PROGRAM is the name of the 
program being executed. 

If N appears anywhere in the output string other than in the 
file name, no output will be generated. 


3.58 PANEL 


PANEL specifies the panel location, panel outline color, 
line type and width separated with slashes. 

The panel location determines the location of the view 
region on the graphics device. It ma y be specified using a 
number or abbreviation as follows: 


NUMBER 

ABBREVIATION 

DESCRIPTION 

0 

ALL 

Ent ire devi ce 

1 

UL 

Upper left quadrant 

2 

UR 

Upper right quadrant 

3 

LL 

Lower left quadrant 

4 

LR 

Lower right quadrant 

5 

L 

Left half 

6 

R 

Right half 

7 

T 

Top half 

8 

B 

B o 1 1 om half 


The view region may also be specified as four numbers 
separated with semicolons, giving the lower left and 
upper right corners in fractions of the graphics display 
area. The origin is in the lower left of the display. 

For example, the lower left quadrant can be given as: 

0 ; 0 ; . 5 ; . 5 

If the panel location is unspecified, the current location 
is unchanged. 

The panel outline color, line type and line width specify 
the values used to dr aw a box around the vi e w region. 

If the color is 0 or unspecified, no box is dr awn . 
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3.59 POINTS 


POINTS determines the grid points at which the data values 
will be plotted. The four numbers, start i, start y, 
increment x and increment y, should be separated with 
semicolons. If any number is not included, the value 1 will 
be used. Note that grid point (1,1) is in the 1 owe r left 
hand corner of the grid. 

Ex amp 1 e s : 

POINTS-1 -- A data value will be plotted at every 

grid point. 

POINTS-2 ; 2 ; 2 ; 2 -- The data value will be plotted at 

grid point (2,2) and at every second 
grid point in each direction. 


3.60 POSN 


POSN is the position number to be used to plot the data. 

If the position number is 0, the data will be plotted 
centered at the station. 

The other position numbers are: 

7 

1 3 

2 0 4 

5 6 

8 


3.61 PROJ 


PROJ is the map projection, projection angles, and margins 
separated by slashes. 

The following simple map projections may be specified: 

MER Me r c a t o r 
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NPS North Polar Stereographic 

SPS South Polar Stereographic 

LCC Northern Hemisphere Lambert Conic Conformal 

SCC Southern Hemisphere Lambert Conic Conformal 

CED Cylindrical Equidistant 

MCD Modified Cylindrical Equidistant 

UTM Universal Transverse Mercator 

NOR North Orthographic 

SOR South Orthographic 

For simple map projections, the lower left and upper right 
corners of the graphics area should be specified in GAREA. 

The following full map projections may also be specified: 


MER 

(CYL) 

Me rcator 


CED 

(CYL) 

Cy 1 indrical Equidistant 


MCD 

(CYL) 

Modified Cylindrical Equidistant 


STR 

(AZM) 

Polar Stereographic 


AED 

(AZM) 

Az imu tha 1 Equ idistant 


ORT 

(AZM) 

Orthographic 


LEA 

(AZM) 

Lambert equal area 


GNO 

(AZM) 

Gnomon i c 


LCC 

(CON) 

Northern Hemisphere Lambert Conic 

Conf o rma 1 

SCC 

(CON) 

Southern Hemisphere Lambert Conic 

Conf o rma 1 

UTM 

(OBM) 

Universal Transverse Mercator 


TVM 

(OBM) 

Transverse Me rcator 



For these projections, GAREA specifies the lower left and 
upper right corners. 

In addition, three angles MUST be specified in PROJ . The 
angles have the following meanings for the different 
projection classes: 


CYL 

AZM 


CON 


OBM 


The angles are not used. A default value for 
angle2, polon is used. 

anglel -- polat must be +90. or -90. 
angle2 -- polon is the central longitude 
ang le3 -- not used 

anglel -- standard latitude 1 

angle2 -- polon is the central longitude 

angle3 -- standard latitude 2 


anglel -- tangential longitude 
angle2 -- not used 
angle3 -- not used 


The angles for the full map projection types are given as 
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three numbers separated with semicolons. Note that THREE 
angles must be entered even if some angles are not used. 

There are two satellite projections available: 

AO I specifies AO I PS/2 navigation where GAREA is an 
AOIPS image name 

NPG specifies Naval Postgraduate School navigation 
where GAREA is the base name of NPGS navigation 
lat/lon files. 

If the projection is DEF , the current map projection will 
be used. 

In the grid programs the following graph projections are also 
available: 


LIN 

linear x , 

linear y 

LOG 

linear x , 

1 oga r i thmi c y 

KAP 

linear x , 

y ** KAPPA 

POL 

polar coordinates ( R, THETA ) 


Margins may be input as four numbers separated with semicolons 
NM wi 1 1 specify that no margins are to be set. 


3.62 PTYPE 


PTYPE is the type of y axis plot to be used, the height-to- 
width ratio of the plot, and the margins, separated by slashes 

The valid inputs for plot type for the y axis are: 


LIN 

1 inear 

LOG 

1 oga r i thmi c 

STUVE 

scaled by raising to KAP PA , 2 / 7 , p owe r 

KAP 

same as STUVE 

SKEWT 

logarithmic y-axis; skewed x-axis 

SKEW 

s ame a s SKEWT 


If the he i gh t - t o -wi d t h ratio is 0, the entire screen will 
be used. This is also the default when ratio is not specified 

The margins are specified as: left;bottom;right;top. The 

values are mu ltiples of the current character size. If the 
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margin specification is missing or incomplete, a default 
appropriate to the program will be used. The default in the 
profile programs is: 10;5;(6 times wind barb);5. 

The right margin allows room for plotting winds and 
typically has a value around 20. 


3.63 RAOBFL 


RAOBFL is the name of the raob data file created by the 
TRAFIC real-time ingest program. Data from this file will 
be extracted and written to a GEMPAK sounding file. 


3.64 REGION 


In GPBOX, REGION may be VIEW, PLOT or DEVICE where: 

VIEW corresponds to the view region (includes margins); 

PLOT corresponds to the plot region (excludes margins); 

DEVICE corresponds to the entire device space. 

In OABOX, the region must be GRID, DATA, or EXTEND where: 

GRID corresponds to the area over which the objective 
analysis is performed; 

DATA corresponds to the area from which data will be 
extracted for the objective analysis; 

EXTEND corresponds to the extended grid area over which 
the first pass of the objective analysis is 
pe r f o rmed . 


3.65 RESPOND 


$RESPOND is a logical variable indicating whether the user 
will respond interactively to GEMPAK programs. If set to 
NO, programs will not wait for input from the user. 
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RESPND does not appear as a parameter in any program. 

In order to change RESPND within the TAE , the global variable, 
SRESPOND must be modified. This can be done either in the 
GEMGLB tutor or with the TAE command LET in command mode. 

In no-TAE mode, simply set SRESPOND to YES or NO in any 
p r o g r am . 


3.66 SAVF I L 


SAVFIL is the name of a file in which graphics from a 
graphics device can be saved. The graphics can only be 
saved on a device with readback capability. This file does 
not save graphics commands as they are executed. Rather, it 
saves a bitmap of the current graphics screen. GPSAVE can 
be used to create the file. GPREST wi 1 1 display a file 
stored using GPSAVE. The only device which currently has 
this capability is the IIS. 


3.67 SCALE 


SCALE is the scaling factor for the data. All data will be 
multiplied by 10 ** SCALE. 

If SCALE is not specified, an appropriate scaling will be 
selected in the grid programs. In the surface and sounding 
programs, a default of 0 will be used. 

If the absolute value of SCALE is greater than 5 in a 
surface or sounding program, a value of 0 will be used. If 
the absolute value in a grid progr am is greater than 20, 
a default will be computed. 


Note that scaling data may create integer overflow problems 
in programs, such as SEMAP, which round data to the nearest 
integer. Valid integers must be in the range -2*10**9 to 
2*10**9, approximately. 
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3.6 8 SEARCH 


SEARCH is used to control the search radius in an objective 
analysis program. The radius is the maximum distance that 
a station may be from a grid point to be used in the analysis 
for that point. The search radius will be set so that stations 
whose weighting factor would be less than EXP (-SEARCH) will 
not be used, SEARCH must have a value in the range 1 - 50. 

If the value is outside this range, a default value of 20 will 
be used. If SEARCH is given a very small value, many grid 
points will not have three stations within the search area, and 
consequently will be set to the missing data value. 


3.69 SFEFIL 


SFEFIL is the name of the surface edit file which will be 
used to update a surface file. 

SFEFIL is a text file which may be created using the program 
SFLIST with F as an output device. A text editor may be 
used to create or change the surface edit file. 

The parameters to be edited must be specified at the 
beginning of the surface edit file. For example: 

PARM = TMPF ; EMTF 

Parameters that have character values, such as WTHR , cannot 
be edited. Instead, the numeric value must be used, e.g., 
W604 in the case of WTHR. All the parameters which are to be 
edited mu st already exist in the surface data file, since no 
parameter transformations will be done. Station information, 
such as latitude, longitude and elevation, cannot be changed 
with SFEFIL. Use GEMPAK program SFSTNS to update station 
information in a surface file. 

Th e data foil ow the par ame ter list. 

Stations and times which do not already exist in the surface 
data file will be added. 


3-37 



GEMPAK VARIABLES 


3.70 SFFILE 


SFFILE is the name of the surface data file to be accessed. 


3.71 SFOUTF 


SFOUTF is the output surface data file. 

SFOUTF is used in programs which create a new surface data 
file. It is also used in programs in which data is moved 
from an input file to an output file, in which case, SFFILE 
is the original file. 


3.72 SFPARM 


SFPARM is a list of surface parameters to be used in a 
surface program. 

The four-character parameter names must be separated by 
semicolons. For example: 

SFPARM - TMPF ; EftVPF ;WTHR 

If a parameter is not present in a data set, the program 
will attempt to derive it from those parameters which are 
present. For example, relative humidity can be computed from 
the temperature and dewpoint temperature. Note that STHA, 
STHE , and SMXR are computed from PALT. 

Functions can be used to mod i f y /qua 1 i f y the output 
parameters. These functions must follow the parameter name 
and are applied in the order given. The function symbols 
must precede the qualifying value and are defined as follows: 


* 

mu 1 t i p 1 y 

- 

equal to 

/ 

divide 

$ 

direction 

+ 

add 

% 

vertical coordinate 

- 

subtract 

i 

depth 

< 

less than 



> 

greater than 
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The GEMPAK surface parameters and the corresponding four- 
character abbrevations are: 


TMPC 

TEMP 

TMPF 

TMPK 

STIIA 

STHE 

TVRK 

TVRC 

TVRF 

THTV 

TMAX 

TMIN 

SSTC 

LTMP 


DWPC 

DWPT 

DWPF 

LWPK 

DPDC 

DPDF 

DPDK 

SMXR 

SMXS 

RELH 

VAPR 

VAPS 


PMSL 
PALT 
ALT I 
ALTM 
SALT 
SMSL 
SAL I 
ZMSL 
Z000 
Z950 
Z900 
Z850 
Z800 


TEMPERATURE 


temperature in Celsius 

temperature in Celsius (alternate name) 
temperature in Fahrenheit 
temperature in Kelvin 

surface potential temperature in Kelvin 

surface equivalent potential temperature in Kelvin 

virtual temperature in Kelvin 

virtual temperature in Celsius 

virtual temperature in Fahrenheit 

virtual potential temperature in Kelvin 

maximum 24 -hour temperature in Celsius 

minimum 24-hour temperature in Celsius 

sea surface temperature in Celsius 

surface lifted temperature 


MOISTURE 


Celsius 

Celsius (alternate 

Fahrenheit 

Kelvin 

Celsius 
Fa r enhe i t 
Kelvin 
in g/kg 


d ewpo i n t 

d ewp oint in Celsius (alternate n ame ) 
d ewp oint 
d ewp oint 
d ewp oint 
d ewp oint 
d ewp oint 
surface 

surface saturated mixing ratio in g/kg 
relative humidity in percent 
vapor pressure in millibars 
saturation vapor pressure in millibars 


1 n 
i n 
i n 
i n 

depression 
depression 
depression 
mixing ratio 


PRESSURE AND ALTIMETER 


mean sea level pressure in millibars 
pressure in millibars from altimeter 
altimeter setting in inches of mercury 
altimeter setting converted to millibars 
abbreviated altimeter in millibars 
abbreviated mean sea level pressure 
abbreviated ALT I 
estimated height at PMSL 
estimated height at 1000 mbs 
estimated height at 950 mbs 
estimated height at 900 mbs 
estimated height at 850 mbs 
estimated height at 800 mbs 
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WIND 

UWND - u component of the wind in me t e r s / s e cond 

VWND - v component of the wind in meters/second 

UKNT - u component of the wind in knots 

VKNT - v component of the wind in knots 

DRCT - wind direction in degrees 

SPED - wind speed in meters/second 

SKNT - wind speed in knots 

PSPD - packed speed and direction (ddfff) in me t e r s / s e c ond 

PKNT - packed speed and direction (ddfff) in knots 

GUST - wind gusts in knots 

GUMS - wind gusts in me t e r s / s e c ond 

CLOUD 

For the following parameters, the character x may be 
replaced by L, M, or H, indicating low, middle or high 
clouds. The character T indicates the value of the parameter 
at the level of maximum cloud coverage. 

xCLD - character cloud coverage 

TCLD - character maximum cloud coverage 

xCLO - fractional cloud coverage 

TCLO - fractional maximum cloud coverage 

CLCx - numeric cloud coverage 

CLCT - numeric maximum cloud cover 

CLDS - combined cloud coverage from three levels 

CMBC - numeric combined cloud coverage from three levels 

CLHx - cloud height in hundreds of feet 

CLDx - combined cloud height and coverage 

CLDT - maximum value of CLDx 

CQMx - numeric combined cloud height and coverage 
COMT - maximum value of COMx 

For the next two cloud parameters the x can be replaced by 
1, 2, or 3 and represents the cloud report number. 

CHCx - numeric combined cloud height and coverage 
CHDx - combined cloud height and short code 

The following numeric cloud parameters are WMO standard 
codes. 

CFRC - fraction of celestial dome covered by cloud fromWMO 
Code 2700 

CLFR - fraction of celestial dome covered by low or mid 
clouds from WO Code 2700 
CTYL - low-level cloud genera from WMO Code 0513 
CTYM - mid-level cloud genera fromWMO Code 0513 
CTYH - high-level cloud genera from WMO Code 0513 
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CBAS 


WTHR 

VCOD 

WNUM 

W604 

WSYM 

WTMO 

WRMO 

WASH 

RWSH 

PWTH 

PVMM 


STID 

STNM 

SLAT 

SLON 

SELV 


VS BY 

VSBK 

P06 I 

P06M 

P24 I 

P24M 

SNOW 

HAIL 

HLSZ 

WHGT 

WPER 


SPAC 

BLNK 


Cloud base height from MO Code 1000 
WEATHER 

character weather code from any numeric code 
character weather code 
numeric weather code 

numeric weather code from the 604 file 
weather symbol number 
character we ather code fr om VMO 
numeric VMO weather code 
Asheville surface data weather code 
numeric Asheville surface data weather code 
character past we ather in VMO code 
numeric past weather in WMO code 

STATION VARIABLES 

character station identifier 
5-digit station identifier 
station latitude in degrees 

station longitude in degrees. West is negative 
station elevation in meters 

MI SCELLANEOUS 


visibility in 
visibility in 
precipitation 
precipitation 
precipitation 
precipitation 
s n ow depth in 
hail flag 

hail size in centimeters 
wa ve height in me t e r s 
wave period in seconds 


statute miles 
k i 1 ome t e r s 
over last 6 hours 
last 
last 
last 


over 
over 
over 
inches 


in hundreths of inches 
6 hours in millimeters 
24 hours in hundreths of inches 
hours in millimeters 


24 


SPACING 

plot a space, accounted for in FILTER 
plot a blank, not accounted for in FILTER 


3.73 SFPRMF 


SFPRMF is the name of the surface parameter packing file or 
the info rma tion itself entered as foil ows : 
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SFPRMF-PRM1 /MINI -MAXI -RES1 ; PRM2 /MIN2 -MAX2 -RES2 ; ... 

where PRMn is the parameter name, MINn is the minimum for 
PRMn , MAXn is the maximum for PRMn, and RESn is the 
resolution. These entries are the same as those placed in 
a packing file as described below. If the data are not to 
be packed, the list of parameter names separated by semicolons 
can be entered: SFPRMF = PRM1 ; PRM2 ; . . . 

IF SFPRMF specifies a file name, then each line of the file 
must contain the following information separated by blanks 
or t abs : 


parameter name CHAR *4 

minimum data value REAL 

maximum data value REAL 

resolution REAL 


The resolution should be an integral power of 10; otherwise 
the next smaller resolution will be used. For example, 
res - .5 will become .1. 

If the data are not to be packed, the minimum and maximum 
data values and the resolution should not be included. Note 
that either all of the parameters or none of them must have 
packing info rma tion. 

GEMTABL : SFPACK . TBL is a typical table used to pack surface 
data. This file may be used as an example when creating a 
packing file. 


3.74 SHIPFL 


SHIPFL is a logical parameter which indicates that the 
surface file contains stations wh ich are not at a fixed 
location, such as moving ships, aircraft, or floating buoys. 

A ship file will store and write data differently from a 
standard GEMPAK surface file. However, these files can be 
read using the surface library and all programs which access 
data in surface files can be used unchanged. 

Note that station information cannot be added to these 
files, since the station information must be stored with the 
data. 
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3.75 SKPMIS 


SKPMIS is a logical flag which indicates whether non- 
reporting stations will be listed. 

If SKPMIS is YES, nonreporting stations will not be 
listed. Otherwise, all stations will be listed. 

In SFMAP , markers can be plotted at nonreporting 
stations when SKPMIS is NO. 


3.76 SNEFIL 


SNEFIL is the name of the sounding edit file which will be 
used to update a sounding file. 

SNEFIL is a text file which may be created using the program 
SNLISTwith F as an output device. A text editor may be 
used to create or change the edit file. 

Th e par ame ters to be edited mu st be specified at the 
beginning of the edit file* For example: 

SNPARM - PRES ;TMPF ;DvVPF 

If the output sounding file exists, the data set parameters 
must be the same as those listed in the edit file. Station 
information, such as latitude, longitude, and elevation, 
cannot be changed with SNEFIL. Use GEMPAK program SNSTNS 
to update station information in a sounding file. 

If the file does not exist, a new, unpacked file with these 
parameters will be created. 

The data follows the parameter list. Only level data will 
be added to the file. Stability indices from STNDEX will be 
i gno r ed . 

Stations and times which do not already exist in the surface 
data file will be added. 
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3.77 SNFILE 


SNFILE is the filename for an upper air data set. 


3.78 SNOUTF 


SNOUTF is the output sounding data file. 

SNOUTF is used in programs which create a new sounding data 
file. It is also used in programs in which data is moved 
from an input file to an output file, in which case SNFILE 
is the original file. 


3.79 SNPARM 


SNPARM is a list of upper air parameters to be used in 
an upper air program. 

The four-character parameter names must be separated by 
semicolons. For example: 

SNPARM - PRES ; TMPC; DWPC;THTA 

If a par ame ter is not present in a data set, the progr am 
wi 11 a 1 1 emp t to derive it fr om those par ame t e r s wh i c h are 
present. For example, if pressure, temperature and 
dewpoint are present in a data set, then relative humidity 
canbederived. 

The layer parameters are computed over a default layer between 
the specified level and the next level upward in the sounding. 
Th e depth of the layer for layer averages ma y be specified 
preceded by ! in the user input. The shear and RICH 
calculations can be done relative to a direction by specifing 
the direction preceded by $ in the user input. 

In-line functions can be use to modify/qualify the output 
parameters. These functions are applied in the order given. 

Th e function s ymb o 1 s mu st precede the qualifying value and 
are defined as foil ows : 
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* multiply - equal to 

/ divide $ direction 

+ a( id % verticalcoordinate 

subtract I depth 

< less than 

> greater than 

The following page lists the GEMPAK sounding parameters and the 
corresponding four-character abbreva t i ons . 

TEMPERATURES 
temperature in Celsius 

temperature in Celsius (alternate name) 
temperature in Fahrenheit 
temperature in Kelvin 
potential temperature in Kelvin 
virtual potential temperature in Kelvin 
virtual temperature in Kelvin 
virtual temperature in Celsius 
virtual temperature in Fahrenheit 

MOISTURE PARAMETERS 

IWPC - d ewp oint in Celsius 

IWPT - dewpoint in Celsius (alternate name) 

EWPF - dewpoint in Fahrenheit 

EWPK - d ewp oint in Kelvin 

DPDC - dewpoint depression Celsius 

DPDF - dewpoint depression Farenheit 

DPDK - dewpoint depression Kelvin 

MIXR - mixing ratio in g/kg 

MIXS - saturation mixing ratio in g/kg 

PWTR - precipitable water in mm 

RELH - relative humidity in percent 

THTE - equivalent potential temperature in Kelvin 
VAPR - vapor pressure in millibars 
VAPS - saturation vapor pressure in millibars 
LHVP - latent heat of vaporization in J/kg 

HEIGHT 

HGHT - height in meters 

HGTM - height in meters (alternate name) 

HGTK - height in kil ome t e r s 

HGTD - height in dek ame t e r s 

HGFT - height in feet 

HGFH - height in hundreds of feet 

HGFK - height in thousands of feet 


TMPC - 
TEMP - 
TMPF - 
TMPK - 
THTA - 
THTV - 
TVRK - 
TVRC - 
TVRF - 
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HGML 

STDZ 

DHGT 

MHGT 


PRES 
PMSL 
ALT I 


UWND 

VWND 

UKNT 

VKNT 

DRCT 

SPED 

SKNT 

PSPD 

PKNT 

VOIP 

WNML 


LCLT 

LCLP 


Not e : 

SHOW 
LIFT 
SWET 
KINX 
CTOT 
VTOT 
TOTL 
CAPE - 
CINS - 
EQLV - 
LFCV - 
BRCH - 
MLTH - 
MLMR - 


- height in miles 

- standard height convention used on upper air charts 

- dry hydrostatic height in meters 

- moist hydrostatic height in meters 

PRESSURE AND ALTIMETER 

- station pressure in millibars 

- mean sea level pressure in millibars 

- altimeter setting in inches from PRES 

WINDS 

- u component of the wind in meters/second 

- v component of the wind in meters/second 

- u component of the wind in knots 

- v component of the wind in knots 

- wind direction in degrees 

- wind speed in me t e r s / s e cond 

- wind speed in knots 

- packed speed and direction (ddfff) in me t e r s / s e cond 

- packed speed and direction (ddfff) in knots 

wind component toward a specific direction 

wind component toward 90 deg from a specific direction 

LIFTED CONDENSATION LEVEL 

temperature in Kelvin of parcel raised to LCL 
pressure in millibars of parcel raised to LCL 

STABILITY INDICES 

Default depths are enclosed in parentheses. 

Showa Iter Index 
Lifted Index (100 mb) 

Sweat Index 
K Index 

Cross Totals Index 
Vertical Totals Index 
Total Totals Index 

Convective Available Potential Energy (500 m) 
Convective Inhibition (500m) 

Equilibrium Level 

Level of Free Convection 

Bulk Richardson Number (500 m, 6000 m) 

Mixed layer mean potential temperature (500 m) 

Mixed layer mean mixing ratio (500 m) 


LAYER QUANTITIES 
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Note: 


RICH 

BVFQ 

BVSQ 

BVPD 

LAPS 

STAB 

STAP 

SHRM 

SHRD 

SEPA 


STID 

STNM 

SLAT 

SLON 

SELV 

RANG 

AZIM 

LATI 

LONG 

DELT 


DDEN 

PSYM 


SPAC 

BLNK 


Default depths are the layer depths in the sounding 
except for SEPA for which the layer depth is 5 K. 

Richardson number 
Brunt-Vai sala frequency 
Brunt -Va i sa 1 a frequency squared 
Brunt -Va i sa 1 a frequency 
Temperature lapse rate 
Potential temperature lapse rate 
Stability with respect to pressure in k/mb 
Wind shear magnitude 
Wind shear direction 

Vertical separation in mb between isentropes 

STATION VARIABLES 

character station identifier 
5-digit station identifier 
station latitude in degrees 

station longitude in degrees, We st is negative 
station elevation in meters 
range (specialized data) 
azimuth (specialized data) 
latitude in degrees 

longitude in degrees, West is negative 
delta time (specialized data) 

MI SCELLANEOUS 

density of dry air in kg/(m**3) 

Montgomery stream function 

SPACING 

plot a space, accounted for in FILTER 
plot a blank, not accounted for in FILTER 


3.80 SNPRMF 


SNPRMF is the name of the sounding parameter packing file or 
the information itself entered as follows: 

PRM1 /MINI -MAXI -RES1 ; PRM2 /MIN2 -MAX2 -RES2 ; . . . 

where PRMn is the parameter name, MINn is the minimum for 
PRMn , MAXn is the maximum for PRMn, and RESn is the 
resolution. These entries are the same as those placed in a 
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packing file as described below. If the data are not to be 
packed, the list of parameter names separated by semicolons 
can be entered: SNPRMF = PRM1 ; PRM2 ; . . . 

IF SNPRMF specifies a file name, then each line of the file 
must contain the following information separated by blanks 
or tabs : 


parameter name CHAR *4 

minimum data value REAL 

maximum data value REAL 

resolution REAL 


The resolution should be an integral power of 10; otherwise 
the next smaller resolution will be used ( e.g., res = .5 
wi 11 bee ome . 1 ) . 

If the data are not to be packed, the minimum and maximum 
data values and the resolution should not be included. Note 
that either all of the parameters or none of them must have 
packing information. 

GEMTABL: SNPACK.TBL is a typical table used to pack surface 
data. This file may be used as an example when creating a 
packing file. 

Note that if MRGDAT is NO, the data will be packed using the 
standard packing for unmerged data, and the value of SNPRMF 
will be i gnor ed . 


3.81 SOURCE 


SOURCE indicates whether the data used to compute the 
average station spacing are to be read fr om a surface or 
sounding file. Use SF to read from a surface file and SN 
to read from a sounding file. 


3.82 STATION 


STATION is the station to use in the meteogram program. 
Either a station identifier or station number may be 
entered. 
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If a single station is entered, all traces will use data 
from that station. If a list of stations is entered, 
trace 1 will be for station 1, trace 2 for station 2, .. 

Stations in the list must be separated using semicolons. 


3.83 STNCOL 


STNCOL specifies the color for the station identifier, 
time and the parameters specified in STNDEX . These 
parameters are written at the top of the plot. 


3.84 STNDEX 


STNDEX is the list of stability indices or station parameters 
for upper air data. The items in the list must be separated 
using semicolons. The depth for the layer averages may 
be specified preceded by a ! in the user input. For example, 
BRCH1 1000 1 8000 instructs the program to average over a mixed 
layer 1000 meters deep and lower tropospheric layer 8000 
meters deep. Similarly MLTH1750 results in an average 
potential temperature over a mixed layer 750 meters deep 
while the default is 500 meters. 


foil owi ng 

is a list of valid names: 

CTOT 

Cross totals index 

KINX 

K index 

LIFT 

Lifted i nd ex 

SHCW 

Showalter index 

SWET 

SWEAT index 

TOTL 

Total totals index 

VTOT 

Vertical totals index 

CAPE 

Convective Available Potential Energy 

CINS 

Convective Inhibition 

EQLV 

Equilibrium Level 

LFCV 

Level of Free Convection 

BRCH 

Bulk Richardson Number 

LCLT 

Temperature in Kelvin of the LCL 

LCLP 

Pressure in millibars of the LCL 

MLTH 

Mean mixed layer potential temperature 

MLMR 

Mean mixed layer mixing ratio 

STID 

Station identifier 

STNM 

Station n umb e r 
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SLAT Station latitude 

SLON Station longitude 

SELV Station elevation 


3.85 STNFIL 


STNFIL is the name of a station file which contains station 
information. This includes the character identifier, number, 
name, state, country, latitude, longitude and elevation for 
each station. 

All this information, except the station name, is stored as 
station header information in surface and sounding data 
files. 

The file GEMTABL : SFSTNS . TBL is a surface station table 
containing US and Canadian stations. GEMTABL : SNSTNS . TBL is 
an upper air table for the Northwest quadrant. 

Info rma tion in the files mu st be stored using the exact 
format used in current tables, since they are read with a 
FORTRAN FORMAT statement. They may be changed using a text 
editor. 


3.86 TAXIS 


TAXIS contains the range, increment and location for labels 
on a time axis input as: 

START- STOP- INC ; LB ; GL ; TM 


START and STOP are GEMPAK date/ 1 ime s wh i ch ma y be abbreviated. 
If the values are omitted, the data range will be used. The 
time axis will be reversed if the TAXIS specification begins 
with R or if the times are input with a later time first. 

INC is the time increment in hours and minutes. The form 
for INC is HHHMM. If one or two digits are entered, hours 
will be assumed. If INC is omitted, a default appropriate 
for the range will be used. If the time range exceeds 
720 days, the increment is ignored, and an appropriate 
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labeling interval is selected automatically. 

LB, GL, and TM are the frequencies for labels, grid lines 
and tick ma rks . 

Ex amp 1 e s : 

5/12-6/18-6 Draw the axis from 5/12 to 6/18 with 

labels at 6-hour increments. 


3.87 TEXT 


TEXT is the size, font, text width and ha rdwa r e / s of twa r e 
flag for graphics text separated with slashes. 

The size is a real number multiplier for the default text 
size. If the size is zero or unspecified, the current size 
wi 11 be used. 

The text width is the integer line width to be used in 
generating software text. If the text size, font or width 
is not specified, the current value is used. 

Th e h a r dwa re/ s o f twa re selector mu s t be MV o r SW t o 
change to hardware- or software-generated text. 

Ex amp 1 e s : 

TEXT - 1/2/ /HW -- text size - 1, hardware text font 2 

TEXT = 2.5 -- text size - 2.5, current text font 


3.88 THTALN 


THTALN specifies the color, line type, line width, minimum, 
maximum, and increment for the background dry adiabats 
(potential t emp erature lines). Th e values should be 
separated by slashes. If THTALN is blank, no lines will 
be dr awn . 
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3.89 THTELN 


THTELN specifies the color, line type, line width, minimum, 
maximum, and increment for the background moist adiabats 
(equivalent potential temperature lines). The values 
must be separated by slashes. If THTELN is blank, no lines 
will be drawn . 


3.90 TIMSTN 


TIMSTN contains the maximum number of times to include in a file 
and the number of stations to be included in addition to 
the stations in STNFIL . 


3.91 TITLE 


TITLE is the title color, title line, and title string 
separated by slashes. 

If the title color is 0, a title is not plotted. 

The title line specifies the line on which the title will 
be written. If the line is not specified, a default of 0 

will be used. The value of the title line has the following 
meaning s : 6 

0 b o 1 1 om line 
* n n lines fr om b o 1 1 om 

+n n lines from top 

The title string is the title to be written. If no title 
string is specified, a default title will be dete rm i n e d by 
the progr am . 
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3.92 TRACE 


The TRACE parameters contain specifications for each trace 
on the me t e o g r am in the foil ow i n g f o rma t : 

par ame ters/colors/range /wi t ne s ! pa r ame ters/ colors/range /wi t ne s 

The parameters before the 1 will be plotted on the left of the 
plot; those after the 1 will be plotted on the right. 

The parameters may be any GEMPAK surface parameter. Real 
valued parameters will be drawn as a graph. Character valued 
parameters will be rotated 90 degrees and written on the plot. 
In addition to the standard surface parameters, the following 
are also valid: 


BARB 

ARRW 

DARR 

WSYM 


wi nd barbs 
wind arrows 
directional arr ows 
we a t h e r s ymb o 1 s 


For GUST and GUMS, the character G wi 1 1 be plotted. 

Up to four parameters may be plotted along each axis. The 
parameters must be separated using semicolons. Character 
and weather symbol data may only be plotted in positions 
1 , 2 or 3 . 

Each parameter name may be followed by an asterisk, the 
size or line type, a second asterisk and the width. For 
example, WSYM*. 5*5 will draw weather symbols half the 
default size with a line width of 5. TMPF*3 will plot 
a t emp erature line using dash pattern 3. 

The colors for the parameters must also be separated using 
semicolons. If a single number is entered, all parameters 
are dr awn in that color. If a zero is entered, the current 
default color is used. 


The range specifies the scaling of the y axis. The format 
is: start ; stop; increment , Note that in this program, the 

parts of range must be separated using semicolons. If no 
range is given, it is selected using the data values. 

Witness lines are specified in witnes. These are horizontal 
dotted lines. A list of y values may be entered separated 
by semicolons. If the value of witnes is YES, a witness line 
will be centered on the plot. 
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3.93 TRACE 1 


See the TRACE documentation. 


3.94 TRACE 2 


See the TRACE documentation. 

3.95 TRACE 3 

See the TRACE documentation. 

3.96 TRACE 4 

See the TRACE documentation. 

3.97 TRACE 5 

See the TRACE doc ume n t a t i on . 

3.98 VCOORD 


VCOORD specifies the vertical coordinate system of the 
levels to process. There are four coordinates: 

NONE — surface data only 

PRES - pressure 

THTA - theta (isentropic) 

HGHT - height 
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A list or range of levels will be based on the vertical 
coordinate type, with the exception of MAN or VAS which 
a 1 wa ys refer to a pressure coordinate syst cm . 

Note that data are ignored within supe r ad i aba t i c layers 
found wo r k i n g u pwa r d s f r om the surface wh en interpolations 
are performed to isentropic coordinates . 


3.99 WIND 


WIND specifies the wind symbol, size, width, type, and 
head size separated by slashes. 

The wind symbol contains a letter for symbol type, a letter 
for symbol units and a color number with no separators. 

The character meanings are: 

TYPE B - BARB A - ARROW 

UNITS K « KNOTS M = Meters/sec 

COLOR Color number 0 » no wind plotted 

The default is BM1 , i.e., barbs in meters/sec plotted in 

color n umb er 1. If a partial specification is given, the 
remaining characteristics will be taken from the default. 

The wind size is a real number which will be used as a 
multiplier for the default wind symbol size. If this number 
is negative, zero, or missing, the current size will be used. 
The sizes for barbs and arrows are independent. 

The wind type defines the arrow or barb type to be used. 

If the wind type is equal to 1, calm winds are plotted. If 
the wind type is 2, calm winds will not be plotted. 

The head size is a real valued multiplier used for the length 
of the arrow head. This variable is not used for wind barbs. 

Note that the units will be ignored in GEWIND. In that 
program, the units are determined by the vector requested in 
GVECT . 
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3.100 WINPOS 


WINPOS specifies the position for plotting winds. Up to three 
separate wind profiles may be plotted. The stability indices 
will also be positioned at the top of the plot according to 
WINPOS. Position 1 is the leftmost position. The value of 
WINPOS will be incremented modulo 3 until the screen is 
cleared or the program is run again. 


3.101 XAXIS 


XAXIS contains the left bound, right bound, label increment, 
and frequency information separated by slashes in the form: 

lef t/right/incr eme nt/lbfq;gdfq; tkfq 

The frequencies for labels, grid lines, and tick marks 
follow the last slash and are separated with semicolons. 
Appropriate defaults are provided if values are not 
specified. If the increment is positive, all the label 
values will be divisible by the increment. If it is 
negative, the label values will begin with the left value 
and be separated by the increment. The bounds and 
increment are for the scaled data. 


3.102 YAXIS 


YAXIS contains the lower bound, upper bound, label increment, 
and frequency information separated by slashes in the form: 

1 owe r /uppe r / inc r emen t / 1 bf q ;gdfq;tkfq 

The frequencies for labels, grid lines, and tick marks 
follow the last slash and are separated with semicolons. 

If the incr eme nt is positive, all the label values will 
be divisible by the increment. If it is negative, 
the label values will begin with the lower value and 
be separated by the increment. The bounds and 
increment are for the scaled data. 

Defaults will be supplied if no specification is given. 
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The following defaults will be used when appropriate 


vcoord 

1 owe r 

iipp e r 

defaults 

PRES 

1020 

100 

Mandatory levels 

THTA 

270 

400 

10 

HGHT 

0 

20000 

1000 
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CHAPTER 4 
GEMPAK PROGRAMS 

This chapter contains documentation on all the GEMPAK programs. 
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4 . 1 PROGRAM GDCF I L 


GDCF IL creates a GEMPAK grid file. 
INPUT PARAMETERS 


GDOUTF 

PROJ 

GRDAREA 

KXKY 

MAXGRD 

CPYF I L 

ANLYSS 


♦Output grid file 

*Ma p p r o j e c t i on/ ang 1 e s /ma r g i n s 

♦Area covered by grid 

* Number of grid points in x;y 

♦Maximum number of grids 

♦File to be copied 

♦Grid analysis block 


PROGRAM DESCRIPTION 

This program creates a GEMPAK grid file. Information about 
the file to be created is received fr om the user, fr om the 
grid navigation table, or from another grid file. 

Each grid file must contain navigation information. It 
must also contain analysis block information which is used 
wh en an objective analysis progr am wr ites to the grid file. 

If the value in CPYFIL begins with a #, the rest of the 
string specifies a grid number or name to be found in the 
grid navigation table. This table currently contains the 
following grid names and numbers: 

LFM 26 

NGMC 101 
SFUS 301 
SFWV 302 
SNUS 311 
SNWV 312 

Any other non-blank entry in CPYFIL is taken as the name of an 
existing grid file, whose navigation and analysis blocks will 
be copied to the new grid file. 

If CPYFIL is blank, the navigation and analysis information 
will be obtained from PROJ, GRDAREA, KXKY, and ANLYSS. KXKY 
must contain the number of points in the x and y directions. 
The two integers must be separated with a semicolon. If 
the grid projection type specified in PROJ is CED, the 
grid spacing consisting of deltax and deltay may be entered 
in KXKY by prefixing the numbers with a #. Deltax and deltay 
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are in degrees longitude and degrees latitude, respectively. 
The analysis information must contain deltan, which is the 
station spacing, and the extend area separated by a slash. 
The extend area is four integers separated by semicolons. I 
ANLYSS is blank, default values will be assigned. If there 
is no analysis block in the file to be copied, ANLYSS is 
used to create one. 

MAXGRD is the maximum number of grids which will be allowed 
in the grid file. 


EXAMPLES : 

1. Create a grid file called LFM.GRD from the information in 
the grid navigation table, allowing up to 100 grids in the 
file. PROJ , GRDAREA and KXKY are ignored since this 
information will be obtained from the grid navigation 
table. 

GDOUTF - lfm.grd 
PROJ 

GRDAREA - 
KXKY 

MAXGRD - 100 
CPYF I L - # 1 f m 
ANLYSS - 

2. Create a grid file called SOUND . GRD wi t h a maximum of 15 
grids. The grid dimensions are 20 by 30; the area 

is US with a Mercator projection. The station 
spacing is 5 degrees of latitude and the grid extension 
for objective analysis is 3 grid points. 

GDOUTF = sound.grd 
PROJ - me r 
GRDAREA = us 
KXKY = 20; 30 
MAXGRD =15 
CPYF I L = 

ANLYSS = 5/3 ; 3 ; 3 ; 3 

3. Create a file for evenly spaced lat/lon grids over area IL- 
using a grid spacing of .75 degrees in each direction. 
Compute reasonable values for the analysis block. 

GDOUTF = i 1 . g r d 
PROJ = ced 
GRDAREA = i 1 - 
KXKY = #.75; .75 
MAXGRD = 100 


4-3 



GEMPAK PROGRAMS 


CPYFIL - 
ANLYSS - 

4. Create a grid file called NEW.GRD from the information 
within grid file OLD . GRD . 

GDOUTF - new.grd 
PROJ 

GRDAREA - 
KXKY 

MAXGRD - 100 
CPYFIL - old.grd 
ANLYSS - 


ERROR MESSAGES 

[GDCFIL -l] Fatal error initializing TAE. 

[GDCFIL -2] Fatal error reading TAE parameters. 

[GDCFIL -3] Fatal error initializing GEMPLT. 

[GDCFIL -4] Invalid navigation. 

[GDCFIL -5] Grid area ... is invalid. 

[GDCFIL -6] Invalid grid size. 

[GDCFIL -7] The grid file name may not be blank. 

[ GDCF IL -8] Error reading table. 

[GDCFIL -9] Grid name ... cannot be found in grid table. 

[GDCFIL -10] Invalid extend region. Try 0;0;0;0. 
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4.2 PROGRAM GDCNTR 


GDCNTR d r aws contour lines through a gridded scalar field. 
INPUT PARAMETERS 


GDATTIM 

♦Grid date/time 

GLEVEL 

♦Grid level 

GVCORD 

♦Grid vertical coordinate 

GFUNC 

♦Gr id function 

GDFILE 

♦Grid file 

PROJ 

♦Ma p projection/angles/ ma r g i n s 

GAREA 

♦Graphics area 

LINE 

♦Line color/dash /wi dth/label 

CINT 

♦Contour interval /mi n /max 

SCALE 

♦Scaling factor 

MAP 

♦Ma p color/dash /wi d t h 

LATLON 

♦Line color/da sh/wi dth/label/inc 

TITLE 

♦Title color/line/title 

DEVICE 

♦Graphi cs devi ce 

PANEL 

♦Panel location/color /da sh/wi d th 

CLEAR 

♦Clear screen flag 

TEXT 

♦Text size/font / h w flag 


PROGRAM DESCRIPTION 

GDCNTR draws contour lines through a scalar grid computed 
using the GEMPAK grid diagnostic functions. 

The attributes of the contour lines, including the color, 
line type, line width, and label frequency are specified 
in LINE. The four attributes are separated with slashes; 
semicolons separate the values for each attribute. If the 
line type is set to a single negative number, negative 
contour values will have the absolute value of the line type 
and positive values will be solid. If the label type is set 
to a single number, n, then every nth value will be labeled. 

Note that color components may be set by appending « NAME to 
a color number. For example, 1=PURPLE will set color number 
1 to purple. 

Using the variable, Cl NT, the contour interval, minimum value, 
and maximum value may be entered. A scaling factor may be 
entered in SCALE. The data in the grid file will be multiplied 
by 10 ** SCALE. If no contour interval is entered, a default 
value which will generate 5 to 10 contour lines will be 
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selected. The scaling will be done before the contour interval 
is computed. 

A map, lat/lon grid, and title may also be included. 


EXAMPLES 

1. Draw contours every 4 degrees through the 850-mb 

temperature field for the 24 -h forecast time. Negative 
values will be dashed using line type 3; every other 
contour line will be labeled; all the contour lines will be 
drawn in color number 1. The graphics area will be 
centered on Colorado. The contours are drawn on a clear 
screen with a dotted map, dotted lat/lon lines every 5 
degrees and no title. 


GDATTIM 

- 

f 24 

GLEVEL 

- 

850 

GVCORD 

- 

pres 

GFUNC 

- 

t emp 

GDFILE 

- 

n gm . g r d 

PROJ 

- 

me r 

GAREA 

- 

c 0 

LINE 

* 

1/ -3/1/2 

CINT 

= 

4 

SCALE 



MAP 

- 

1/10 

LATLON 

= 

l/10/l/l/5;5 

TITLE 

act 

0 

DEVICE 

- 

vt 

PANEL 


0 

CLEAR 

« 

ye s 

TEXT 

** 



2. Using the values of the variables supplied above, the 
dewpoint lines are to be drawn without erasing the 
temperature lines. The contour interval is set to 5. The 
lines should all be drawn in green using color number 2, 
with labeling and line types as above. Only parameters 


wh i c h mu s t 

be changed fr om the above definitions are sh own 

GFUNC 

- 

dwp t 

LINE 

- 

2 = g r een/ -3/1/2 

CINT 

- 

5 

CLEAR 

- 

no 

MAP 

- 

0 

LATLON 


0 


3. Now clear the screen, scale the data by 10**5, draw 
contour lines of the divergence of the observed wind. 


4-6 



GEMPAK PROGRAMS 


alternating between colors 3 and 7. Use a default 
contour interval. Use heavy solid lines, labeling 
every contour line. 


GFUNC 

CLEAR 

SCALE 

CINT 

LINE 

MAP 


d i v ( obs ) 
ye s 
5 

3;7/l/7/l 

1 


ERROR MESSAGES 


[ GDCNTR 

+ 1] 

[GDCNTR 

-1] 

[GDCNTR 

-2] 

[GDCNTR 

-3] 

[GDCNTR 

-4] 

[GDCNTR 

-5] 

[GDCNTR 

-6] 


WARNING. There are no contour levels. 
Fatal error initializing TAE . 

Fatal error reading TAE parameters. 
Fatal error initializing GEMPLT . 

Grid requested is not available. 

Error setting grid navigation for file 
There are no grids in grid file. 
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4.3 PROGRAM GDCNVRT 


GDCNVRT converts GEMPAK3 grid files to new grid files. 

INPUT PARAMETERS 

OLDGRD * GEMPAK 3 grid file 

GDOUTF “Output grid file 

MAXGRD ‘Maximum number of grids 


PROGRAM DESCRIPTION 

GDCNVRT converts GEMPAK3 grid files to the new GEMPAK format. 
It is necessary to use this program since the GEMPAK programs 
are incompatible with GEMPAK3 files. 

Note that GEMPAK4 grid files do not need to be converted to 
be used in GEMPAK5 programs. However, because grid packing 
options have been added, GEMPAK5 grids cannot be used in 
GEMPAK 4 p r o g r ams . 

MAXGRD is the maximum number of grids that can be stored in 
the new file. If MAXGRD is less than the number of grids in 
OLDGRD, the new file will be created with MAXGRD equal to the 
number of grids in OLDGRD. 


EXAMPLES 

1. Convert the GEMPAK3 grid file, SFC.GRD to the new GEMPAK 
grid file, NEWSFC.GRD, which can contain a maximum of 500 
grids . 

OLDGRD - s f c . g r d 
GDOUTF - newsfc.grd 
MAXGRD - 500 


ERROR MESSAGES 

[ GDCNVRT -l] Fatal error initializing TAE . 

[ GDCNVRT - 2 j Fatal error reading TAE par ame t e r s . 

[GDCNVRT -3] Old file ... cannot be opened. 

[GDCNVRT -4] New file ... cannot be created. 
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4.4 PROGRAM GDCROSS 


GDCROSS draws vertical cross sections through a grid. 
INPUT PARAMETERS 


CXSTNS 

GDATTIM 

GVCORD 

GFUNC 

GVECT 

WIND 

GDFILE 

PTYPE 

YAXIS 

CINT 

SCALE 

LINE 

BORDER 

TITLE 

CLEAR 

DEVICE 

TEXT 

PANEL 


♦Cross-section station line 

*Gr id date/t ime 

♦Grid vertical coordinate 

♦Grid function 

♦Grid vector function 

♦Wind symbol / s i z/wd th/ typ/hdsz 

♦Gr id file 

♦Plot type/h:w ratio/margins 

♦Ystrt/ystop/yinc/lbl ;gln; tck 

♦Contour i n t e r va 1 /mi n/max 

♦Scaling factor 

♦Line c o 1 o r /d a sh/wi d th/ 1 abe 1 

♦Background c o 1 o r /d a sh/wi d th 

♦Title c o 1 o r / 1 i ne / t i 1 1 e 

♦Clear screen flag 

♦Graphics device 

♦Text size/font / h w flag 

♦Panel 1 oc a t i on/ co 1 o r /dash/wi d th 


PROGRAM DESCRIPTION 

GDCROSS draws a vertical cross section between two points 
in a grid field. The cross-section path is a line segment 
on the grid projection plane joining the two points. Gridded 
data are interpolated to the cross-section plane at intervals 
corresponding to approximately one grid increment. 

CXSTNS specifies the cross-section end points which must be 
separated by a >. Endpoints may be entered as grid coordinates 
separated by a semicolon and preceded with an a. Latitude 
and longitude specifications must be separated with a semicolon. 
Stations may be entered using either character identifiers or 
station numbe r s . 

The vertical coordinate, set in GVCORD, can be PRES, THTA, or 
HGHT. The vertical axis scaling, set in PTYPE, can be LIN, 

LOG, KAP or STUVE . STUVE and KAP are the same; SKEWT may 
not be entered. The plot aspect ratio and margins may also be 
entered in PTYPE. 

Both scalar and vector fields can be displayed in the cross- 
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section plane. Contour lines are drawn through the scalar 
field. Vector fields are depicted using arrows or barbs. 

Circulations can be displayed in the cross-section plane by 
specifying GFUNC = CIRC ( V, W ) , where V is a vector field 
and W is the scalar vertical motion in the GVCORD coordinate. 
The horizontal component of the circulation is the tangential 
component of V. The vertical component is W scaled up to 
account for the exaggerated aspect ratio of the display 
relative to that of the real atmosphere. For the CIRC 
operator, W is assumed to be pressure velocity in mb/s 
for PRES and THTA coordinates and cm/ s for the HGHT coordinate. 
If the vertical component does not require scaling, the 
circulation can be specified explicitly in the fo rm 
GVECT = VECR ( TANG ( V ), W )/VERT, where VERT is a flag 
indicating that the vector components are not to be north 
relative. 


EXAMPLES 

1. Plot temperature in Celsius on a log P chart along the 
cross section fr om BWI t o LAX . Plot the contours in 
color number 2, using thick, solid lines. Plot wind 
barbs for the observed wind. Plot the chart to 50 
millibars using the default labels. Accept default 
values for the other parameters. 

CXSTNS 
GDATTIM 
GVCORD 
GFUNC 
GVECT 
WIND 
GDF I LE 
PTYPE 
YAXIS 
CINT 
SCALE 
LINE 
BORDER 
TITLE 
CLEAR 
DEVICE 
TEXT 
PANEL 


2 . N ow , using the s ame file, plot a log pressure cross 
section of the vorticity advection, scaled by 10**5. 
Also, plot the ageostrophic circulation using arrows. 


- bwi>lax 

- last 

- pres 
= t emp 
* obs 

- bml 

- ngm.grd 

- log 

- /50 


- 2/1/7 


= yes 
= qm 


4-10 



GEMPAK PROGRAMS 


Use a cross-section line from grid point (1,1) to grid 
point (10,20). Use a he i ght - to-width plot ratio of 2. 


GFUNC 

GVECT 

GPOINT 

PTYPE 

SCALE 

WIND 


adv(avor(obs) jobs) 
c i r c ( age ; omeg ) 

® 1 ; 1 > 10 ; 20 

log/2 

5 

ami 


ERROR MESSAGES 

[GDCROSS +2] 
[GDCROSS +1] 
[GDCROSS -1] 
[GDCROSS -2] 
[GDCROSS -3] 
[GDCROSS -4] 
[GDCROSS -5] 
[GDCROSS -6] 
[GDCROSS -7] 
[GDCROSS -8] 
[GDCROSS -9] 
[GDCROSS -10] 
[GDCROSS -11] 
[GDCROSS -12] 
[GDCROSS -13] 
[GDCROSS -14] 
[GDCROSS -18] 
[GDCROSS -19] 
[GDCROSS -20] 


Surface value of GVCORD is not available. 

WARNING. There are no contour levels. 

Fatal error initializing TAE. 

Fatal error reading TAE par ame t e r s . 

Fatal error initializing GEMPLT . 

Input for CXSTNS is invalid. 

Input for GDATTIM is invalid. 

Input for GVCORD is invalid. 

Input for PTYPE is invalid. 

Error defining graph coordinates. 

No points found for cross section. 

LIN is not possible for this vertical coordinate. 
Invalid cross-section coordinates. 

No 1 eve Is-- check input for GDATTIM and GVCORD. 
©LEVEL in-line parameter is not allowed. 

%VCORD in-line parameter is not allowed. 

GVECT cannot be evaluated. 

GFUNC cannot be evaluated. 

Input for YAXIS is invalid. 
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4.5 PROGRAM GDDELT 


GDDELT deletes grids from GEMPAK grid files. 


INPUT PARAMETERS 

GDFILE *Grid file 

GDNUM *Grid numbers 


PROGRAM DESCRIPTION 

GDDELT deletes grids from GEMPAK grid files. 

The variable GDNUM allows the user to specify the grids to 
be deleted. GDNUM can be either a range of grid numbers or 
a list of individual grids. If GDNUM -ALL, all the grids 
in the file will be deleted. If GDNUM - LIST, all the grids 
in the file wi 11 be listed and the user wi 11 be p r omp ted to 
select the grids to be deleted. Note that if some grids 
are deleted, the grid n umb ers for the r ema ining grids ma y 
change. 


EXAMPLES 


1. Delete grids 3, 4, and 

5 f r om a grid file. 

GDFILE 

GDNUM 

ngm. grd 
3-5 


2. Delete grids 3, 6, and 

8 from SFC.GRD. 

GDNUM 

GDFILE 

3 ; 6 ; 8 
s f c . g rd 


3 . De 1 e t e all 

the grids from the file TODAY . GRD . 

GDNUM 

GDFILE 

all 

today . g rd 


ERROR MESSAGES 



[GDDELT -1] 
[GDDELT -2] 
[GDDELT -3] 
[GDDELT -4] 

Fatal error initializing TAE . 

Fatal error reading TAE par ame t e r s . 
Invalid grid range entered. 

Grid number is invalid. 
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[ GDDELT -5] There are no grids in file .... 
[GDDELT -6] No valid grids in list. 
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4.6 PROGRAM GDDIAG 


GDDIAG computes a scalar diagnostic grid and adds it to the 
grid file. 

INPUT PARAMETERS 


GDFILE 

GDOUTF 

GFUNC 

GDATTIM 

GLEVEL 

GVCORD 

GRDNAM 

GPACK 


♦Grid file 
♦Output grid file 
♦Gr id function 
♦Grid date/t ime 
♦Gr id level 

♦Grid vertical coordinate 
♦Output parameter name 
♦Number of bits/packing type 


PROGRAM DESCRIPTION 

GDDIAG computes a diagnostic grid and saves the result in 
a GEMPAK grid file. 

The input grids for the computation must be in GDFILE. The 
resulting grid is added to the output file, GDOUTF. The 
input and output files may be the same. 

The name of the output grid is given in GRDNAM. If this 
value is blank, the name generated by the grid diagnostic 
package will be used. GPACK defines grid packing. 

If GPACK is blank, no packing will be done. Grid packing 
details are documented in the GPACK variable documentation. 


EXAMPLES 

1. Compute the divergence of the observed wind at 850 mb 
from the 48-h forecast for the last model run and 
save the result in MYFILE.GRD. Use the default name 
for the grid. Do not pack the grid. 


GDFILE 

GDOUTF 

GFUNC 

GDATTIM - 

GLEVEL 

GVCORD 

GRDNAM 


ngm. grd 

my f i 1 e . g r d 

d i v fobs] 

f 4 8 

850 

pres 
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GPACK 


2. Compute the magnitude of the gradient of the dewpoint 
depression at 500 mb for the last time in the grid 
file and store the result in the input file. Name 
the output grid DEW3RAD. Use GRIB packing with 
16 bits per grid point . 

GDFILE - ngm.grd 
GDOUTF = ngm.grd 

GFUNC - mag ( grad ( sub (temp.dwpt) ) ) 

GDATTIM - last 
GLEVEL - 500 

GVCORD = pres 
GRDNAM = dewgrad 
GPACK = 1 6 / g r i b 


ERROR MESSAGES 

[GDDIAG -1] Fatal error initializing TAE . 

[GDDIAG -2] Fatal error reading TAE parameters. 

[GDDIAG -3] Error initializing GEMPLT . 

[GDDIAG -4] There is no input file specified. 

[GDDIAG -5] Different navigation in input, output files. 

[GDDIAG -6] Error opening input files. 

[GDDIAG -7] Error writing output grid. 

[GDDIAG -8] Only one output file is permitted. 

[GDDIAG -9] Output file open failure. 

[GDDIAG -10] Grid files have different navigations. 
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4.7 PROGRAM GDEDIT 


GDEDIT reads grids from a sequential edit file and adds them 
to a GEMPAK grid file. 

INPUT PARAMETERS 

GDEF IL *Grid edit file 

GDFILE *Gr i d file 

GPACK *Number of bits/packing type 


PROGRAM DESCRIPTION 

GDEDIT reads grids from a sequential text edit file and adds 
them to a GEMPAK grid file. The input edit file can be 
created using a text editor or using the OUTPUT = F option 
in GDLIST. 

Certain header information must be included before the grid 
data. Required information includes the grid size and the 
grid identifier. Scaling information may also be included. 

The grid size must be included on a line: 

COLUMNS: 1 33 ROWS: 1 14 

The edit grid must begin at column - 1 and row - 1 and 
include the entire grid. WARNING: the edit grid size must 
match the size specified in the grid file. Otherwise, a 
fatal error is generated. 

Th e grid identifier is found on a line such as: 

841227/1200F48 850 PRES DIVOBS 

The first line containing a slash ( / ) is assumed to be a 
gr id i den t i f i e r . The order of values is: time 1, time 2 
(optional), level 1, level 2 (optional), vertical 
coordinate, parameter name. 

If the data were scaled before being listed, scaling 
information must be included. This is of the form: 

Scale factor: 10**5 

This information is found by searching for "factor” and 
••**» ° 
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Data follow this header information. The rows are assumed 
to read from the last (top) to the first. Each row of 
data c on t a i n s : 

ROW 14 0.27 ... 

The word ROW and the row number must be included. The row 
number is NOT checked to see that this is the correct row. 
Data for each row may wrap to the next line. 

An e x amp le of a grid edit file foil ows : 

Grid file: SMALL . GRD 

GRID: TIME1 TIME2 LEVEL 1 LEVEL2 VCORD PARM 

841225/0000 0 PRES TEMP 


AREA: DSET 
COLUMNS : 1 

7 

ROWS: 

1 

GRID 

6 

SIZE: 

7 

Scale factor : 

10**0 






COLUMN: 

1 

2 

3 

4 

5 

6 

7 

ROW 6 

-12.81 

-12.40 

-11.98 

-11.30 

-10.82 

-11.18 

-12.25 

ROW 5 

-9.89 

-9.90 

-9.41 

-8.26 

-7.00 

-6.43 

-6 . 80 

ROW 4 

-2 . 14 

-2 . 94 

-3.66 

-3.43 

-2.36 

-1 .29 

-0.92 

ROW 3 

4.12 

4.25 

3.23 

2 . 35 

2.63 

3 .51 

3.98 

ROW 2 

6.39 

7 . 35 

7.03 

6.29 

6 .31 

6.85 

7 . 14 

ROW 1 

7.26 

7 . 95 

7 . 80 

7 .57 

7.87 

8.49 

8 . 96 

EXAMPLES : 









1 . Ad d the grids in the edit file GDL I ST .FIL to the grid 
file 25DEC12Z .GRD. Do not pack the data. 


GDEFIL - gdlist.fil 
GDFILE - 25decl2z.grd 
GPACK 

2. Add the grid in TMPC. FIL to the grid file NMC.GRD. Pack 
the data using the GRIB format. Select the number of 
bits so that two digits after the decimal place will be 
saved . 
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GDEF I L - tmp c . f i 1 
GDFILE - nmc.grd 
GPACK = +2/de c 


ERROR MESSAGES 


[GDEDIT 

-1] 

[GDEDIT 

-2] 

[GDEDIT 

-3] 

[GDEDIT 

-4] 

[GDEDIT 

-5] 

[GDEDIT 

-6] 

[GDEDIT 

-7] 

[GDEDIT 

-8] 


Fatal error initializing TAE. 

Fatal error reading TAE parameters. 

Valid grid identifier not found. 

Grid size cannot be determined. 

The grid size in the edit file is incorrect. 
Error reading grid data. 

The edit file ... is invalid. 

The output grid file ... is invalid. 
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4.8 PROGRAM GDINFO 


GDINFO lists information about GEMPAK grid files. 


INPUT PARAMETERS 

GDFILE ‘Grid file 

L STALL ‘Full list flag 

OUTPUT ‘Output device (TPF ) 


PROGRAM DESCRIPTION 

GDINFO lists information about GEMPAK grid files. 

The navigation information and grid analysis information will 
be listed. If requested, the grids in the file will also be 
listed. If LSTALL - YES, all the grids in the file will be 
listed. If LSTALL - NO, only the navigation and analysis 
information will be displayed. 


EXAMPLES 

1. List information about grid file NGM.GRD. All of the 
grids are to be listed. The listing is to be at the 
user’s terminal and to a file named ngm.list. 


GDFILE 

LSTALL 

OUTPUT 


ngm. g rd 
ye s 

tf /ngm. list 


ERROR MESSAGES 

[GDINFO -1] 
[GDINFO -2] 


Fatal error initializing TAE. 

Fatal error reading TAE parameters. 
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4.9 PROGRAM GDLIST 


GDL I ST lists data fr om a scalar grid field. 


INPUT PARAMETERS 


GDATTIM 

GLEVEL 

GVCORD 

GFUNC 

GDFILE 

GAREA 

PROJ 

SCALE 

OUTPUT 


•Grid date/time 
*Gr id level 

•Grid vertical coordinate 
*Gr id f unc t i on 
‘Grid file 
•Graph i c s area 

•Map p r o j e c t i on/ ang 1 e s /ma r g i ns 
•Scaling factor 
•Output device (TPF) 


PROGRAM DESCRIPTION 


GDLIST lists a diagnostic grid computed from the grids in a 
GEMPAK gridfile. 


The data may be listed over a subgrid. The variables GAREA 
and PROJ specify the data subset area. This area is only 
approximate; the subgrid selected will cover the subset area 
GAREA must be specified as a geographic area, as an area 
^l ere d on a station, as 1 a t i t ude / 1 ong i tude bounds or as 
DSET or GRID. If DSET or GRID^s chosen, the entire grid 
will be printed. 


If output is sent to the printer, the full 132 columns will be 
used. Output to the terminal or to a file will be 80 columns 

rnriQT is Sent t0 a file ’ thc file wi 1 1 be named 

GDLIST. FIL by default. 


EXAMPLES 


1 . 


List surface temperature from the April 5 00Z NGM cycle for 
the O-h forecast time. The output will be to the terminal. 
Ihe data subset area is a zoomed area centered on Colorado. 


GDATTIM - 

GLEVEL 

GVCORD 

GFUNC 

GDFILE 

GAREA 


0405/00f00 

0 

pres 
tmpf 
n gm . g r d 
co** 
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2 . 


PROJ 

- 

me r 

SCALE 

- 


OUTPUT 

“ 

t 

Compute 

the 

divergen 

Scale the data by 10 

called 

d i v . 

us . 

GFUNC 


div(obs) 

GAREA 

- 

us 

PROJ 


c ed 

SCALE 

- 

5 

OUTPUT 

- 

f /d i v . us 


and list the output over 
and write the output to 


the US. 
a file 


ERROR MESSAGES 

[GDLIST -1] 
[GDLIST -2] 


Invalid area type. 

Fatal error reading TAE par ame t e r s . 
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4.10 PROGRAM GDMAP 


GDMAP plots data from a GEMPAK grid data file. 


INPUT PARAMETERS 


GDATTIM 

GLEVEL 

GVCORD 

GFUNC 

GDFILE 

GAREA 

POINTS 

POSN 

COLORS 

MARKER 

MAP 

LATLON 

PANEL 

TITLE 

SCALE 

DEVICE 

PROJ 

CLEAR 

TEXT 


•Grid date/time 
•Gr id level 

•Grid vertical coordinate 

•Gr id function 

•Grid file 

•Graphics area 

•XI ;yl ;xinc ; y i nc 

• Po s i t i on 

•Co lor list 

•Marker c o 1 o r / type / s i ze /wid th/hw 

•Map c o 1 o r /d a sh/wi d th 

•Line c o 1 o r/da sh/wi d th/ 1 abe 1 / inc 

•Panel 1 oca t i on/ co 1 o r/da sh/wi d th 

•Title c o 1 o r / 1 i ne / t i t 1 e 

•Seal ing factor 

•Graphics device 

•Map p r o j e c t i on/ ang 1 e s /ma rg ins 

*C1 ear screen flag 

•Text size/font/hw flag 


PROGRAM DESCRIPTION 


GDMAP plots data computed from GEMPAK grid files on a map 


The variable, POSN, is used to select 
data relative to the grid point using 
numbe r s : 


the position for the 
the following position 


7 

1 3 

2 0 4 

5 6 

8 


Position Owill plot data centered at the station, 
invalid position number is entered, position 0 will 


If an 
be used. 


Th* J ariab ^?I * s used to scale the data by 10 ** SCALE 

The data will be scaled and then rounded to the nearest 
integer before it is plotted. 
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POINTS allows the user to select points for plotting. Four 
numbers, separated with semicolons, may be entered: t e irs 

grid point in the x direction, the first grid point in the y 
direction, the increment in the x direction and the increment 
in the y direction. The points are numbered from the start of 
the grid, not the graphics area. If any number is not entered, 
a value of 1 is assumed. Thus, if POINTS = 1, all points wi 
be plotted. Note the grid point (1,1) is the lower left corner 

of the grid. 


A map and title may 


also be included. 


EXAMPLES 


The temperature field at 850 mb will be plotted in color 
at position number 7 which is directly above the station. 
The data will be plotted at every point for an area 


centered on California, 
plotted at the station. 


Marker number 5, a cross 
The title, "Temperature 


will be 
wi 11 be 


2 . 


wr i 1 1 en 

i n 

color 6 . 

GDATTIM 

* 

303/12 

GLEVEL 

= 

850 

GVCORD 

= 

pres 

GFUNC 

- 

t emp 

GDF I LE 

- 

ngm. grd 

GAREA 

« 

c a 

POINTS 

- 

1 

POSN 

- 

7 

COLORS 

- 

3 

MARKER 

- 

1/5 

MAP 

- 

2/2 

LATLON 

- 


PANEL 

- 

0 

TITLE 

- 

6 / / T emp e r a t u r e 

SCALE 

- 

0 

DEVICE 

- 

sg 

PROJ 

- 

me r 

CLEAR 

- 

yes 

TEXT 

* 


Ad d the 

d e wp oint depressio 

wh i c h i 

s directly bel ow t h 

and title 

are turned off s 

GFUNC 


sub [ t emp ; dwp t ] 

POSN 

- 

8 

MARKER 

= 

0 

MAP 

- 

0 

TITLE 

- 

0 


to the plot 
grid point 


in position 8, 

The map . markers 


since they are already plotted. 
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CLEAR — no 

3. Now clear the screen and plot the relative vorticity of 

the observed wind, scaled by 10**5, centered at the station 
in color 1. Add a map and title in color 1. Omit the 
markers since the data will overlay them. 

GFUNC 
COLORS 
MAP 
SCALE 
MARKER 
CLEAR 
POSN 

ERROR MESSAGES 

[GEMAP -1] Fatal error initializing TAE. 

[GDMAP -2] Fatal error reading TAE parameters. 

[GEMAP -3] Fatal error initializing GEMPLT. 


- vor [ obs ] 

- 1 
- 1 

- 5 

- 0 

- yes 

- 0 
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4.11 PROGRAM GEMOD 


GEMOD moves data from one GEMPAK grid file to another. 


INPUT PARAMETERS 

GDFILE *Gr i d file 

GDOUTF *Output grid file 

GDNUM *Gr i d numbers 

GPACK *Numbe r of bits/packing type 


PROGRAM DESCRIPTION 

GEMOD moves grids from an input grid file to an output 
file. The input grid file name is specified in GDFILE. The 
output grid file name is specified in GDOUTF . The input and 
output grid files must have the same navigation information. 

The numbers of the grids to be moved are specified in GDNUM. 
Either a list of grid numbers, with semicolons separating the 
numbers, or a range of grid numbers, with the first and last 
n umb ers separated by a dash, ma y be entered. 

The output grids will be packed using the information in 
GPACK. 


EXAMPLES 

1. Add grids 1, 4 and 7 from file, OLD.GRD, to file NEW. GRD. 
Pack the data using 16 bits for each data value. 

GDFILE - old.grd 
GDOUTF - new.grd 
GDNUM - 1 ; 4 ; 7 

GPACK - 1 6 / g r i b 


2. Copy all the grids in TODAY. GRD to DEC. GRD. Do not 
pack any data. 

GDFILE - today. grd 
GDOUTF = dec. grd 
GDNUM - all 
GPACK 
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ERROR MESSAGES 

[GDMOD -1] 
[GDMOD -2] 
[GDMOD -3] 
[GDMOD -4] 
[GDMOD -5] 
[GDMOD -6] 
[GDMOD -7] 
[GDMOD -8] 
[GDMOD -9] 
[GDMOD -10] 
[GDMOD -11] 


Fatal error initializing TAE. 

Fatal error reading TAE parameters. 
Invalid grid range entered. 

Grid number ... is invalid. 

There are no grids in file .... 

No valid grids in list. 

Error reading input grid. 

Error writing output grid. 

Error opening files. 

Files contain different navigations. 
Too many grids to transfer. 
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4.12 PROGRAM GDPROF 


GDPROF draws profiles of gridded data. 
INPUT PARAMETERS 


GDATTIM 

GVCORD 

GFUNC 

GPOINT 

GDFILE 

GVECT 

LINE 

MARKER 

BORDER 

PTYPE 

SCALE 

XAXIS 

YAXIS 

WIND 

WINPOS 

FILTER 

TITLE 

PANEL 

CLEAR 

TEXT 

DEVICE 

OUTPUT 

THTALN 

THTELN 

MIXRLN 


♦Grid date/time 

♦Grid vertical coordinate 

♦Gr id f unc t i on 

♦Grid point 

♦Grid file 

♦Grid vector function 
♦Line c o 1 o r /d a sh/wi d th/ 1 abe 1 
♦Marker co 1 o r / type / s i z e /wi d th/hw 
♦Background c o 1 o r /da sh/wi d th 
♦Plot type/h:w ratio/margins 
♦ Sc a 1 i ng factor 

♦Xstrt/xstop/xinc/lbl ;gln; tck 

♦Ystrt/ystop/yinc/lbl ; g In ; t ck 

♦Wind symbo 1 / s i z /wd t h/ t yp/hd s z 

♦Wi nd position 

♦Filter data flag 

♦Title c o 1 o r / 1 i ne / t i t 1 e 

♦Panel 1 oc a t i on/ co 1 o r /da sh/wi d th 

♦Clear screen flag 

♦Text size/font/hw flag 

♦Graph ics device 

♦Output device (TPF) 

♦THTA color/dash/width/mn/mx/inc 
♦THTE color/dash/width/mn/mx/inc 
♦MIXR color/dash/width/mn/mx/inc 


PROGRAM DESCRIPTION 

GDPROF draws vertical profiles at a point. Data from all 
levels in the grid file are interpolated to the point. 

GPOINT specifies the point for the profile. It may be 
entered as a set of x and y grid points separated with a 
semicolon and preceded with an @, as a latitude and 
longitude separated with a semicolon, or as a station 
character identifier or station number. 

Th e vertical coordinate, set in GVCORD , can be P RE S , THTA , o r 
HGHT . The vertical axis type, set in PTYPE, can be LIN, 

LOG, KAP, or SKEW. Only temperatures may be plotted on a 
skew T plot. The he i gh t - 1 o -wi d th ratio of the plot may 
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be specified in PTYPE following a /. If no value is entered 
a value of 0 will be used. 

The value in GFUNC is computed at every level in the data set. 
No errors are returned if the function cannot be computed at 
any level, except at the surface. If data at the surface 
cannot be computed, an error message is written and 
plotting continues. GVECT specifies a vector to be 
plotted in the right margin. 

If GFUNC is a temperature or dewpoint temperature, background 
lines of potential temperature, equivalent potential 
t emp erature and mixing ratio can be displayed by specifying 
THTALN , THTELN , and MIXRLN . If these background lines are 
requested and GFUNC is not a temperature, an error message 
is generated. 


EXAMPLES 

1. Plot temperature in Celsius on a log P chart at IAD for 

the latest time in the file. Use the default line, marker, 
title and background colors. Label the Y axis from 1000 mb 
to 100 mb in increments of 100 mb. Label the X axis from 
-40 degrees to +40 degrees in increments of 10 degrees. Add 
all the background lines as dotted lines. Plot wind barbs in 
the right ma r g i n . 


GDATTIM - 

GVCORD 

GFUNC 

GPOINT 

GDFILE 

GVECT 

LINE 

MARKER 

BORDER 

PTYPE 

SCALE 

XAXIS 

YAXIS 

WIND 

WINPOS 

FILTER 

TITLE 

PANEL 

CLEAR 

TEXT 

DEVICE 

OUTPUT 

THTALN 


last 
pres 
t emp 
IAD 

ngm. g rd 
obs 
1 
1 
1 

LOG 

0 

-40/40/10 

1000 / 100/100 

bml 

YES 

1 

0 

ye s 

vt 

no 

1/10 
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THTELN - 1/10 

MIXRLN - 1/10 

2. Using the same file, plot the same parameters at HAT. 

The winds will be plotted to the right of the winds 
plotted above. 

GPOINT - hat 

BORDER - 0 

CLEAR - no 

THTALN - 0 

THTELN - 0 

MIXRLN - 0 

3. Clear the screen and plot the relative humidity at BWI in 
color number 2 using a dotted line. Do not plot any 
winds. 

GFUNC - r e 1 h 

GVECT 

LINE = 2/10 

GPOINT - bwi 


ERROR MESSAGES 


[GDPROF 

+ 1] 

Inva 1 i 

d parameter for background lines 

[ GDPROF 

-1] 

Fatal 

error initializing TAE. 

[GDPROF 

-2] 

Fatal 

error reading TAE parameters. 

[GDPROF 

-3] 

Fatal 

error initializing GEMPLT . 

[GDPROF 

-4] 

Input 

for GPOINT is invalid. 

[GDPROF 

-5] 

Input 

for GDATT IM is invalid. 

[GDPROF 

-6] 

Input 

for GVCORD is invalid. 

[GDPROF 

-7] 

Input 

for PTYPE is invalid. 

[GDPROF 

-8] 

Error 

defining graph coordinates. 

[GDPROF 

-9] 

No po i 

nts found for profile. 

[GDPROF 

-10] 

There 

are no levels at this time. 

[GDPROF 

-11] 

Input 

... for y axis is invalid. 

[GDPROF 

-12] 

I npu t 

... for x axis is invalid. 
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4.13 PROGRAM GDSTAT 


GDSTAT computes statistics on a time series of grids 


INPUT PARAMETERS 


GDFILE 

♦Grid file 

GDOUTF 

♦Output grid file 

GDATTIM 

♦Grid date/time 

GLEVEL 

*Gr id level 

GVCORD 

♦Grid vertical coordinate 

GFUNC 

♦Grid function 

GRDNAM 

♦Output parameter name 


PROGRAM DESCRIPTION 

GDSTAT computes the average, standard deviation, and number 
of reporting times at each grid point for a time series of 
grids. The output grids will be written to the output grid 
file. 

The input grid files are given in GDFILE. The output 
file is given in GDOUTF . There may be one or two input 
files. The output file may also be an input file. 

GFUNC specifies the grid diagnostic function to be computed. 
GDATTIM must be a time range. The output grid names will 
be the name in GRDNAM prefixed by AVG, STD, and CNT. 

If GRDNAM is blank, the default name from the grid diagnostic 
computation will be used. 


EXAMPLES 


1. Compute the grid statistics for the 850 mb vorticity 
field using all the times in the grid file. Name 
the output grids AVGVOR, STDVOR and CNTVOR and write 
them into the same file. 


Write the output grids into the same file. 


GDFILE 

GDOUTF 

GFUNC 

GVCORD 

GLEVEL 

GDATTIM 


ngm. g rd 
ngm. grd 
vor ( obs ) 
pres 
850 
all 
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GRDNAM = vo r 


2. Compute the grid statistics for the surface temperature 
data using the data for every 12 hours from Jan 1 to 
Jan 15. The input and output files are different. The 
default names will be used. 


GDF I LE 

GDOUTF 

GFUNC 

GVCORD 

GLEVEL 

GDATTIM 

GRDNAM 


jan90sfc.grd 

s t a t . g r d 

tmpc 

none 

0 

900101/0000-900115/0000-12 


ERROR MESSAGES 


[GDSTAT 

-1] 

[ GDSTAT 

-2] 

[GDSTAT 

-3] 

[GDSTAT 

-4] 

[GDSTAT 

-6] 

[GDSTAT 

-7] 


Fatal error initializing TAE . 

Fatal error reading TAE parameters. 
Fatal error initializing GEMPLT . 

No valid times found in grid file. 
There are fewer than four grids. 

Error writing grid ... to output file. 
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4.14 PROGRAM GDSTREAM 


GDSTREAM draws streamlines through a GEMPAK vector grid. 


INPUT PARAMETERS 


GDATT IM 

*Gr i d 

date/ t ime 

GLEVEL 

*Gr i d 

level 

GVCORD 

*Gr i d 

vertical coordinate 

GVECT 

*Gr i d 

vector function 

GDFILE 

*Gr id 

file 

GAREA 

♦Graph 

i c s area 

LINE 

♦Line 

color/dash /wi dth/label 

MAP 

♦Map col or/dash/width 

LATLON 

♦Line 

co 1 o r /da sh/wi dth/label/inc 

TITLE 

♦Title 

color/line /title 

PANEL 

♦Panel 

location/ color/ da sh/wi d t h 

TEXT 

♦Text 

size/font / hw flag 

DEVICE 

♦Graphics device 

PROJ 

♦Map p r o j e c t i on/ang 1 e s /mare i ns 

CLEAR 

♦Cl ear 

screen flag 


PROGRAM DESCRIPTION 

GDSTREAM draws streamlines through any vector grid computed 
by the grid diagnostics package. 

Th e vector field is specified in GVECT . S t r e ami i n e s are dr awn 
using a modification of the NCAR streamline program. 

The line color, dash pattern and line width are set in LINE. 

A map, title and 1 a t i tude/ 1 ong i tude lines may also be included. 


EXAMPLES 

1. Draw streamlines through the observed wind field at 850 
mb using color 3. Add a Mercator map in color 1 and 
linetype 2 and a title in color 2. Center the graphics 
area on Wisconsin. 

GDATT IM — last 

GLEVEL - 850 

GVCORD - pres 

GVECT - o b s 

GDFILE — ngm.grd 

GAREA - wi 
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LINE - 3 

MAP - 1/2 

LATLON 

TITLE - 2 

PANEL - 0 

TEXT 

DEVICE - vt 

PROJ - me r 

CLEAR - yes 

2. Clear the screen and draw the streamlines of the gradient 
of the t emp e r a t u r e . 

GVECT - g r a d ( t emp ) 

CLEAR - yes 


ERROR MESSAGES 

[ GDSTREAM -1] Fatal error initializing TAE. 
[GDSTREAM -2] Fatal error reading TAE parameters. 
[GDSTREAM -3] Fatal error initializing GEMPLT . 
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4.15 PROGRAM GDTSER 


GDTSER draws a time series of a scalar at a single level. 
INPUT PARAMETERS 


GDATTIM 

GLEVEL 

GVCORD 

GFUNC 

GPOINT 

GDFILE 

PTYPE 

TAXIS 

YAXIS 

BORDER 

LINE 

MARKER 

TITLE 

CLEAR 

SCALE 

PANEL 

DEVICE 

TEXT 


♦Grid date/time 
♦Grid level 

♦Grid vertical coordinate 
♦Grid function 
♦Grid point 
♦Gr id file 

♦Plot type/h:w ratio/margins 
♦Time 1 - t ime 2-tinc;lbl;gln;tck 
♦Ystrt/ystop/yinc/lbl ;gln; tck 
♦Background c o 1 or /da sh/wid th 
♦Line co 1 o r /da sh/wi d th/ 1 abe 1 
♦Marker co 1 o r / type/ s i z e /wi d th/hw 
♦Title co 1 o r / 1 i ne / t i t 1 e 
♦Clear screen flag 
♦Scaling factor 

♦Panel 1 oca t i on/ c o 1 o r /da sh/wi d th 

♦Graphics device 

♦Text size/font / hw flag 


PROGRAM DESCRIPTION 

GDTSER draws a time series of a parameter at a fixed level. 

The variable may be any scalar computed by the grid diagnostics 
package . 

GPOINT specifies the location for the time series. It may 
be entered as a set of x and y grid points separated with a 
semicolon and preceded with an @, as a latitude and 
longitude separated with a semicolon, or as a station character 
identifier or station number. 


The times to plot are specified in GDATTIM. Only those 
times specified will be plotted. Lines will be drawn 
connecting the data points provided no more than two 
points are missing between segments. Note that the times 
may represent a series of base times with the same or no 
forecast time, or a series with the same base time and 
a list of forecast times. 

The time axis is specified in TAXIS using the usual GEMPAK 
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date/time conventions, including FIRST and LAST. If 
TAXIS has no increment, a reasonable value is chosen. 

If TAXIS is blank, the limits are taken from the first and 
last valid times set in GDATTIM. 

If the user desires mu ltiple lines on the s ame graph , 

TAXIS and YAXIS should be explicitly set; then LINE, MARKER, 
and TITLE can be varied for successive combinations of GLEVEL , 
GVCORD , GPOINT , and GFUNC . 

The height-to-width ratio of the plot and the margins may 
be specified in PTYPE following a /. Note that the first 
value in PTYPE, the axis type, is ignored. 


EXAMPLES 


1. Plot the 00 h forecast temperatures in Celsius at the 
500 mb level at BWI initialized at times between 26/12 
and 31/00 in the file. Use the default line, marker, 
title and background colors. Label the Y axis from 
-40 to 0 in increments of 10 degrees. Label the time 
axis from 26/00 to 31/12 in increments of 12 h. Make 
the plot half as high as wide. 


GDATTIM 

GLEVEL 

GVCORD 

GFUNC 

GPOINT 

GDFILE 

PTYPE 

TAXIS 

YAXIS 

BORDER 

LINE 

MARKER 

TITLE 

CLEAR 

SCALE 

PANEL 

DEVICE 


26/12f00-31/00f00 

500 

pres 

tmpc 

bwi 

ngm. g rd 
/ .5 

26/00-31/12-12 

-40/0/10 

1 

1 

1 

1 

ye s 
999 
0 

vt 


2. Using the same file, plot a time series of the 800 mb 
geostrophic wind speed at latitude and longitude of 39 
and -77 for all forecasts initialized at 29/00. Plot 
the graph with the height/width ratio of 0.3 in the top 
half of the page and use internal scaling for both axes. 

GDATTIM - 29/OOf all 
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TAXIS 

GLEVEL - 800 

GFUNC - ma g ( g e o ) 

GPOINT - 39; -77 

PTYPE = / . 3 

YAXIS 

PANEL - t 


ERROR MESSAGES 

[GDTSER -1] Fatal error initializing TAE. 

[GDTSER -2] Fatal error read ing TAE parameters. 

[GDTSER -3] Fatal error initializing GEMPLT. 

[GDTSER -4] Input for GPOINT is invalid. 

[GDTSER -5] Input for GLEVEL is invalid. 

[GDTSER -6] Input for GVCORD is invalid. 

[GDTSER -7] Input for PTYPE is invalid. 

[GDTSER -8] Error defining gr aph coordinates. 
[GDTSER -9] No points found for plot of .... 

[GDTSER -10] Some points missing for plot of .... 

[GDTSER -11] Input for ... in TAXIS i s inva lid. 

[GDTSER -12] Time range is size zero. 

[GDTSER -13] No grids found with times in GDATTIM. 
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4.16 PROGRAM GEWIND 


GEWIND draws wind barbs or arrows at each grid point of a vector 
field. 


INPUT PARAMETERS 


GDATTIM 

*Gr id date/t ime 

GLEVEL 

*Gr id level 

GVCORD 

♦Grid vertical coordinate 

GVECT 

*Grid vector function 

GDFILE 

♦Grid file 

GAREA 

♦Graphics area 

POINTS 

♦XI ; y 1 ; x i nc ;yinc 

WIND 

♦Wind symbo 1 / s i z /wd t h/ t yp/hd s z 

MAP 

♦Map co lor /dash/width 

LATLON 

♦Line c o 1 o r /d a sh/wi d th/ 1 abe 1 / inc 

PANEL 

♦Panel 1 oc a t i on/ c o 1 o r /da sh/wid th 

TITLE 

♦Title c o 1 o r / 1 i ne / t i t 1 e 

DEVICE 

♦Graphics device 

PROJ 

♦Map p r o j e c t i on/ ang 1 e s /ma r g i ns 

CLEAR 

♦Clear screen flag 

TEXT 

♦Text size/font / hw flag 


PROGRAM DESCRIPTION 

GEWIND draws wind barbs or arrows at each grid point of a 
vector field. 

POINTS allows the user to select points for plotting. Four 
numbers, separated with semicolons, may be entered: the first 

grid point in the x direction, the first grid point in the y 
direction, the increment in the x direction and the increment 
in the y direction. The points are numbered from the start of 
the grid, not the graphics area. If any number is not entered, 
a value of 1 is ass ume d . Thus, if POINTS “ 1, all points will 
be plotted. The grid point (1,1) is the lower left corner of 
the grid. 

Note that in this progr am , a M or K entered in WI ND to indicate 
meters/sec or knots will be ignored. By default, wind vectors 
will be computed in meters/sec. The operand KNOTV will convert 
the output to knots. Thus OBS will be computed in meters/sec 
while KNOTV (OBS) will be converted to knots. 


EXAMPLES 


4-37 



GEMPAK PROGRAMS 


Draw wind barbs at every grid point. The 
the observed wind at 850 mb on March 3 at 
is centered on Colorado. 


vector field is 
1 2Z . Th e plot 


GDATTIM - 

GLEVEL 

GVCORD 

GVECT 

GDF I LE 

GAREA 

POINTS 

WIND 

MAP 

LATLON 

PANEL 

TITLE 

DEVICE 

PROJ 

CLEAR 

TEXT 


303/12f00 
850 
pres 
ob s 

ngm. grd 
c o 
1 

bml 

1/10 

1/10 

0 

1 

vt 
me r 
ye s 


2. Using the values of the variables 
arrows at every other grid point, 
is the temperature gradient. The 
cleared. The arrows are in color 

GVECT — g r a d ( t emp ) 

POINTS - 1 ; 1 ; 2 ; 2 

WIND - am? 

CLEAR — no 

3. Plot wind barbs of the observed wind in knots. Note that 
WIND — BK1 and WIND ** BM1 give the same results. The 
units can only be changed in GVECT. 


supplied above, draw 
The vector to be used 
screen should not be 
7 . 


GVECT 

WIND 

kno t v( 
bkl 

obs ) 


ERROR MESSAGES 




[GEWIND -1] 
[GEWIND -2] 
[GEWIND -3] 

Fatal 

Fatal 

Fatal 

error 

error 

error 

initializing TAE. 
reading TAE par ame t e r s 
initializing GEMPLT. 
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4.17 PROGRAM GPAOIM 


GPAOIM sets up a view region for an AOIPS image. 

INPUT PARAMETERS 

GAREA ‘Graphics area 

PROGRAM DESCRIPTION: - 

This program adjusts the view region to correspond to an image 
size on an International Imaging Systems (II) display device. 
This program must be run for images which are less than 512 x 
512. 

GAREA is the file name of the image to be displayed. 

Note that after this program is run, PANEL must be set to 
blank, so that the view region will not be reset when plotting 
on the image . 

ERROR MESSAGES 


[GPAOIM 

-1] 

Fatal 

error 

initializing TAE. 

[GPAOIM 

-21 

Fatal 

error 

reading the TAE parameters 

[GPAOIM 

-31 

Fatal 

error 

initializing GEMPLT . 

[GPAOIM 

-4] 

Error 

reading image label. 
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4.18 PROGRAM GPAREA 


GPAREA updates the user’s geographic area table. 


INPUT PARAMETERS 

AREA *Data area 

GNAME ‘Geographic abbreviation/name 

PROGRAM DESCRIPTION 

This program allows a user to create and update custom 
geographic files and to query the bounds of a particular 
area . 

The user’s geographic area table is a file call GEOG.TBL. 

If this file does not exist, it will be created. In order 
to access areas from this file in any GEMPAK program, a # 
must be prefixed to the area n ame . 

In this program, AREA is entered as in any GEMPAK program. 
GNAME is the name to be associated with the area. This 
may be followed by a slash and a longer name which is only 
used for reference. 

If GNAME is blank, this program will simply display the 
information about the area defined in AREA. 

Note that the central latitude and longitude and the 
latitude and longitude range are the numbers actually stored 
inthefile. 


EXAMPLES 

1. Get the area called EAST from the GEMPAK table and add 
this area to the individual table with the name EASTUS . 

GNAME »* e a s t u s 
AREA - east 

2. Add the area BIGUS to the custom file. 

GNAME » bigus/large usa area 
AREA - 2 3 ; - 1 24 ; 52 ; - 60 

3. Query the bounds of BIGUS in the custom file. 
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GNAME - 
AREA - #b i gus 


ERROR MESSAGES 

[GPAREA -1] 
[GPAREA -2] 


Fatal error initializin 
Fatal error reading TAE 


TAE. 

par ame t e r s . 


4-41 



GEMPAK PROGRAMS 


4.19 PROGRAM GPBOX 


GPBOX draws a box around a region. 


INPUT PARAMETERS 

LINE *Line co 1 o r /da sh/wi d th/ 1 abe 1 

REGION *Type of region 

DEVICE *Graphics device 


PROGRAM DESCRIPTION 

This program draws a box around a region on the current 
graphics device. REGION can be specified as VIEW, PLOT, or 
DEVICE. VIEW is the view region of the graphics device. 

The view region may be changed in GEMPAK programs using the 
parameter PANEL. The PLOT region is the area used for the 
data plot, excluding margins. DEVICE is the entire device 
space . 


EXAMPLES 

1. Draw a solid line of width 1 around the view region on the 
VT devi c e . 

LINE - 1 

REG I ON — v i ew 
DEVICE - vt 

2. Draw a line around the device region in color 2, using line 
type 3 and width 5 . 

LINE - 2/3/5 

REGION — device 


ERROR MESSAGES 

[GPBOX -1] Fatal error initializing TAE . 

[GPBOX -2] Fatal error in reading the TAE parameters. 

[GPBOX -3] Fatal error initializing GEMPLT. 
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4.20 PROGRAM GPCLEAR 


GPCLEAR clears the current graphics device. 


INPUT PARAMETERS 

DEVICE ’Graphics device 


PROGRAM DESCRIPTION 

This program clears the graphics device. On an interactive 
device, the screen is erased. On a hardcopy device, the 
next plot will be on a new page. 


ERROR MESSAGES 


[GPCLEAR 

-1] 

Fatal 

error 

initializing TAE. 

[GPCLEAR 

-2] 

Fatal 

error 

reading TAE parameters. 

[GPCLEAR 

-3] 

Fatal 

error 

initializing GEMPLT . 
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4.21 PROGRAM GPCLOSE 


GPCLOSE closes the current graphics output file. 

INPUT PARAMETERS 
None 


PROGRAM DESCRIPTION 

This program closes the graphics output file, if the 
device produces an intermediate plot file. The plot 
commands may be sent to the device using GPOUT. It 
necessary to call GPCLOSE before calling GPOUT. 


ERROR MESSAGES 


[GPCLOSE -1] 
[GPCLOSE -3] 


Fatal error initializing TAE. 
Fatal error initializing GEMPLT . 


current 
file 
s not 
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4.22 PROGRAM GPCOLOR 


GPCOLOR changes the colors on a color device. 


INPUT PARAMETERS 

COLORS *Co lor list 

CTYPE *Co lor type 

CCCMP ‘Color components or name 

DEVICE ‘Graphics device 


PROGRAM DESCRIPTION 

This program changes the color components of color numbers 
on color devices where the components can be set. 

COLORS is the color number. CTYPE gives the method to use 
when the components are changed. If CTYPE -IN IT, then all 
the colors are set to default values. In this case, the values 
of COLORS and CCCMP are ignored. If CTYPE = RGS , the red 
green and blue components must be specified. If CTYPE = IILS , 
the hue, lightness and saturation are required. If 
CTYPE - NAME, the color components will be read from the 
GEMPAK color tables. Note that only the first three letters 
of the n ame will be used. 

CCCMP gives the color components in one of the following ways: 

NAME: The color is given by name, where only the first 

3 letters are used . 

RGB: The red, green and blue components are given; the 

values mn st be separated wi t h s emi colons . 

HLS : The hue, lightness and saturation are given; the 

values must be separated with semicolons. 

The RGB and HLS numbers must be in the range 0 to 1. 


The list of colors in the GEMPAK table with the corresponding 
red, green, blue, hue, lightness, and saturation values follows 


RED 

RED 

GREEN 

GRE 

BLUE 

BLU 

YELLCW 

YEL 


1.000 .000 .000 

.000 1.000 .000 

.000 .000 1.000 

1.000 1.000 .000 


.000 

.333 

1 .000 

.333 

. 333 

1 . 000 

.667 

.333 

1 .000 

. 167 

.667 

1 . 000 
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CYAN CYA 

MAGENTA MAG 

WHITE WHI 

BLACK BLA 

I 

! Of f wh i t e s 
GRAY GRA 

VANILLA VAN 

I 

I Oranges 

ORANGE ORA 

DKORANGE DKO 

I 

! Browns 

GBRCWN GBR 

BRCWN BRO 

BEIGE BE I 

SAND SAN 

MUD MUD 

! 

! Blues 

AQUA AQU 

SKY SKY 

LTBLUE LTB 

NAVY NAV 

DKBLUE DKB 

INDIGO IND 

! 

1 Greens 

LTGREEN LTG 

AVOCADO AVO 

DKGREEN DKG 

1 

1 Yell ows 

TAN TAN 

GOLD GOL 

BLOND BLO 

LEMON LEM 

I 

i Reds 

PINK PIN 

DKPINK DKP 

MAROON MAR 

APRICOT APR 

VIOLET VIO 

PURPLE PUR 


.000 

1.000 

1 .000 

1 .000 

.000 

1 .000 

1 .000 

1 .000 

1.000 

.000 

.000 

.000 


. 500 

.500 

.500 

1 .000 

.900 

.750 


1 .000 

.500 

.000 

. 863 

.471 

.000 


.500 

.281 

.281 

.500 

.294 

. 196 

1.000 

.672 

.500 

. 851 

.600 

.298 

.395 

. 100 

. 160 


.000 

. 718 

.718 

.000 

.640 

1 .000 

.000 

.500 

1.000 

.000 

.000 

.500 

.000 

.000 

.750 

. 578 

.000 

1 .000 


.600 

1.000 

.000 

. 250 

.500 

.125 

.000 

.562 

.000 


. 750 

.500 

.000 

1.000 

.750 

.000 

. 900 

. 820 

. 390 

.900 

1 .000 

.250 


1 .000 

.600 

.600 

1 .000 

.500 

. 500 

.578 

.000 

.000 

1 .000 

. 700 

. 375 

.625 

.250 

.625 

. 500 

.000 

.500 


.500 

.667 

1 .000 

. 833 

.667 

1.000 

.000 

1.000 

.000 

.000 

.000 

.000 


.000 

.500 

.000 

. 102 

. 883 

. 320 


.083 

.500 

.816 

.092 

.445 

.710 


.000 

.354 

.216 

.051 

.330 

. 269 

.055 

.724 

.464 

.092 

.583 

.456 

.970 

.218 

. 308 


.500 

.479 

. 678 

.558 

.547 

. 831 

.583 

.500 

.816 

.667 

. 167 

.632 

.667 

.250 

.816 

.764 

.526 

. 821 


.232 

. 533 

. 824 

.280 

.292 

. 343 

.333 

.187 

.679 


.113 

.416 

.632 

. 129 

.583 

. 862 

. 143 

. 703 

.490 

.186 

.717 

. 738 


.000 

.733 

.426 

.000 

.667 

.500 

.000 

.193 

.691 

.000 

.692 

.553 

. 833 

.500 

. 354 

. 833 

.333 

. 500 


EXAMPLES 

1. Set color number 4 to violet. 
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COLORS - 4 

CCOMP - violet 
CTYPE - name 

2. Set color number 2 to yellow by specifying the red, green, 
and blue components. 

COLORS - 2 

CCOMP - 1 ; 1 ; 0 

CTYPE - rgb 

3. Set color number 3 to DKPINK by specifying the hue, 
lightness, and saturation components. 


COLORS - 3 

CCOMP - 0 ; . 667 ; .5 

CTYPE - his 


ERROR MESSAGES 


[ GPCOLOR 

-1] 

[ GPCOLOR 

-2] 

[GPCOLOR 

-3] 

[GPCOLOR 

-4] 

[GPCOLOR 

-5] 

[GPCOLOR 

-6] 

[GPCOLOR 

-7] 


Fatal error initializing TAE. 

Fatal error reading TAE parameters. 

Fatal error initializing GEMPLT . 

... is an invalid input for CTYPE. 

Invalid input for COLORS. Enter one value. 
... has fewer than 3 color components. 

... has more than 3 color c omp o n e n t s . 
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4.23 PROGRAM GPEND 

GPEND terminates the GEMPLT subprocesses. 

INPUT PARAMETERS 
None 


PROGRAM DESCRIPTION 


This program terminates the GEMPLT plotting package 
be used to stop the GEMPLT subprocesses. 


ERROR MESSAGES 

[GPEND -1] 
[GPEND -3] 


Fatal error 
Fatal error 


initializing TAE. 
initializing GEMPLT . 


It can 
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4.24 PROGRAM GPMAP 


GPMAP defines a graphics area and draws a map and 1 a t i tude / 1 ong i t ud e 
lines. 


INPUT PARAMETERS 

MA.P ‘Map co 1 o r /da sh/wi d th 

DEVICE ‘Graphics device 

GAREA ‘Graphics area 

PRO J ‘Map pr o j ec t i on/ang 1 e s /marg ins 

PANEL ‘Panel 1 oca t i on/ c o 1 o r /da sh/wid th 

TITLE ‘Title col or/ 1 ine/ t i 1 1 e 

LATLON ‘Line color/dash/width/label/inc 

CLEAR ‘Clear screen flag 


PROGRAM DESCRIPTION 

This program defines a graphics area and device and draws 
a map and/or 1 a t i tude/ 1 ong i tude lines. 


EXAMPLES 

1. Define a graphics area centered on WV on the SG device. 

Set the projection to North Polar Stereographic. Draw 
a map with solid lines in color 3. Draw latitude/longitude 
lines every 5 degrees in color 2 using a dotted dash pattern. 
Do not write a title on the map. 

GAREA - wv 

PROJ - nps 

MAP - 3 

LATLON - 2/10///5;5 

PANEL - 0 

TITLE - 0 

CLEAR - yes 

DEVICE - sg 


ERROR MESSAGES 


[GPMAP 

-1] 

Fatal 

error 

initializing TAE . 

[GPMAP 

-2] 

Fatal 

error 

reading TAE par ame ters 

[GPMAP 

-3] 

Fatal 

error 

initializing GEMPLT . 
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4.25 PROGRAM GPOUT 

GPOUT sends GEMPLT graphics files to the output device. 

INPUT PARAMETERS 
None 

PROGRAM DESCRIPTION 

This program sends intermediate plot files to the graphics 
device. If a plot file is open, it will be closed and sent 
to the plotter. Otherwise, the most recent file for this 
device will be plotted. 

ERROR MESSAGES 

[GPOUT -1] Fatal error initializing TAE. 

I GPOUT -3] Fatal error initializing GEMPLT. 


4-50 



GEMPAK PROGRAMS 


4.26 PROGRAM GPREST 


GPREST restores the graphics plot, saved by GPSAVE, to a graphics 
device having readback capability. 


INPUT PARAMETERS 

SAVFIL *Gr aph i c s save file name 
DEVICE *Gr aph i c s device 


PROGRAM DESCRIPTION 

This program restores the graphics plot to a graphics device 
with readback capability. The file must have been created 
using GPSAVE . 


EXAMPLES 

1. Restore a surface map from a specified disk and directory. 

SAVF I L - s f cma p . s a v 
DEVICE - i i 

2. Restore a skew-t diagram from the default disk and 
directory. 

SAVFIL - skewt.sav 


ERROR MESSAGES 

[GPREST -1] 
[GPREST -2] 
[GPREST -3] 
[GPREST -4] 


Fatal error initializing TAE. 

Fatal error in reading the TAE par ame t e r s . 
Fatal error initializing GEMPLT . 

Error restoring graphics from a file. 
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4.27 PROGRAM GPSAVE 


GPSAVE saves the graphics plot currently displayed to a disk file. 


INPUT PARAMETERS 

SAVFIL *Graphics save file name 

DEVICE *Graphi c s device 


PROGRAM DESCRIPTION 

This program saves the graphics plot displayed on DEVICE to a 
disk file. The graphics information saved can be restored 
using GPREST. These programs will only work on devices having 
a readback capability. 6 


EXAMPLES 

1. Save a surface map to a specified disk and directory. 
SAVF I L — s f cma p . s a v 

2. Save a skew-t diagram to the default disk and directory. 
SAVF I L - s k ewt . s a v 


ERROR MESSAGES 

Fatal error initializing TAE. 

Fatal error in reading TAE parameters. 
Fatal error initializing GEMPLT. 

Error saving graphics to a file. 


[GPSAVE -1] 
[GPSAVE -2] 
[GPSAVE -3] 
[GPSAVE -4] 
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4.2 8 PROGRAM OABOX 


OABOX draws a box around an objective analysis region. 


INPUT PARAMETERS 

LINE *Line c o 1 o r /da sh/wi d th/ 1 abe 1 

DEVICE ^Graphics device 

REGION *Type of region 

GAREA *Graphics area 

PRO J *Map projection/angles/margins 

GDFILE *Gr i d file 


PROGRAM DESCRIPTION 

This program draws a box around a region defined for the 
objective analysis programs. REGION can be specified as 
GRID, DATA, or EXTEND. The default is GRID. 

The GRID region is the area used for the objective analysis. 
The DATA region is the area over which data for the 
analysis is extracted. The EXTEND region is the grid 
area extended in each direction for the first pass analysis. 


EXAMPLES 

1. Draw a solid line around the grid area in color 1 using 
a solid line of width 1. 

LINE - 1 

REGION - grid 

GDF I LE - s amp 1 e . g r d 

DEVICE » vt 

GAREA = us 

PROJ = me r 


2. Draw a line around the data region using color 2, line type 
3, and line width 5. 

LINE - 2/3/S 

REGION = data 


ERROR MESSAGES 
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[OABOX 

-1] 

[OABOX 

-2] 

[OABOX 

-3] 

[OABOX 

-4] 

[OABOX 

-5] 


Fatal error initializing TAE. 

Fatal error reading TAE par ame t e r s . 

Fatal error initializing GEMPLT. 

No box will be drawn since the color is 0. 
Invalid input for REGION. 
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4.29 PROGRAM OABSFC 


OABSFC performs a Barnes objective analysis on surface data. 
INPUT PARAMETERS 


SFFILE 

♦Surface data file 

GDFILE 

♦Grid file 

SFPARM 

♦Surface parameter list 

DATTIM 

♦Da t e / t ime 

DTAAREA 

♦Data area for OA 

GUESS 

♦Guess f i 1 e * t ime 

GAMMA 

♦Convergence parameter 

SEARCH 

♦Search radius 

NPASS 

♦Number of passes 


PROGRAM DESCRIPTION 

OABSFC performs a Barnes objective analysis on surface 
data. Multiple parameters may be analyzed at the same 
time. If more than one time is entered, the analyses 
will be performed sequentially. 

The input surface file is specified in SFFILE. Two files 
may be entered by separating the names with a plus, +. 

The projection and area for the grid are read from the 
navigation block in the output grid file. The extend grid 
area, the station spacing and the data subset area are read 
from the analysis block. The extend area is used to define 
a larger grid for the first pass, allowing data to be 
interpolated back to stations outside the grid area. The 
average station spacing is used to compute the weighting 
functions. The data subset area specifies the data which 
will be used in the analysis. If DTAAREA is blank or set 
to DATA, the data area defined in the analysis block will 
be used. This is the recommended procedure. However, 
DTAAREA may be specified in the same way as the AREA 
parameter used in other programs. This is especially useful 
to eliminate stations which are known to have bad data. 

A first-guess field can be used in the analysis. The name 
of the grid file containing the guess field mu st be entered 
in GUESS. In addition, the time to be used to extract the 
guess field from the file must be entered as the second 
value in GUESS following a *. If a guess field is entered, 
the navigation info rma tion in the guess grid file mu s t be 
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identical to that in the output grid file. When a guess 
field is used, it is inserted into the grid as the zeroth 
pass. 


Each pass of the analysis interpolates data from stations to 
grid points using the weighting function: 

WTFUNC - EXP [ -DIST ** 2 / WEIGHT ] 

where: 


DIST ** 2 - [ ( lat (grid) - lat (stn) ) 
( Ion (grid) - Ion (stn) ) ** 2 
COSLSQ - COS ( lat (grid) ) ** 2 
WEIGHT - 5.051457 * ( DELTAN * 2 / PI ) 
DELTAN — Station spacing read fr om grid 


** 2 + 

* coslsq (grid) ] 
** 2 

file analys is block 


GAMMA, the convergence parameter, is a multiplier for the 
weight for passes after the first pass. GAMMA must be within 
the range 0-1. Any value outside this range will default to 
a value of .3. If GAMMA is 0, the number of passes will be set 
to 1. The recommended value for GAMMA is .3. 


SEARCH is used to control the search radius, which is the 
maximum distance that a station may be from a grid point to be 
used in the analysis for that point. The search radius will be 
£v£ ! fl ^l! Utions whose weighting factor would be less than 
EXP (-SEARCH) will not be used. SEARCH must be in the range 
1 - 50. If it is not, a default value of 20 wi 1 1 be used. If 
a very small value is used, many grid points will not have 3 
stations within the search area and will be set to the missing 
data value. 6 


NPASS controls the number of passes. Valid values are in the 
range 1-5. Note that two passes are STRONGLY RECOMMENDED . 


EXAMPLES 


1. Analyze temperature, dewpoint and the wind components 
from the latest time in the surface file, SURFDAT. Add 
the computed grids to the file, TEST.GRD. 


SFFILE 

GDFILE 

SFPARM 

DATTIM 

DTAAREA 

GUESS 

GAMMA 


su r f da t 
test .grd 

tmpf ; dwpf ; uwnd , vwnd 

1 a 8 t 

data 

.3 
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SEARCH - 10 

NPASS - 2 

2. Repeat the dewpoint analysis setting the convergence 
parameter to .25 and eliminating the data at stations DCA 
and BUF . 

SFPARM - dwpf 
DTAAREA - data/-®dca;buf 
GAMMA - .25 

3. Repeat the wind analysis using the 24-hour forecast from 
the model data stored in model. ngm. 

SFPARM = uwnd ; vwnd 
GUESS - mod e 1 . ngm* f 24 


ERROR MESSAGES 

[OABSFC +9] 
[OABSFC +8] 
[OABSFC +7] 
[OABSFC +6] 
[OABSFC +5] 
[OABSFC +4] 
[OABSFC +3] 
[OABSFC +2] 
[OABSFC +1] 
[OABSFC -1] 
[OABSFC -2] 
[OABSFC -3] 
[OABSFC -4] 
[OABSFC -5] 
[OABSFC -6] 
[OABSFC -7] 
[OABSFC -8] 
[OABSFC -9] 
[OABSFC -10] 
[OABSFC -11] 


No data from ... will be used. 

Parameter ... is a character parameter. 
Parameter ... cannot be computed, 

WARNING: Area is not DATA area in file. 

WARNING: The recommended number of passes is 2 

... is invalid for NPASS. It is set to 2. 

... is invalid for search. It is set to 20. 
WARNING: GAMMA is 0. There will be only 1 pas 

... is invalid for GANMA. It is set to .3. 
Fatal error initializing TAE . 

Fatal error reading TAE parameters. 

Fatal error initializing GEMPLT . 

Grid size is too large. 

Not enough buffer space; reduce number of grids 
There are too many stations in data subset area 
There are too few stations in data subset area. 
There are no times in the surface file. 

No valid parameters have been entered. 

Surface file ... could not be opened. 

Time cannot be found in .... 
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4.30 PROGRAM OABSND 


OABSND performs a Barnes objective analysis on upper air data. 
INPUT PARAMETERS 


SNFILE 

♦Sounding data file 

GDF I LE 

*Gr id file 

SNPARM 

♦Sounding parameter list 

STNDEX 

♦Stability indices 

LEVELS 

♦Ve rtical levels 

VCOORD 

♦Vertical coordinate type 

DATTIM 

♦Da t e/ t ime 

DTAAREA 

♦Data area for OA 

GUESS 

♦Gue ss f i 1 e ♦ t ime 

GAMMA 

♦Conve rgence par ame ter 

SEARCH 

♦Search radius 

NPASS 

♦Number of passes 


PROGRAM DESCRIPTION 

OABSND performs a Barnes objective analysis on upper air 
data. Multiple parameters and levels may be analyzed at the 
same time. If more than one time is entered, the analyses 
will be performed sequentially. 

The input sounding file is specified in SNFILE . Two files may 
be entered by separating the names with a plus, +. 

The projection and area for the grid are read from the 
navigation block in the output grid file. The extend grid 
area, the station spacing and the data subset area are read 
from the analysis block. The extend area is used to define 
a larger grid for the first pass, allowing data to be 
interpolated back to stations outside the grid area. 

The average station spacing is used to compute the weighting 
functions. The data subset area specifies the data which will 
be used in the analysis. If DTAAREA is blank or set to DATA, 
the data area defined in the analysis block will be used. 

This is the recommended procedure. However, DTAAREA may be 
specified in the same way as the AREA parameter used in other 
programs. This is especially useful to eliminate stations which 
are known to have bad data. 

A first-guess field can be used in the analyses. The name 
of the grid file containing the guess field must be entered 
in GUESS. In addition, the time to be used to extract the 
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guess field from the file must be entered as the second 
value in GUESS following a *. If a guess field is entered, 
the navigation info rma tion in the guess grid file mu s t be 
identical to the that in the output grid file. 

Each pass of the analysis interpolates data from stations to 
grid points using the weighting function: 

WTFUNC - EXP [ -DIST ** 2 / WEIGHT ] 

wh ere: 

DIST ** 2 - [ ( 1 a t (grid) - lat (stn) ) ** 2 + 

( Ion (grid) - Ion (stn) ) ** 2 * coslsq (grid) ] 
COSLSQ - COS ( lat (grid) ) ** 2 
WEIGHT - 5.051457 * ( DELTAN * 2 / PI ) ** 2 

DELTAN - Station spacing read from grid file analysis block 


GAMvlA, the convergence parameter, is a multiplier for the 
weight for passes after the first pass. GANWA. mu s t be within 
the range 0-1. Any value outside this range will default to 
a value of .3. If GAMMA is 0, the number of passes will be set 
to 1. The recommended value for GANMA is .3. 

SEARCH is used to control the search radius, which is the 
maximum distance that a station may be from a grid point to be 
used in the analysis for that point. The search radius will be 
set so that stations whose weighting factor would be less than 
EXP (-SEARCH) will not be used. SEARCH must be in the range 
1 - 50. If it is not, a default value of 20 will be used. If 
a very small value is used, many grid points will not have 3 
stations within the search area and will be set to the missing 
data value. 

NPASS controls the number of passes. Valid values are in the 
range 1-5. Note that two passes are STRONGLY RECOMMENDED. 


EXAMPLES 


1. Analyze temperature, dewpoint and the wind components 

at 850 and 500 mb from the file, SNDDAT . Add the output 
grids to the file, TEST.GRD. Use default values in the 
analyses . 


SNFILE 

GDFILE 

SNPARM 

STNDEX 

LEVELS 


sndda t 
t e s t . g r d 

tmp c ; dwp c ; uwnd , vwnd 
850;500 
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VCOORD — pres 
DATTIM - last 
DTAAREA - data 
GUESS 

GAMMA = . 3 

SEARCH - 10 

NPASS - 2 

2. Do the surface analysis for the same parameters. Include 
a pressure analysis also. 

SNPARM - tmp c ; dwpc ; uwnd ; vwnd ; p r e s 
LEVELS = 0 

3. Repeat the wind analysis using the 24-hour forecast from 
the model data stored in model. ngm. 

SNPARM - uwnd ; vwnd 
GUESS — mod e 1 . ngm* f 24 


ERROR MESSAGES 

[OABSND +10] 
[ OABSND +9] 
[OABSND +8] 
[OABSND +7] 
[OABSND +6] 
[OABSND +5] 
[OABSND +4] 
[OABSND +3] 
[OABSND +2] 
[OABSND +1] 
[OABSND -1] 
[OABSND -2] 
[OABSND -3] 
[OABSND -4] 
[OABSND -5] 
[OABSND -6] 
[OABSND -7] 
[OABSND -8] 
[OABSND -9] 
[OABSND -10] 
[OABSND -11] 
[OABSND -12] 


No data from ... will be used. 

Station par ame ter ... cannot be c ompu ted. 
Parameter ... is a character parameter. 
Parameter ... cannot be computed. 

WARNING: Area is not DATA area in file. 

WARNING: The recommended number of passes is 2 

... is invalid for NPASS. It is set to 2. 

... is invalid for search. It is set to 20. 
WARNING: GAMMA is 0. There will be only 1 pas 

... is invalid for GAMMA. It is set to .3 . 
Fatal error initializing TAE. 

Fatal error reading TAE par ame t e r s . 

Fatal error initializing GEMPLT. 

Grid size is too large. 

Not enough buffer space; reduce number of grids 
There are too many stations in data subset area 
There are too few stations in data subset area. 
There are no times in the surface file. 

No valid parameters have been entered. 

LEVELS cannot contain a range. 

Sounding file could not be opened. 

Time cannot be found in . . . . 
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4.31 PROGRAM OAGRID 


OAGRID creates a GEMPAK grid file which can be used in a Barnes 
objective analysis program. 

INPUT PARAMETERS 


GDFILE 

♦Grid file 

DELTAN 

♦Station spac ing 

DELTAX 

♦X spacing 

DELTAY 

♦Y spacing 

GRDAREA 

♦Area covered by grid 

EXTEND 

♦Points to extend grid 

DTAAREA 

♦Data area for OA 

SOURCE 

♦Data source (SN or SF) 

SNFILE 

♦Sounding data file 

SFFILE 

♦Surface data file 

SNPARM 

♦Sounding parameter list 

SFPARM 

♦Surface parameter list 

DATTIM 

♦Da t e / t ime 

LEVELS 

♦Ve rtical levels 

MAXGRD 

♦Maximum number of grids 


PROGRAM DESCRIPTION 

Th is p r o g r am all ows the user to create a GEMPAK grid file 
which contains the information required to perform a Barnes 
objective analysis. 

The output grids will be evenly spaced 1 a t i tude / 1 ong i tud e 
(CED) grids. Three areas used by the objective analysis 
p r o g r ams are defined in this progr am . GRDAREA defines a 
region for the output grid. The upper right corner 
specified will be moved toward the lower left in order 
to align it on a grid point. 

The second area, defined by EXTEND, is used to extend the 
grid area outward by some number of grid points. This area 
is used as a first-pass grid area in the objective analysis. 
Th e default EXTEND values are 2, 2, 2, 2. 

DTAAREA defines the area over which station data will be 
input to the analysis. Only data within the EXTEND area are 
used for the second pass. If a value for DTAAREA is not 
specified, the EXTEND area will be used. 

If values for the station spacing and grid spacings, DELTAN , 
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DELTAX, and DELTAY, are all specified by the user, they are 
stored in the grid file for use by the analysis programs. If 
any of these numbers is 0, a suggested station spacing is 
computed using the station data from the file specified. This 
station spacing will be used to compute values for DELTAX and 
DELTAY. The computed station spacing is the average of the 
average minimum station spacing and the uniform station 
spacing. These station spacing values are computed using the 
stations reporting data for the first DATTIM. The average 
minimum station spacing is the average of the distances from 
each station to its closest station. The uniform spacing is 
the spacing that would be found between stations if they were 
evenly spaced over the data area. 

GDCF IL can also be used to create grid files. OAGR ID creates 
only CED grids, but values for the grid spacing and station 
spacing are estimated from the input surface or upper air file 
GDCF I L creates grid files for grids in any projection, but 
the grid spacing and station spacing must be input directly. 


EXAMPLES 


1. Create a grid file called TEST.GRD. Use the surface 

temperature data at the latest time from stations in the 
surface file SURFDAT to compute the station spacing. The 
grid area is WV; use defaults for the remaining parameters 
Create a file wh ich can contain up to 100 grids. 


2 . 


GDF I LE 

- 

test . g r d 

DELTAN 

- 

0 . 

DELTAX 

- 

0 . 

DELTAY 

am 

0 . 

GRDAREA 

- 

wv 

EXTEND 

- 

2 ; 2 ; 2 ; 2 

DTAAREA 

- 


SOURCE 

am 

sf 

SFFILE 

- 

surf da t 

SFPARM 

- 

tmpf 

DATTIM 

■» 

last 

LEVELS 

- 

0 

MAXGRD 


100 


Create a grid file called GRID.GRD with a station spacing 
of 5, x spacing of 3 and y spacing of 4. The grid area 
is over PA; the extended grid area should include 1 grid 
point in each direction; the data area is PA- . 

GDF I LE - grid.grd 

DELTAN - 5 . 

DELTAX - 3 . 
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DELTAY - 4 . 

GRDAREA - pa 
EXTEND - 1 ; 1 ; 1 ; l 

DTAAREA - pa- 


ERROR MESSAGES 

[OAGRID -1] 
[OAGRID -2] 
[OAGRID -3] 
[OAGRID -4] 
[OAGRID -5] 
[OAGRID -6] 
[OAGRID -7] 
[OAGRID -8] 
[OAGRID -9] 
[OAGRID -10] 
[OAGRID -11] 
[OAGRID -12] 
[OAGRID -13] 


Fatal error initializing TAE. 

Fatal error reading TAE par ame t e r s . 
... is invalid for GRDAREA. 

... is invalid for DTAAREA. 

No data file name specified. 

Parameter input is invalid. 

Parameter ... cannot be calculated. 
Invalid value for DELTAX or DELTAY. 
Too few stations to calculate DELTAN . 
Source must be set to SN or SF . 
Station data file is invalid. 

Invalid time requested. 

Invalid level requested. 
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4.32 PROGRAM SFCFIL 


SFCFIL creates a new GEMPAK surface file. 


INPUT PARAMETERS 

SFOUTF ‘Output surface file 

SFPRMF ‘Surface parameter packing file 

STNFIL ‘Station information file 

SHIPFL ‘Ship data file flag 

TIMSTN ‘Time s / add i t i ona 1 stations 


PROGRAM DESCRIPTION 

SFCFIL creates a GEMPAK surface file. The file may be a 
standard file or a ship format file. 

If SHIPFL is set to YES, a ship format file is created. In 
this case, both station location and time are stored with 
each data entry. This capability is useful wh en the station 
location varies in time, such as for moving ships, aircraft 
or free-floating buoys. 

If a standard file is to be created, the maximum number of 
times to be included in the file must be entered as the first 
value in TIMSTN. If a ship format file is being created, the 
maximum number of entries in the file is given by the first 
value in TIMSTN. 

If STNFIL is not blank, information about all the stations in 
STNFIL will be added to the data set. Space will be left in the 
file for the additional number of stations specified as the 
second parameter in TIMSTN. Note that an error will result if 
STNFIL is blank and TIMSTN does not request more stations. 

SFPRMF contains information about the parameters to be included 
in the file. SFPRMF may be either a list of parameters or the 
n ame of a packing file. If a list is entered, the par ame t e r s 
must be separated with semicolons; packing information may also 
be included after a slash with the minimum and maximum values 
and the resolution separated by dashes. For example, to 
include temperature and dewpoint in a file without packing, 
SFPRMF may be entered as : 

SFPRMF - TMPC ; EWPC 

To pack the data, use; 
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SFPRMF = TMPC/ -127-63- . 1 ; LWPC/ -63-63-1 

SFPRMF may also contain the name of a packing file. A packing 
file for data which is not to be packed contains a list of 
parameters with one parameter per line. In a file for packed 
data, each line must include the parameter name, the minimum 
and maximum data values and the resolution, all separated with 
spaces. The default packing file for surface data is 
GEMTABL : SFPACK . TBL . 


EXAMPLES 

1. Create a surface file called SURF . DAT with a maximum of 15 
times using the default station and parameter files. Leave 
room in the file for 100 stations in addition to the 
stations in STNFIL. Use the GEMPAK standard packing file 
and table file. 

SFOUTF - surf.dat 

SFPRMF - g emt ab 1 : s f pack . t bl 

STNFIL - g emt ab 1 : s t a t i ons . t bl 

SHIPFL - NO 

TIMSTN - 15/100 

2. Create a ship format file which can contain a maximum of 
1000 ship reports. Use a locally developed parameter file 
named SHIP.PRM for packing. 

SFOUTF - ship.dat 
SFPRMF = ship.prm 
STNFIL - 
SHIPFL - YES 
TIMSTN - 1000 


ERROR MESSAGES 

[SFCFIL +2] WARNING! ADDSTN was negative -- set to 0. 
[SFCFIL +lj Cannot add stn to a ship file; STNFIL ignored. 

[SFCFIL -l] Fatal error initializing TAE. 

[SFCFIL -2] Fatal error reading TAE parameters. 

[SFCFIL -3] Error opening station file ... . 

[SFCFIL -4] File does not include room for any stations. 

[SFCFIL -5] File does not include room for any times. 

[SFCFIL -6] SFPRMF is incorrectly specified. 

[SFCFIL -7] The output file name is blank. 
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4.33 PROGRAM SFCNVRT 


SFCNVRT converts a GEMPAK3 surface file to a GEMPAK 5 surface file. 


INPUT PARAMETERS 

OLDSFC *GEMPAK3 surface file 

SFOUTF ^Output surface file 

DATTIM *Date/time 


PROGRAM DESCRIPTION 

SFCNVRT creates a new GEMPAK surface data file and adds all 
the data at the requested times to the new file. All the 
stations in the original file will be included in the output 
file. No AREA subset is available. Since this program will 
create a new file each time that it is run, all the times 
to be converted must be specified in DATTIM. This program 
does not include a dynamic tutor, so new values cannot be 
input once SFCFIL is run. 

There is no difference between GEMPAK4 and GEMPAK5 surface 
files. This program is used to convert GEMPAK3 files since 
later GEMPAK programs cannot be run with the GEMPAK3 files. 


ERROR MESSAGES 


[ SFCNVRT 

-1] 

Fatal 

[SFCNVRT 

-2] 

Fatal 

[ SFCNVRT 

-3] 

Error 

[ SFCNVRT 

-4] 

Error 

[SFCNVRT 

-5] 

Error 

[ SFCNVRT 

-6] 

Error 


error initializing TAE. 
error reading TAE parameters, 
opening old file .... 
creating new file .... 
reading times from old file, 
finding times to process. 
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4 . 34 PROGRAM SFDELT 


SFDELT deletes data from a surface data file. 


INPUT PARAMETERS 

SFFILE ‘Surface data file 

DATTIM ‘Date/ time 

AREA “Data area 


PROGRAM DESCRIPTION 

SFDELT deletes data from a surface data file. The stations 
to be deleted may be specified in AREA; the times to be 
deleted are given in DATT IM . 

If AREA is set to DSET or ALL, all the data present at the 
times in DATTIM will be deleted, along with the headers for 
those t ime s . 

If AREA specifies a list of stations to be deleted, only the 
data for those stations at the specified times will be deleted. 
The station and time information in the file will remain 
unchanged . 


EXAMPLES 

1. Delete all the data at time 850722/1200. 

SFFILE - 85 jul . sf c 

DATTIM - 850722/1200 

AREA - d s e t 

2. Delete data at BWI for the most recent time in the file. 

SFFILE - realtime. sfc 
DATTIM - last 
AREA - ®bwi 


ERROR MESSAGES 

[SFDELT -l] Fatal error initializing TAE. 
[SFDELT -2] Fatal error reading TAE parameters. 

[SFDELT -3] Error deleting data at time .... 
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4.35 PROGRAM SFEDIT 


SFEDIT adds or changes data in a surface file using a sequential 
text file. 


INPUT PARAMETERS 

SFEFIL ‘Surface edit file 

SFFILE ‘Surface data file 


PROGRAM DESCRIPTION 

SFEDIT adds information from a sequential text file to an 
existing GEMPAK surface file. The program may be used to add 
new data to a file or to change existing values. The output 
file mu st already exist. 

The data to be added must reside in a text file, SFEFIL. This 
file can be created by specifying F as the output device in 
GEMPAK program SFLIST or by generating the file with a text 
editor. A text editor may be used to change the file. 

The parameters to be edited must be specified at the beginning 
of the edit file, for example: PARM-TMPF ; EWPF . Only 

parameters stored in the output file can be edited. Parameters 
that have character values, such as WTHR, are stored as real 
numbers. Thus, the numeric weather code must be used. 

Station information, such as the latitude, SLAT, cannot be 
changed using this program. Rather, these changes can be 
made using the GEMPAK program, SFSTNS. 

The data follow the parameter list. Each data listing must 
include the time, followed by the station character or 
numeric identifier followed by the data. If data are 
missing, the current missing data value, -9999., must be 
used. Blank fields will not be recognized. 

If a station or time does not exist in the output file it will 
be added to the file. If stations are added, GEMPAK program 
SFSTNS should be used to add station latitudes, longitudes 
and e 1 e va t i on s . 


EXAMPLES 

1. Change values of TMPF and EWPF in the file SFCDATA.DAT. 
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SFEFIL - sflist.fil 
SFFILE - sfcdata.dat 

The file SFLIST.FIL follows. Note that the second line 
containing titles is ignored. 

PARM - TMPF ;EWPF 

STN YYMMDD/HHMM TMPF DWPF 

BWI 850426/0400 55.40 42.80 

SBY 850426/0400 62.60 55.40 

NHK 850426/0400 55.40 42.80 

The program replaces the temperature and dewpoint for 
stations BWI, SBY and NHK at the given time values in 
output file SFCDATA.DAT. 


ERROR MESSAGES 

[ SFEDIT +2] Station ... cannot be added to the file. 

[ SFEDIT -l] Fatal error initializing TAE. 

[SFEDIT -2] Fatal error reading TAE parameters. 

[SFEDIT -3] Edit file ... cannot be opened. 

[SFEDIT -5] PARM keyword not found. 

[SFEDIT -6] Too many parameters to edit. 

[SFEDIT -7] Time ... cannot be added to data set. 

[SFEDIT -8] Parameter ... is not in the data set. 

[SFEDIT -9] SFEDIT cannot continue with invalid parameters. 
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4.36 PROGRAM SFGRAM 


SFGRAM draws a meteogram for surface data. 
INPUT PARAMETERS 


SFFILE 

DATTIM 

STATION 

TRACE 1 

TRACE 2 

TRACE 3 

TRACE 4 

TRACE 5 

NTRACE 

TAXIS 

BORDER 

MARKER 

TITLE 

CLEAR 

DEVICE 

TEXT 


•Surface data file 
•Da t e / t ime 
•Stations 

•Parms/ colors/range/wi tness 
* Pa rms / colors/ range /wi tness 
•Parms/ colors/range /wi tness 
•Pa rms / colors/range /wi tness 
•Parms / colors/range /wi tness 
•Number of traces 
•Time 1 - t ime 2- t inc ; 1 bl ; g In ; tck 
•Background c o 1 o r /d a sh/wi d th 
•Marker co 1 o r / t ype/ s i z e /wi d th/hw 
•Title co 1 or / 1 i ne / t i 1 1 e 
•Clear screen flag 
•Graphics device 
•Text size/font/hw flag 


PROGRAM DESCRIPTION 


SFGRAM draws a time series plot for surface data. Up to 
five traces may be drawn in a single plot. NTRACE determines 
the number of traces to draw. If NTRACE is 5, the plot area 
will be divided into fifths; if NTRACE is 4, the plot area will 
be divided into fourths, etc. If STATION is a single station, 
all the traces will be for that station. If STATION is a 

st * ti( > ns ' TRACE1 will plot data from the first station, 
I RACE 2 will plot data from the second station, etc. 


Each TRACE parameter contains specifications for the parameters 
to be plotted in the corresponding trace. The format for 
each trace is: 


pa rms / colors/range/wi tnes ! pa rms /colors/range/wi tnes 

The parameters before the ! will be plotted on the left of the 
P lot ! thos e after the 1 will be plotted on the right. 

Any GEMPAK surface parameter may be entered. Real valued 
parameters will be drawn as a graph. Character valued 
parameters will be rotated 90 degrees and written on the plot. 

n addition to the standard surface parameters, the following 
are also valid: 6 
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BARB 

ARKW 

DARR 

WSYM 


wind barbs 
wind arrows 
directional arrows 
weather symbols 


For GUST and GUMS, the character, G, will be plotted. 


Up to four parameters may be plotted along each axis. The 
parameters must be separated using semicolons. Character 
and weather symbol data may only be plotted in positions 
1 , 2 or 3 , wh ere position 1 is at the bott om of the trace, 
position 2 is in the middle and position 3 is at the top. 

Each parameter name may be followed by an asterisk, the 
size or line dash type, a second asterisk, and the width. 
For example, WSYM* . 5 * 5 wi 1 1 draw weather symbols half the 
default size with a line width of 5. TMPF*3 will plot 
a t emp erature line using dash pattern 3. 

The colors for the parameters must also be separated using 
semicolons. If a single number is entered, all parameters 
are dr awn in that color. If a zero is entered, the current 
default color is used. 


The range specifies the scaling of the y axis, 
i s : 

start ;stop; incr erne n t 


The format 


Note that, in this program, the parts of range must be 
separated using semicolons. If no range is given, it 
is selected using the data values. 

Witness lines are specified inwitnes. These are horizontal 
dotted lines. A list of y values, separated by semicolons, 
may be entered. Alternatively, if the value of witnes is YES, 
a witness line will be centered on the plot. 


The time axis is specified in TAXIS as a minimum time, a 
maximum time, and a time increment separated with dashes 
If any or all of the parts of TAXIS is blank, reasonable 
values will be selected by the program. 


EXAMPLES 

1. Plot a meteogram for BWI using all the times in the file. 
Use the default time axis and background colors. Do not 
plot markers. TRACE 1 plots TMPC and DWPC in colors 2 
and 3 with a data range from -10 to 30 degrees and a 
witness line at 0 degrees. TRACE 2 will plot VSBY using a 
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solid line connecting points with labels on the left and 
PMSL using a dashed line with labels on the right. TRACE3 
plots directional arrows at the bottom of the plot along 
with the wind speed in knots and wind gusts. TRACE 4 
plots cloud cover codes at the bottom and weather symbols 
with size - .6 at the top. TRACES plots the ALTM. 


SFFILE 

DATTIM 

STATION 

TRACE 1 

TRACE 2 

TRACE 3 

TRACE 4 

TRACES 

NTRACE 

TAXIS 

BORDER 

MARKER 

TITLE 

CLEAR 

DEVICE 

TEXT 


nws $ d a t a :sfdata.sfc 

all 

abwi 

tmpC ;dwpc/2 ; 3/- 10 ; 30 ; 10/0 
vs by 1 pms 1*3 
d a r r ; skn t ; gu s t 
elds; ;wsym* . 6 
a 1 tm 
5 

1 

0 

1 

ye s 
qm 


ERROR MESSAGES 

[SFGRAM +3] 
[ SFGRAM - 1 ] 
[ SFGRAM - 2 ] 
[SFGRAM -3] 
[SFGRAM -4] 
[ SFGRAM - 5 ] 
[ SFGRAM - 6 j 
[SFGRAM -7] 
[ SFGRAM - 8 ] 
[ SFGRAM - 9 ] 
[SFGRAM -10] 
[SFGRAM -11] 
[SFGRAM -12] 
[SFGRAM -14] 
[SFGRAM -15] 
[SFGRAM -16] 


Hardware text cannot generally be rotated. 
Fatal error initializing TAE. 

Fatal error reading TAE parameters. 

Fatal error initializing GEMPLT. 

No stations have been entered. 

... cannot be plotted in position 4. 

The time range along the x axis is 0. 

The parameter ... cannot be computed. 

The graph coordinates are invalid. 

There is no data at station .... 

Station ... is invalid. 

There are no times in the file. 

The time ... is invalid. 

The file ... cannot be opened. 

There are no parameters specified. 

Error in specifying TAXIS. 
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4.37 PROGRAM SFL604 


SFL604 lists data from a GEMPAK surface file in a fixed format 


INPUT PARAMETERS 


SFFILE 

♦Surface data 

file 


AREA 

♦Data area 



DATTIM 

♦Da t e / t ime 



OUTPUT 

♦Output device 

(TPF ) 

SKPMIS 

♦Skip missing 

data 

flag 

IDNTYP 

♦STNM or STID 



SFPARM 

♦Surface parameter 

list 


PROGRAM DESCRIPTION 

SFL604 lists data from a surface file in a format designed to 
display AIRWAYS data. 

Data will be listed for the stations and times requested in 
AREA and DATTIM. If SKPMIS is YES, stations which are within 
the area requested, but wh ich have not reported data, will not 
be listed. If SKPMIS is set to NO, then only the date and 
station identifier will be listed for nonreporting stations. 

Stations may be identified in the output listing by character 
or numeric identifier. This is accomplished by setting IDNTYP 
to STID or STNM. 

This program will list predefined parameters only. However, 
stations may be selected using conditions in SFPARM. For 
example, if SFPARM - WTHR-S+ , then only those stations 
reporting heavy snow will be listed. The conditions specified 
in SFPARM can be for any computable parameter, not just 
those which will be listed. 


EXAMPLES 

1. List the data to a file for all the stations in Maryland 
at the most recent time. 


SFFILE - 

nws $ d a t a : n ews f c . d a t 

AREA 

@md 

DATTIM - 

last 

OUTPUT - 

f 

SKPMI S - 

ye s 
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IDNTYP = st id 
SFPARM - 

2. List the data at 12Z and 13Z at BWI and DCA. The following 
parameters must be changed from example 1: 

DATUM = / 1 2 - / 1 3 

AREA - Obwi ; d c a 

3. List all stations in Delaware, including nonreporting 
stations. Identify the stations with a station number. 

The following parameters must be changed from example 1: 

AREA - @de 
SKPMIS - no 
IDNTYP - stnm 

4. List all stations reporting temperatures greater than or 
equal to 75 degrees. The following parameters must be changed 
from ex amp lei: 

SFPARM — tmpf>75 


ERROR MESSAGES 

[ SFL604 

+ 1] 

[ SFL604 

-1] 

[ SFL604 

-2] 

[ SFL604 

-3] 

[ SFL604 

-4] 


Parameter . . . cannot be computed. 
Fatal error initializing TAE. 

Fatal error reading TAE par ame t e r s . 
No valid stations could be found. 

No parameters can be computed. 
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4.38 PROGRAM SFLIST 


SFLIST lists surface data from a GEMPAK surface data file 


INPUT PARAMETERS 

SFFILE ’Surface data file 

AREA *Da t a area 

DATTIM *Da t e / t ime 

SFPARM ’Surface parameter list 
OUTPUT ’Output device (TPF) 

IDNTYP ’ STNM or ST ID 


PROGRAM DESCRIPTION 

SFLIST lists any parameters which can be derived from the data 
in a surface data file. The stations and times to be included 
are specified in AREA and DATTIM. If no data are reported for 
a station, that station will not be listed. The listings will 
be grouped by time. 

Conditions can be specified for the parameters. The conditions 
are documented in the SFPARM documentation. 


EXAMPLES 


1. List the air and dewpoint temperatures in Fahrenheit and 

the mean sea-level pressure of the stations in Maryland for 
the latest t ime : 


SFFILE 

AREA 

DATTIM 

SFPARM 

OUTPUT 

IDNTYP 


nws $ d a t a : n ews f c . d a t 


@md 

last 

tmpf ; dwpf ; pms 1 


t 

s t i d 


2. List the same parameters for DCA, and ORD at 0900 and 
1000 GMT at the terminal and line printer. Identify 
stations by their station numbers. 


AREA 

DATTIM 

OUTPUT 

IDNTYP 


@d c a ; o r d 
09-10 

tp 

s tnm 
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3. 


List the temperature, dewpoint, weather and cloud 
for stations in Maryland which are reporting any 


AREA 

«md 

DATTIM - 

last 

SFPARM - 

tmpf ; dwpf ;wthr-s-/s/s+; elds 

OUTPUT - 

t 

IDNTYP - 

s t i d 


ERROR MESSAGES 

[SFLIST 

+ 1] 

[SFLIST 

-1] 

[SFLIST 

-2] 

[SFLIST 

-3] 

[SFLIST 

-4] 


Par ame ter ... is not c omp u t a b 1 e . 
Fatal error initializing TAE. 

Fatal error reading TAE parameters. 
No stations reporting data. 

No valid computable parameters. 


cove r 
snow. 
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4.39 PROGRAM SFMAP 


SFMAP plots surface station data on a map. 


INPUT PARAMETERS 

AREA *Da ta area 

GAREA ‘Graphics area 

SFPARM ‘Surface parameter list 

DATTIM ‘Date/ time 

SFFILE ‘Surface data file 

COLORS ‘Color list 

WIND "Wind symbo 1 / s i z/wd t h/ typ/hd s z 

MAP ‘Map col or /dash/width 

LATLON ‘Line co 1 o r /da sh/wi d th/ 1 abe 1 / i nc 
MARKER ‘Marker co 1 or / type / s i z e /wi d th/hw 
TITLE ‘Title co 1 o r / 1 i ne / t i t 1 e 

CLEAR ‘Clear screen flag 

PANEL ‘Panel 1 oca t i on/ c o 1 o r /da sh/wid th 

DEVICE ‘Graphics device 

PROJ ‘Map pr o j e c t i on/ ang 1 e s /ma r g ins 

FILTER ‘Filter data flag 

TEXT ‘Text size/font/hw flag 

SKPMIS ‘Skip missing data flag 


PROGRAM DESCRIPTION 

SFMAP plots data at station locations on a map. Any 
parameter that can be computed from the parameters in 
the data set can be displayed. 

The order of the input parameters determines their location 
on the plot. Parameters are separated by semicolons. A 
position may be skipped by entering two consecutive semicolons 
or entering the name SPAC or BLNK. A parameter will be 
plotted centered at the station if it is the last parameter 
in the list and the symbol @ is appended to the name. The 
following chart shows the location on the model of each 
position. 

7 

1 9 3 

2 + 4 

5 10 6 

8 

If FILTER is set to YES, the stations will be filtered so that 
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?v?ni a n Pln ? 8tatlons wiU not b « Plotted. If a position is 

d w !! n * two semicolons or if the parameter is BLNK, the 

ui!d ^ n ° a 3 locate the s Pace- The parameter SPAC may be 

Plots will S ni r ot I?* SPaCC With thC fUtCr option - Then later 
U P 6 SamC statlons aft er filtering, provided 

that the same parameter locations are specified. 

J;; d !;! 0ns ‘ a “ be specified for the parameters. The conditions 
are documented in the SFPARM documentation. Weather symbols 

: d p ottcd * size and width ° f th « 
usina * 1 d by appendin 8 the numbers to the parameter name 
using a as a separator. For example, WSYM*2*S will plot 
weather symbols with size of 2. and line width of 5 P 


EXAMPLES 


1 . 


Plot a polar stereographic map with stations in the area 
cha r » e on Ma * y 1 and wi th temperature, dewpoint, and the 

""nd 4 WC a a nS e f r if? de ^ UtCSt * im ^ Us * Positions 
mart i' ? fllter overlapping stations. Do not plot 

markers where data are missing. Plot - P 

us ing color 2 . 


wind barbs in knots 


AREA 


md 

GAREA 

as 

md 

SFPARM 

- 

tmpf ; dwpf ; ;wthr 

DATTIM 

- 

last 

SFFILE 

« 

nws $ d a t a : n ews f c . d a t 

COLORS 

- 

2 ; 3 ; 4 

WIND 

■a 

bk2 

MAP 


2 

MARKER 

_ 

3 

TITLE 

«x 

5 

SKPMI S 


ye s 

CLEAR 

- 

ye s 

DEVICE 

- 

c t 

PROJ 

- 

nps 

FILTER 

- 

ye s 


Plot a Mercator map of the Washington DC area with 
emperetures to the left of the station and mining ratio 
he right. Do not clear the screen prior to plotting 

?^lbof he “S US * ng dash pattern 3. Do not plot any wind 
ymb o 1 s . Th e changes fr om e x amp 1 e 1 a r e : 


AREA 

GAREA 

SFPARM 

WIND 

MAP 


d c a * * 
d c a * * 

; tmpf ; ;mi xr 
0 

1/3 
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CLEAR - no 

PROJ - me r 

3. Plot only weather s ymb ols centered at the station, 
the size of the weather symbol to half the default 
the line width to 5. 

MARKER - 0 
SFPARM -wsym®*.5*5 


ERROR MESSAGES 

[ SFMAP 

-1] 

[ SFMAP 

-2] 

[SFMAP 

-3] 

[SFMAP 

-4] 

[SFMAP 

-5] 


Fatal error initializing TAE. 

Fatal error reading TAE parameters. 
Fatal error initializing GEMPLT . 
Parameter ... is not calculable. 
Winds are not calculable. 


Change 
size and 
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4.40 PROGRAM SFMOD 


SFMOD moves selected surface data from an input surface file 
output surface file. 


to an 


INPUT PARAMETERS 

SFFILE ‘Surface data file 

SFOUTF ‘Output surface file 

DATTIM ‘Date/ time 

AREA ‘Data area 


PROGRAM DESCRIPTION 


data from an listing GEMPAK surface file, 

fill SFOUTF Wri rJ S thC data int ° another existing surface 
tile, SFOUTF. This program can be used to subset the 

original data set by time and/or stations. 


This program will not create 
program, SFCFIL, can be used 
parameters in the output file 
parameters in the input file, 
conversions will be done. 


a new surface file. The GEMPAK 
to create a surface file. If the 
are not the same as the 
the required parameter 


If a requested station or time is not in the 
Will be added to the file if there is room, 
times and stations that can be included in a 
by TIMSTN in SFCFIL. 


output file, it 
The number of 
file is specified 


EXAMPLES 


1 . 


TEST* SFP^* 1 t f0r /! 1 stations at the latest time in file 
TEST.SFC into a file called LAST.SFC. 


SFFILE — test.sfc 
SFOUTF — last.sfc 
DATTIM - last 
AREA - d s e t 


2 . 


Put the data for stations 
all times into the output 


in Maryland and Virginia for 
dataset. 


SFFILE ■■ test.sfc 
SFOUTF - last.sfc 
DATTIM - all 
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AREA - ®md/®va 


ERROR MESSAGES 

[ SFMOD 

+ 3] 

[ SFMOD 

+2] 

[SFMOD 

+ 1] 

[SFMOD 

-1] 

[SFMOD 

-2] 

[SFMOD 

-3] 

[SFMOD 

-4] 


Some stations were not added to file. 
Character parameter cannot be used: 
Parameter ... cannot be calculated. 
Fatal error initializing TAE . 

Fatal error reading TAE parameters. 
There are no parameters to be computed 
Time . . . cannot be added. 
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4.41 PROGRAM SFRAWD 


SFRAWD decodes raw surface data files and writes the output to a 
GEMPAK surface file. 


INPUT PARAMETERS 

SFFILE ‘Surface data file 

SFEFIL “Surface edit file 

DATTIM “Date/ time 

ADDSTN “Add station flag 

CTLFLG PARM NAME-CTLFLG TYPE- ( STRING , 4 ) DEFAULT-® $CTLFLG 

“Control character flag 
OUTPUT “Output device (TPF) 


PROGRAM DESCRIPTION 


decodes raw surface data stored in the text file, 
SFEFIL. The decoded data is stored in a GEMPAK surface 
ata file, SFFILE. This file must already exist. It can 

e created using SFCFIL. The following parameters will be 
decoded : 


PMSL ALT I TMPF DWPF SKNT DRCT GUST WNUM 
CHC1 CHC2 CHC3 VSBY 


If any of these parameters has not been included in the file 
an error message will be written. A packing file for these ’ 
parameters can be found in GEMTABL : SFPACK . TBL . 


The r aw 
received 
station 
records 
mu st beg 
r epor t s , 
CTLFLG i 
characte 
contain 


a s 


input file must contain either whole bulletins 
from the Domestic Data Service or a series of 
reports. If the data set contains bulletins, the 
in the file must be 80 bytes each, and each bulletin 
in with a <CNTL>A in column 1. If the input contains 
each report must be written on a single line, 
s a logical variable indicating whether control 
rs are present. If CTLFLG - YES, the input file must 
bulletins. If CTLFLG - NO, the file must have reports 


ADDSTN is a logical variable which determines whether 
not already in the file will be added to the file. 


s t a t i on s 


DATTIM is a range of times which will be 
The range should consist of a start time 
separated with a dash. If only one time 


included in the file, 
and an end t ime 
is included, 24 hours 
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of data will be saved. These times must be entered as full 
GEMPAK times. Note that if the input file has reports, rather 
than bulletins, the year, month and day will all be obtaine 
from the last time. Thus the times cannot span days. If the 
input file contains bulletins, the day will be obtained from 
the bulletin header and the year and month from the last time. 

OUTPUT specifies the output device for optional messages. 

These messages will include times and stations added to the 
file and stations which are not included in the file. If 
OUTPUT - N, no output will be written. 

Note that this program will execute slowly. This is to be 
expe c t ed . 


EXAMPLES 

1. Decode the data from a surface archive file containing 
bulletins from the Domestic Data Service. Do not add 
any stations to the file. Process only the times from 
05 UTC to 12 UTC . Write the output to a file called 
ma r22.dat. 

SFFILE - ma r22 . sf c 

SFEFIL = ma r 2 2 . surf 

DATTIM = 900322/0500-900322/1200 

CTLFLG - yes 

ADDSTN = no 

OUTPUT - f /ma r22.dat 

2. Decode the reports from BWI which were entered by hand 
in a text file. Do not write any optional output. 

SFFILE - bwi 22 . s f c 

SFEFIL - bwi 22 . s u r f 

DATTIM - 900322/0000-900322/2300 

CTLFLG - no 

ADDSTN - no 

OUTPUT = n 


ERROR MESSAGES 

[SFRAWD -1] Fatal error initializing TAE. 

[SFRAWD -2] Fatal error reading TAE parameters. 

[SFRAWD -3] Parameter ... is not in the data set. 

[SFRAWD -4] The last time in ... is invalid. 

[SFRAWD -5] The first time in ... is invalid. 
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4.42 PROGRAM SFSTNS 


SFSTNS modifies the station information in a surface file. 


INPUT PARAMETERS 

SFFILE *Surf ace data file 

STNFIL ’Station information file 

ADDSTN ’Add station flag 

IDNTYP * STNM or STID 


PROGRAM DESCRIPTION 

SFSTNS updates the station information in a GEMPAK surface 
fiie. The station information includes the character station 
identifier, STID, the station number, STNM, the latitude, 

SLAT, the longitude, SLON, the elevation, SELV, the state, 

STAT, and the country, COUN. 

This information must be stored in a fixed format in the table 
file specified in STNFIL. See the default table for an example 
ot the required format. Station names may be included, but are 
not used. The current default GEMPAK surface station table 
tor US and Canadian stations is GEMTABL: SFSTNS .TBL 

ADDSTN is a logical parameter that indicates whether stations 
which are in the table file but not already in the surface file 
will also be added to the surface file. 

IDNTYP governs whether station numbers, STNM, or character 

identifiers, STID, will be used to identify stations in the 
table. 


EXAMPLES 


1. Update surface file called SURF.DAT with station 

information fromMYSTN.DAT, adding stations that are not 
already in the surface file. 

SFFILE - surf.dat 
STNFIL - mystn.dat 
ADDSTN - yes 
IDNTYP - st id 


ERROR MESSAGES 
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[ SFSTNS +1] WARNING! No stations were updated. 

[SFSTNS -l] Fatal error initializing TAE . 

[SFSTNS -2] Fatal error reading TAE parameters. 

[SFSTNS -3] Invalid input for IDNTYP; must be STID or STNM. 

[SFSTNS -4] STNFIL ... cannot be opened. 
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4.43 PROGRAM SNCFIL 


SNCFIL creates a new GEMPAK sounding file. 


INPUT PARAMETERS 

SNOUTF *Output sounding file 

SNPRMF *Sound i ng parameter packing file 

STNFIL *Station information file 

MRGDAT *Merge data file flag/part type 

TIMSTN *Time s / add i t i on a 1 stations 


PROGRAM DESCRIPTION 

SNCFIL creates a GEMPAK sounding file. The file may be 
created to store merged or unmerged data. Merged data files 
may contain any meteorological parameters with all parameters 
included at each level. Unmerged files store the mandatory 
data and s i gn i f i c an t - 1 e ve 1 temperature and wind data as 
separate parts. 

MRGDAT must be set to YES or NO. If YES, a merged file will 
be created. The parameters to be included in the file must be 
specified in SNPRMF. If MRGDAT is NO, an unmerged file 
will be created. In this case, the type of unmerged file 
file should be included following a slash. Type 1 will 
create a file with only mandatory data below 100 mb. Type 2 
specifies a file with mandatory and s i gn i f i c an t - 1 e ve 1 data 
below 100 mb, and type 3, which is the default, specifies 
a file with mandatory and s i gn i f i c an t - 1 eve 1 data below and 
above 100 mb . 

SNPRMF contains information about the parameters to be included 
in a merged file. Either a list of parameters or the name of a 
packing file may be entered. If a list is entered, the parameters 
must be s epar a t ed wi th semicolons; packing inf ormat i on may also 
be included after a slash with the minimum and maximum values 
and the resolution separated using dashes. For example, to 

tem P erature and dewpoint in a file without packing, 

SNPRMF maybeenteredas: 5 

SNPRMF - TMPC;DWPC 
To pack the data, use: 

SNPRMF - TMPC/-127-63- . 1 ; DWPC/ - 63 - 63 - 1 
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SNPRMF may also contain the name of a packing file. A packing 
file for data which is not to be packed contains a list of 
parameters with one parameter per line. In a file for packed 
data, each line must include the parameter name, the minimum 
and maximum data values and the resolution, all separated with 
spaces. The default packing file for merged sounding data is 
GEMTABL : SNPACK . TBL . 

If STNFIL is not blank, information about all the stations in 
STNFIL will be added to the data set. Space will be left in the 
file for the additional number of stations specified as the 
second parameter in TIMSTN. Note that an error will result if 
STNFIL is blank and TIMSTN does not request additional stations. 


EXAMPLES 

1. Create a sounding file called SOUND.DAT which will contain 
merged data. Use the standard station and parameter files 
Leave r o om in the file for 100 stations in addition to the 
stations in STNFIL and allow 15 times to be added to the 


file. 



SNOUTF 

a* 

s ound . d a t 

SNPRMF 

« 

g emt abl : snpa ck . tbl 

STNFIL 

* 

g emt abl : snstns . tbl 

MRGDAT 

= 

yes 

TIMSTN 


15/100 

Cr e a t e 

an 

unmerged sounding file called WEEKLY . SND which. 

will store mandatory and significant-level data below and 

above 

100 

mb. Do not add any stations now. Leave room for 

14 t ime s 

and 300 stations. 

SNOUTF 

5= 

we ekly . snd 

SNPRMF 

- 


STNFIL 

— 


MRGDAT 

— 

no/3 

TIMSTN 

as 

14/300 


ERROR MESSAGES 


[SNCFIL 

+ 2 ] 

t SNCF I L 

-i] 

[SNCFIL 

-2] 

[SNCFIL 

-3] 

[ SNCF I L 

-4] 

[ SNCF I L 

-5] 

[ SNCF I L 

-6] 

[SNCFIL 

-7] 


Negative number of stations entered -- 0 used. 
Fatal error initializing TAE . 

Fatal error reading TAE par ame t e r s . 

Error opening station file ... 

File does not include ro om for any stations. 
File does not include room for any times. 
SNPRMF is incorrectly specified. 

The file name for the new file is blank. 
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4.44 PROGRAM SNCNVRT 


SNCNVRT converts GEMPAK3 sounding files to GEMPAK 5 sounding files. 


INPUT PARAMETERS 

OLDSND *GEMPAK3 sounding file 

SNOUTF *Output sounding file 

TIMSTN *Time s / add i t i ona 1 stations 


PROGRAM DESCRIPTION 

SNCNVRT adds the data from a GEMPAK3 sounding file to a 
GEMPAK5 file. All the stations in the old file will be 
included in the output file. No AREA subset is available. 

If the GEMPAK file is an existing file, data will be added 
to the file. If the file does not exist, a new file will be 
created. When a new file is created, the maximum number of 
times and stations that may be added to the file can be 
specified in TIMSTN. If the number of times and/or stations 
requested is smaller than the number of t ime s / s t a t i ons in 
the GEMPAK3 file, the numbers from the original file will 
be used. 

There is no difference between GEMPAK4 and GEMPAK5 sounding 
files. This program is used to convert GEMPAK 3 files since 
later GEMPAK programs cannot be run with these files. 


EXAMPLES 

1. Create a new file, GEMPAK. SND which can contain a maximum 
of 4 times and 175 stations. Convert the data in the 
GEMPAK3 file, 84DEC12Z.SND , which contains 1 time and 
137 stations and write the converted data to GEMPAK. SND . 

OLDSND — 84de c 1 2z . snd 
SNOUTF - gempak.snd 
TIMSTN - 4/175 

2. Add the data in the GEMPAK3 file called 84DEC0Z . SND to 
the file GEMPAK. SND whi ch was created in example 1. 

Note that the values in TIMSTN will not be used, since 
the file already exists. 

OLDSND — 84de cOz . snd 
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SNOUTF = gempak.snd 
TIMSTN - 0 


ERROR MESSAGES 


[ SNCNVRT 

+ 1] 

[ SNCNVRT 

-1] 

[ SNCNVRT 

-2] 

[ SNCNVRT 

-3] 

[ SNCNVRT 

-4] 

[ SNCNVRT 

-5] 

[ SNCNVRT 

-6] 

[ SNCNVRT 

-7] 

[ SNCNVRT 

-8] 

[ SNCNVRT 

-9] 

[ SNCNVRT 

-10] 

[ SNCNVRT 

-11] 


Station ... cannot be added to file. 

Fatal error initializing TAE. 

Fatal error reading TAE parameters. 

Error opening GEMPAK3 file ... . 

Invalid vertical coordinate in data set. 
GEMPAK5 output file ... cannot be opened. 
Invalid parameters in existing output file. 
Unmerged output file cannot be used. 

Invalid number of times to create file. 
Invalid number of stations to create file. 
Time ... cannot be added to file. 

Time ... is already in output file. 
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4.45 PROGRAM SNCROSS 


SNCROSS draws cross sections through sounding data. 
INPUT PARAMETERS 


CXSTNS 

SNPARM 

SNFILE 

DATTIM 

VCOORD 

PTYPE 

YAXIS 

TAXIS 

LINE 

BORDER 

CINT 

WIND 

TITLE 

PANEL 

DEVICE 

CLEAR 

FILTER 

CURVE 

TEXT 


‘Cross-section station line 
“Sounding parameter list 
“Sounding data file 
“Da t e / t ime 

“Vertical coordinate type 

“Plot type/h:w ratio/margins 

“Ystrt/ystop/yinc/lbl ;gln; tck 

*Timel-time2-tinc; lbl ;gln; tck 

“Line c o 1 o r /d a s h/wi d t h/ 1 abe 1 

“Background c o 1 o r /da sh/wi d th 

“Contour i n t e r va 1 /mi n/max 

“Wind symbo 1 / s i z /wd th/ typ/hd s z 

“Title c o 1 o r / 1 i ne / t i 1 1 e 

“Panel 1 o c a t i on/ c o 1 o r /da sh/wi d th 

“Graphics device 

“Clear screen flag 

“Filter data flag 

“Curve fit type 

“Text size/font/hw flag 


PROGRAM DESCRIPTION 

SNCROSS draws cross sections using upper air sounding data. 
Th e cross-section line, CXSTNS , mu st be specified as a list 
of stations separated with s emi colons. Alternatively, a 
single station may be entered in CXSTNS and a list of’times 
in DATTIM. In this case, a time section will be drawn. 

Any parameter that can be computed from the data set 
parameters can be displayed. The parameter to be displayed 
is specified in SNPARM. If the value of SNPARM is I SEN, 
isentropes will be drawn. Note that SNPARM may also be 
THTA , i n wh ich case, potential t emp erature will be gridded 
and contoured. I SEN wi 1 1 fit splines to the station data 
in plot space and then check for tangled lines and untangle 
them if necessary. When SNPARM - ISEN, the results should 
be the same as in the G EMP AKA version of SNCROSS. 

The color, line dash and width of the contour lines are set 
in LINE. The contour interval or spacing of the isentropes 
is given in CINT. The curve type to be used for generating 


4-90 



GEMPAK PROGRAMS 


the i sentropes is set in CURVE. The current curve types are 
piecewise straight lines, cubic splines with linear ends, 
cubic splines with parabolic ends, and cubic splines wit 
extrapolated ends. It is recommended that the curve type be 
set to 2 to use a general spline. 

The plot background consists of a pressure axis, a horizontal 
axis with the station names, a hatched region indicating the 
part of the plot below the Earth’s surface, and vertical lines 
at each station. The station lines are solid up to the level 
a t wh ichwind data stop, and are dashed fr om there to the 
level at which temperature data stop. The color for the entire 
background is given in BORDER. 

The horizontal axis represents a straight line between the 
first and last stations. The positions of intervening stations 
are proportional to the perpendicular projections of the actual 
positions onto the section line. All of these calculations 
are done in lat/lon coordinates. If the plot is a time section, 
the times will be displayed on the x axis with the earliest 
time at the left. If the first character in TAXIS is an R, the 
earliest time will appear on the right. The rest of the TAXIS 
i npu t is ignored. 

The vertical coordinate may be specified as LIN, LOG, or STUVE; 
SKEWT is not valid in this program. The bottom and top limits 
for the y axis are specified in YAXIS, but the axis labelling 
specifications are ignored. 


EXAMPLES 

1 . Dr aw isentropes in color 5, line type 2, labelling 
every line. Plot the background in color 1. Use a 
logarithmic scale for the y axis and a 4 Kelvin 
interval for the isentropes. Drawwind barbs in color 
number 2 at the stations. 


CXSTNS 

SNPARM 

SNFILE 

DATTIM 

VCOORD 

PTYPE 

YAXIS 

TAXIS 

LINE 

BORDER 

CINT 

WIND 

TITLE 


sep ; ggg ; jan; ckl ;ahn;chs 


i s e n 

s ound . s nd 
last 
pres 
log 


5/2/1 

1 

4 

bm2 
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PANEL = 0 

DEVICE - vt 
CLEAR = yes 
FILTER - yes 
CURVE - 2 

TEXT 

2. Draw a time section using all the data in the file for 
IAD. 

CXSTNS = iad 
DATTIM = all 


ERROR MESSAGES 

[ SNCROSS -1] 
[ SNCROSS -2] 
[SNCROSS -3] 
[SNCROSS -4] 
[SNCROSS -5] 
[SNCROSS -6] 
[SNCROSS -7] 
[SNCROSS -8] 
[SNCROSS -9] 
[SNCROSS -10] 
[SNCROSS -11] 
[SNCROSS -12] 
[SNCROSS -13] 
[SNCROSS -14] 
[SNCROSS -IS] 
[SNCROSS -16] 
[SNCROSS -17] 
[SNCROSS -18] 


Fatal error initializing TAE . 

Fatal error reading TAE parameters. 

GEMPLT initialization error. 

Vertical coordinate for isentropes must be PRES 
There are no times in the file. 

Fewer than four stations/ times were selected. 
Data buffer is too small. 

Temperature or pressure data not available. 

The station ... cannot be found in the data set 
Error setting up graph; check invalid LOG axis. 
Input ... for PTYPE is invalid. 

Input ... for YAXIS is invalid. 

Parameter ... is not computable. 

Parameter ... is a character. 

The grid coordinates cannot be defined. 

The station has no data at time .... 

Data at time ... is not in the file. 

Wind data cannot be computed. 
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4.46 PROGRAM SNDELT 


SNDELT deletes data from a sounding data file. 


INPUT PARAMETERS 

SNFILE * Sound i ng data file 

DATTIM *Da t e / t ime 

AREA *Data area 


PROGRAM DESCRIPTION 

SNDELT deletes data from a sounding data file. The stations 
to be deleted are specified in AREA; the times to be deleted 
are given in DATTIM. 

If AREA is set to DSET or ALL, all the data present at the 
times in DATTIM will be deleted, along with the headers for 
those t ime s . 

If AREA specifies a list of stations to be deleted, only 
the data for those stations at the specified times will be 
deleted. In this case, the station and time headers in 
the file remain unchanged. 


EXAMPLES 

1. Delete all the data at time 850722/1200. 

SNFILE - 85 j u 1 . snd 

DATTIM - 850722/1200 

AREA = d s e t 

2. Delete data at IAD for the most recent time in the file. 

SNFILE - realtime. snd 
DATTIM - last 
AREA « @ i a d 


ERROR MESSAGES 


t SNDELT 

-1] 

Fatal 

error initial 

i z i n g TAE . 

[ SNDELT 

-2] 

Fatal 

error reading 

TAE par ame t e r s 

[ SNDELT 

-3] 

Error 

deleting data 

at t ime .... 
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4.4 7 PROGRAM SNEDIT 


SNEDIT adds data in a sequential edit file to a sounding file. 


INPUT PARAMETERS 

SNEFIL *Sounding edit file 

SNFILE *Sounding data file 

TIMSTN *Time s / add i t i ona 1 stations 


PROGRAM DESCRIPTION 

SNEDIT adds information from a sequential file to a 
GEMPAK sounding file. It can be used to add or replace 
data at a station. The program can now be used to add 
data to either merged or unmerged data sets. 

The data to be added must be in a text file, SNEFIL. This 
file can be created using F as the output device in SNLIST. 
MRGDAT , which can be set in SNLIST, is a flag indicating 
whether the data to be written will be unmerged or written 
as separate parts. If the edit file contains unmerged 
data, the part name must be included and the parameters 
must be in the order expected by GEMPAK. 

If SNFILE exists, the data will be added to the file. If it 
does not already exist, a new file will be created. A new 
file will be a merged or unmerged data set, depending on the 
type of data in SNEFIL. The maximum number of stations and 
times allocated in a new file will be read from TIMSTN. 

For merged edit data files, the parameters to be edited must 
be specified at the beginning of the edit file. For example: 

PARM=PRES ; TMPC ; EWPC ; DRCT ; SPED ; HGHT 

The parameter line must contain the string, PARM, and =. 

Note that SNL 1ST will write SNPARM «•= xxxx;yyyy;... wh i ch 
is valid. If the parameter list must be continued on the 
next line, the last character on the current line mu s t be a 
semicolon. If SNFILE is an existing file, the parameters 
listed in the edit file mu st be exactly those in the data set. 
For unmerged data, information about the parameters is not 
necessary. 

The information to be added to the sounding file consists of 
station information for a station, followed by station data. 
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followed by station information for the next station, etc. 
Station information must be listed as KEYWORD = value, where 
the valid keywords are STID, STNM, SLAT, SLON , SELV, and TIME. 
The TIME keyword must be found for each station. Either STID 
or STNM is also required. If the station is not already in the 
file, the station identifier, station n umb er, latitude, 
longitude, and elevation found will be added to the file. If 
the station is already in the file, the station location and 
elevation will not be changed. Note that program SNSTNS can 
be used to modify station header values. 


EXAMPLES 

1. Add the data in the edit file, SNLIST.FIL, to the file, 
SOUND.DAT, which does not exist. 

SNEFIL = snlist.fil 
SNOUTF = sound.dat 
TIMSTN =1/ 10 

The file SNLIST.FIL follows: 

SNPARM = PRES ; TEMP; DWPT; DRCT; SPED; HGHT 
STNPRM = SHCW ; L I FT ; K I NX 

STID = I AD STNM = 72403 TIME = 841227/1200 


SLAT = 

38.98 I 

SLON = 

-77.46 ; 

SELV - 

85.0 

SHCW - 

14.11 

LIFT = 

25.60 

KINX = 

11.50 

PRES 

TEMP 

DWPT 

DRCT 

SPED 

HGHT 

1024 . 00 

2 . 40 

-1.40 

0.00 

0.00 

85 . 00 

1000.00 

2 . 20 

-3.80 

-9999.00 

-9999.00 

277 . 00 

850.00 

5 . 60 

-9.40 

-9999.00 

-9999.00 

1592 . 00 

700.00 

-1.90 

-2.30 

-9999.00 

-9999.00 

3167.00 

500.00 

-15.70 

-22.70 

-9999.00 

-9999.00 

5790 . 00 

400.00 

-28.90 

-33.40 

-9999.00 

-9999.00 

7430 . 00 

300 . 00 

-44 . 30 

-93.30 

-9999.00 

-9999.00 

9420 . 00 

250 . 00 

-55.30 

-104.30 

-9999.00 

-9999.00 

10620 . 00 

200.00 

-66.50 

-115.50 

-9999.00 

-9999 . 00 

12000 . 00 

150.00 

-65.90 

-114.90 

-9999.00 

-9999 . 00 

13750 . 00 

100 . 00 

-66.90 

-115.90 

-9999.00 

-9999.00 

16220 . 00 

The file 

SOUND . DAT 

will be 

created. 

1 1 wi 1 1 i 

contain the 


parameters PRES ; TEMP ; EWPT ; DRCT ; SPED ; HGHT . Space for a 
maximum of 1 time and 10 stations will be allocated. The 
time 841227/1200 and station IAD will be added to the 
file and then the data for that time and station will be 
added. Note that the stability indices will be ignored. 
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2 . Ad d the data in the edit file, I AD .FIL, to the file, 
REALT IME . SND , which already exists. 

SNEFIL - realtime. snd 
SNOUTF - iad.fil 
TIMSTN - 0 

The file IAD.FIL follows: 


STID - I AD STNM - 72403 TIME - 900823/0000 


SLAT - 

38 . 

98 ; 

SLON - -77 

.46 SELV 

oo 

0 

TTAA 

0 






PRES 


TMPC 

DWPC 

DRCT 

SPED 

HGHT 

1008.00 


19.20 

18.90 

0.00 

3.00 

-9999.00 

1000.00 


19.00 

18.20 

-9999.00 - 

9999.00 

150.00 

850.00 


14.40 

13.30 

160.00 

8.00 

1538.00 

700.00 


5.40 

4 . 30 

240.00 

6.00 

3157.00 

500.00 


-8.10 

-9 . 80 

215.00 

7.00 

5850.00 

400.00 


-19.90 

-49 . 90 

225.00 

11.00 

7540.00 

300.00 


-34.10 

-64.10 

235.00 

18.00 

9620.00 

250.00 


-44.30 

-9999.00 

245.00 

16.00 

10870.00 

200.00 


-54.30 

-9999.00 

250.00 

22.00 

12330 . 00 

150.00 


-60.50 

-9999.00 

265.00 

14.00 

14140.00 

100.00 


-64.50 

-9999.00 

235.00 

2.00 

16620.00 

TTBB 

0 






PRES 


TMPC 

DWPC 




1008.00 


19 . 20 

18.90 




850.00 


14.40 

13.30 




684.00 


3.60 

1 . 20 




602.00 


0.00 

-1.30 




481 .00 


-9 . 90 

-12.50 




PPBB 

0 






HGHT 


DRCT 

SPED 




0.00 


0.00 

3.00 




914.00 


120.00 

12.00 




1219.00 


130.00 

10.00 




1829.00 


180.00 

7.00 




2134.00 


195.00 

6.00 




2438.00 


210.00 

6.00 





ERROR MESSAGES 


[SNEDIT 

+ 1] 

Data for . 

[SNEDIT 

-1] 

Fatal erro 

[SNEDIT 

-2] 

Fatal erro 

[SNEDIT 

-3] 

The edit f 


.. has been added to the file, 
r initializing TAE. 
r reading the TAE parameters, 
ile cannot be opened. 
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[SNEDIT -4] 
[SNEDIT -5] 
[SNEDIT -6] 
[SNEDIT -7] 
[SNEDIT -8] 
[SNEDIT -9] 
[SNEDIT -10] 
[SNEDIT -11] 
[SNEDIT -12] 
[SNEDIT -13] 
[SNEDIT -14] 
[SNEDIT -15] 
[SNEDIT -16] 
[SNEDIT -17] 
[SNEDIT -18] 
[SNEDIT -19] 
[SNEDIT -20] 


List of parameters not found in edit file. 

Too many parameters in edit file. 

Edit file ... cannot be opened. 

Error creating new sounding file .... 

Edit file parms don’t match those in data set. 
First parameter not a valid vertical coordinate 
Error opening existing sounding file .... 

Cannot write merged data to unmerged data set. 
Sounding file can’t be created with MAXTIM <*= 0 
Sounding file can’t be created with MAXSTN <- 0 
The time ... was not added to the file. 

Th e station ... wa s not added to the file. 

Error writing data to file. 

The part ... has invalid parameters. 

Edit file ... has an invalid format. 

Station ... has an invalid format. 

Cannot write unmerged data to merged data set. 
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4.4 8 PROGRAM SNHODO 


SNHODO draws a hodograph of upper air data. 
INPUT PARAMETERS 


SNFILE *Sounding data file 

AREA ‘Data area 

LINE “Line c o 1 o r /d a s h/wi d t h/ 1 a be 1 

MARKER “Marker c o 1 o r / t yp e / s i z e /wi d t h/hw 

BORDER “Background c o 1 o r / d a s h/wi d t h 

TITLE “Title color/line/title 

XAXI S *Xstrt/xstop/xinc/lbl;gln;tck 

YAXI S *Ystrt/ystop/yinc/lbl;gln;tck 

LEVELS “Vertical levels 

VCOORD “Vertical coordinate type 

DATTIM “Date/ time 

CLEAR “Clear screen flag 

DEVICE “Graphics device 

PANEL “Panel 1 oc a t i on/ c o 1 o r /da sh/wid th 

TEXT “Text size/font / hw flag 


PROGRAM DESCRIPTION 

SNHODO draws a hodograph which represents the vertical 
distribution of the horizontal wind at a station. The winds 
are plotted in meters/ sec. 

The line color, line type, and width are specified in LINE. 
The marker color, type, and size are specified in MARKER. 

Th e levels at wh ich the line will be labeled are specified 
in LEVELS; VCOORD specifies the vertical coordinate for 
LEVELS . 

The x and y axis limits are specified in XAXIS and YAXIS. If 
these values are blank, the axes will be scaled to the actual 
data. 


EXAMPLES 

1. Plot hodographs for stations in Texas at 1200 UTC of 

the last day in the data set. Label the hodograph every 
50 mb from 1000 to 50 mb. Draw a solid line in color 3 
with a width of 5. Plot markers in color 2; draw the 
background in color 1 and the title in color 4. 
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SNFILE 

= 

[ nws . snf i 1 < 

AREA 

« 

@ t X 

LINE 

= 

3//5 

MARKER 

* 

2 

BORDER 

- 

1 

TITLE 

= 

4 

XAXIS 



YAXIS 

= 


LEVELS 

* 

1000-50-50 

VCOORD 

= 

PRES 

DATTIM 

— 

/ 1 2 

CLEAR 

= 

ye s 

DEVICE 

- 

c t 

PANEL 

- 

0 

TEXT 

- 

1 


2. Now plot the Denver hodograph for 1800 UTC on May 25 with 
labels every 2000 meters. Draw x and y axes from -50 to 
50 me ters/sec. 

AREA = ©den 

LEVELS - 1000-24000-2000 

VCOORD - hght 

DATTIM - 0525/18 

XAXIS - -50/50/10 

YAXIS - -50/50/10 


ERROR MESSAGES 

[ SNHODO -1] Fatal error initializing TAE. 

[ SNHODO -2] Fatal error reading TAE parameters. 

[SNHODO -3] Fatal error initializing GEMPLT . 

[SNHODO -4] The input for XAXIS is invalid. 

[SNHODO -5] The input for YAXIS is invalid. 

[SNHODO -6] Winds cannot be computed. 

[SNHODO -7] The vertical coordinate ... cannot be computed. 

[ SNHODO -8] The range along the x or y axis is invalid. 

[SNHODO -9] No winds can be plotted. 

[SNHODO -10] No valid stations were found. 
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4.49 PROGRAM SNLIST 


SNLIST lists upper air data from a sounding file for specified 
vertical levels and stations. 


INPUT PARAMETERS 


SNFILE 

AREA 

DATTIM 

SNPARM 

STNDEX 

LEVELS 

VCOORD 

OUTPUT 

MRGDAT 


♦Sounding data file 
♦Da t a area 
♦Da t e / t ime 

♦Sounding parameter list 
♦Stability indices 
♦Ve rtical levels 
♦Vertical coordinate type 
♦Output device (TPF) 

♦Merge data file flag/part type 


PROGRAM DESCRIPTION 

SNLIST lists parameters derived from an upper air data set for 
the requested stations and times. 

Parameters which can be computed at various levels in the 
data set should be specified in SNPARM. Parameters which have 
a single value at the station, such as stability indices, 
should be specified in STNDEX. 

Any parameters which can be computed from the data set 
parameters may be listed. If the vertical levels requested 
are not present in the data set, the data will be interpolated 
between existing levels. Data will also be interpolated to a 
new vertical coordinate system, if requested. 

Data from an unmerged data set can now be listed without 
being merged if MRGDAT is set to NO. If unmerged data 
are listed, SNPARM, VCOORD and LEVELS wi 1 1 be ignored. 


EXAMPLES 

1. List at the terminal the parameters in the data set for the 
station I AD at the mandatory levels. Also list the 
Showa Iter, lifted and K indices. 

SNFILE - g emd a t a : 84d e c . snd 
AREA = @iad 
DATTIM = 27/12 
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SNPARM 

STNDEX 

LEVELS 

VCOORD 

OUTPUT 


d s e t 

s h ow ;lift;kinx 

man 

pres 

t 


2. List the pressure, mixing ratio, u and v wind components 
and the Montgomery stream function on the isentropic 
surfaces from 300 to 400 degrees in 10-degree increments. 
Change the station to SLO. 


AREA 

SNPARM 

STNDEX 

LEVELS 

VCOORD 


@s 1 o 

p r e s ;mi xr ; uwnd ; vwnd ; psym 

300-400-10 
th t a 


3. List the unmerged data from I AD. 


AREA - @ i a d 

MRGDAT = no 


ERROR MESSAGES 


[SNLIST 

-1] 

[SNLIST 

-2] 

[SNLIST 

-3] 

[SNLIST 

-4] 

[SNLIST 

-5] 

[SNLIST 

-6] 

[SNLIST 

-7] 

[SNLIST 

-8] 

[SNLIST 

-9] 

[SNLIST 

-10] 

[SNLIST 

-11] 

[SNLIST 

-12] 


Fatal error initializing TAE . 

Fatal error reading TAE parameters. 

There are no parameters to be listed. 

Level parameter ... cannot be computed. 

Parm ... is character; cannot be listed. 
Vertical coordinate ... cannot be computed. 
Station parameter ... cannot be computed. 
Parm ... is character; cannot be listed. 

No stations reported at these times. 

Too many parameters; ... cannot be added. 
The data set is merged. 

No inputs for levels. 
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4.50 PROGRAM SNMAP 


SNMAP plots sounding data on a map. 


INPUT PARAMETERS 


AREA 

GAREA 

SNPARM 

DATTIM 

LEVELS 

VCOORD 

SNFILE 

COLORS 

WIND 

MAP 

LATLON 

MARKER 

TITLE 

CLEAR 

PANEL 

DEVICE 

PROJ 

FILTER 

TEXT 


♦Da t a area 
*Gr aph i c s area 
♦Sounding parameter list 
♦Da t e / t ime 
*Ve rtical levels 
♦Vertical coordinate type 
♦Sounding data file 
♦Co lor list 

♦Wind symbo 1 / s i z/wd th/ 1 yp/hds z 

♦Map c o 1 o r /da sh/wi d th 

♦Line c o 1 o r /d a sh/wi d t h/ 1 a be 1 / i nc 

♦Marker c o 1 o r / type / s i z e /wi d th/hw 

♦Title col or/ 1 ine/ t i 1 1 e 

♦Cl ear screen flag 

♦Panel 1 o c a t i on/ c o 1 o r /d a sh/wi d th 

♦Graphics device 

♦Map projection/angles/margins 

♦Filter data flag 

♦Text size/font/ h w flag 


PROGRAM DESCRIPTION 

SNMAP plots sounding data parameters at station locations 
on a map. Any level or station parameter that can be 
c omp uted can be displayed. 

The order of the input parameters determines the output 
location of the data. Parameters are separated by 
semicolons. A position may be skipped by entering two 
consecutive semicolons or entering parameter SPAC or BLNK. 
The following chart shows the placement of the data around 
the station. The number indicates the position of the 
parameter in the SNPARM list: 

7 

1 9 3 

2 + 4 

5 10 6 

8 
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Station data will be filtered; i.e., overlapping stations 
will not be plotted, if FILTER is set to YES. If a 
parameter is BLNK , the filter will not allocate any space 
for that parameter. The parameter SPAC may be used to reserve 
the space with the FILTER option so that later calls will plot 
the same stations after filtering, provided that the same 
number of parameters is specified. Data will be centered on 
the station location for the last parameter listed if an @ 
is appended to the name. 

Using WIND, either wind barbs or wind arrows can be plotted, 
the base being at the station location. Any valid marker type 
may be plotted at the stations using MARKER. 

Conditions can be specified for the parameters. The 
conditions are documented in the SNPARM documentation. Note 
that individual parameters can be scaled using these 
conditional functions. For example, TMPC*10will plot 
temperature multiplied by 10. 


EXAMPLES 

X. Plot a polar stereographic map of stations in the Maryland 
area with temperature, dewpoint, and height in positions 
1, 2, and 3. Plot the data at 850 mb for the latest time 
in the file, MARCH. SND: 


AREA 

- 

md 

GAREA 

— 

md 

SNPARM 

= 

tmp c ; dwp c 

DATTIM 


last 

LEVELS 

- 

850 

VCOORD 

* 

pres 

SNFILE 


ma rch. snd 

COLORS 

« 

2 ; 3 ; 4 

WIND 


bk2 

MAP 

* 

2 

LATLON 



MARKER 

- 

3 

TITLE 


5 

CLEAR 

= 

yes 

PANEL 

* 

0 

DEVICE 


v t 

PROJ 

= 

np s 

FILTER 

= 

no 

TEXT 




2. Plot a Mercator map of the USA with the 500-mb temperatures 
to the left of the station and height to the right. Plot 
the lifted index center bel o w the station in position 8. 
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The data plotted is from February 12, 1985 at 1200 GMT. 
Th e VT 100 terminal screen will not be cleared prior to 
plotting the map. Turn the wind barbs off. 

AREA 
GAREA 
SNPARM - 
DATTIM - 
SNFILE - 
WIND 
CLEAR 
PROJ 


u s 
us 

; tmpc ; ; hgh t ; ; ; ; 1 if t 
850225/1200 
25febl2z . snd 
0 

no 
me r 


ERROR MESSAGES 

[ SNMAP 

+ 1] 

[SNMAP 

-1] 

[ SNMAP 

-2] 

[SNMAP 

-3] 

[SNMAP 

-4] 

[SNMAP 

-5] 

[SNMAP 

-6] 


Parameter ... cannot be computed. 

Fatal error initializing TAE. 

Fatal error reading TAE parameters. 
Fatal error initializing GEMPLT . 

No valid levels have been input. 

A range of levels is invalid in SNMAP . 
Winds cannot be computed. 


4-104 



GEMPAK PROGRAMS 


4.51 PROGRAM SNMOD 


SNMOD moves selected sounding data from an input sounding file to 
an output sounding file. 


INPUT PARAMETERS 


SNFILE 

SNOUTF 

SNPARM 

AREA 

DATTIM 

LEVELS 

VCOORD 

TIMSTN 

MRGDAT 

IDNTYP 


♦Sounding data file 

♦Output sounding file 

♦Sounding parameter list 

♦Data area 

♦Da t e / t ime 

♦Ve rtical levels 

♦Vertical coordinate type 

♦Times/additional stations 

♦Merge data file flag/part type 

* STNM or STID 


PROGRAM DESCRIPTION 

SNMOD takes data from a GEMPAK sounding file, SNFILE, and 
writes it to an output sounding file, SNOUTF. This program 
can be used to subset the original data set by time and/or 
stations and to change the levels and vertical coordinate. 

The output data set may be either a merged or unmerged data 
set. If the output file does not exist, a new file will be 
created by this progr am . I f a n e w file is created, MRGDAT 
will determine whether the output file will be a merged or 
unmerged data set. Note that an error will result if the 
input data set is merged and the output data set is unmerged. 

The value of SNPARM wi 1 1 be ignored if the output data set is 
unmerged. If a merged data file is to be created, SNPARM 
will specify the parameters to be included. In this case, 
if SNPARM - DSET, the parameters in the input data set will be 
used. If the output file is an existing merged data set, the 
par ame ters specified in SNPARM mu st be in the data set. In 
this case, if SNPARM = DSET, the parameters in the output file 
wi 11 be used. 


If the output file is to be created, the maximum number of 
times and stations which can be stored in the file must be 
specified in TIMSTN. 


IDNTYP has been added to specify whether station numbers or 
station identifiers will be used in referring to stations. 
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Generally, IDNTYP should be set to STNM to use station 
numbers. However, in data sets which do not have station 
numbers, or which have a single invalid number, character 
station identifiers must be used and IDNTYP must be set 
to STID. 


EXAMPLES 


Put the mandatory data for stations in an area centered 
on IL at time 841227/1200 into a new unmerged sounding 
file, NEW.SND, which may contain 2 times and 150 
stations. 


SNFILE 

SNOUTF 

SNPARM 

AREA 

DATTIM 

LEVELS 

VCOORD 

TIMSTN 

MRGDAT 

IDNTYP 


gemdata : 84dec . snd 
new. snd 
d s e t 
i 1 

27/12 

man 

pres 

2/150 

no 

s t nm 


2 . 


Add the stations at time 841227/00 to the file created in 
exampl e 1 . 


DATTIM = 27/00 


3. Create an isentropic data set containing the specified 

parameters with levels every 5 degrees. Include all the 
US stations. 


SNFILE 

= 

gemdata : 84dec . snd 

SNOUTF 

SC 

tht a . snd 

SNPARM 

AREA 


thta;pres;mixr; uwnd ; vwnd ; psym 
@us 

DATTIM 


all 

LEVELS 

= 

250-450-5 

VCOORD 


th t a 

TIMSTN 

= 

4/150 


ERROR MESSAGES 

[ SNMOD 

-1] 

[ SNMOD 

-2] 

[SNMOD 

-3] 

[ SNMOD 

-4] 

[SNMOD 

-5] 


Fatal error initializing TAE . 

Fatal error reading TAE parameters. 

Output file cannot be unmerged type. 

A new file cannot be created with no times. 

A new file cannot be created with no stations. 
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[ SNMOD - 6 ] 
[ SNMOD - 7 ] 
[ SNMOD - 8 ] 
[ SNMOD - 9 ] 
[SNMOD -10] 
[SNMOD -11] 
[SNMOD -12] 
[SNMOD -13] 


Parameter ... cannot be computed. 

... is character type; cannot be added to file 
Output vertical coordinate ... cannot be used. 
No valid parameters were specified. 

Time ... cannot be added to the output file. 
Station ... cannot be added to output file. 
Parms do not match those in output data set. 
Output file name is blank. 
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4.52 PROGRAM SNPROF 


SNPROF draws profiles of upper air data. 
INPUT PARAMETERS 


SNFILE 

DATTIM 

AREA 

SNPARM 

LINE 

PTYPE 

VCOORD 

STNDEX 

STNCOL 

WIND 

WINPOS 

MARKER 

BORDER 

TITLE 

DEVICE 

YAXIS 

XAXIS 

FILTER 

CLEAR 

PANEL 

TEXT 

THTALN 

THTELN 

MIXRLN 


♦Sound ing data file 
♦Da t e / t ime 
♦Da t a area 

♦Sounding parameter list 
♦Line c o 1 o r /da sh/wi d th/ 1 a be 1 
♦Plot type/h:w ratio/margins 
♦Vertical coordinate type 
♦Stability indices 
♦Stability index color 
♦Wind symbo 1 / s i z /wd th/ typ/hds z 
♦Wi nd position 

♦Marker c o 1 o r / 1 yp e / s i z e /wi d t h/hw 

♦Background c o 1 o r /d a sh/wi d t h 

♦Title c o 1 o r / 1 i ne / t i t 1 e 

♦Graphics device 

*Ys tr t/ys top/yinc/lbl ;gln; tck 

♦Xstrt/xstop/xinc/lbl ;gln;tck 

♦Filter data flag 

♦Cl ear screen flag 

♦Panel 1 o c a t i on/ c o 1 o r /da s h/wi d th 

♦Text size/font/hw flag 

♦ THTA color/ dash/width/ mn / mx / i n c 

♦THTE color/da sh/wi d th/mn/mx/ i nc 

♦MIXR c o 1 o r / d a s h/wi d t h/mn/mx/ i nc 


PROGRAM DESCRIPTION 

SNPROF draws profiles of upper air data. 

Any two parameters can be plotted; e.g., temperature and 
dewpoint temperature. These parameters are specified in SNPARM. 

he color, line type, and width are specified in LINE. Markers 
can be plotted by setting a color in MARKER. 

Th * ° f y aXiS iS s P ecified in PTYPE. LIN, LOG, STUVE , 

and SKEW are all valid. Note that SKEW is only valid when 
plotting temperature vs. pressure. The bottom and top of 

vayic X 1 S ; r al ° ng With an 1 nc r ement for labels, are set in 
YAAib. If the panel is changed to be less than the full 
screen, the default margins will often be too large. The 
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margins can be set explicitly in PTYPE. 

Winds are plotted in the right margin and can be specified in 
WIND. There are three positions available for plotting winds 
They are numbered 1 to 3, with 1 being the leftmost position. 
The location for the first wind to be plotted is specified 
inWINPOS. The winds in subsequent plots will be plotted 
in the next position, modulo 3, provided the screen is not 
cleared. These positions also apply to the station 
identifier, date/time and stability indices which are 
plotted above the diagram. 

Dry adiabats, moist adiabats and mixing ratio background 
lines can be add to the diagram. The attributes of these 
lines are specified in THTALN , THTELN , and MIXRLN. The 
format for the specification is: 

color / line type / width / start / stop / increment 

If the color is set to 0, no lines will be plotted. If the 
start, stop and increment are not set, defaults will be 
supplied. Note that these lines can only be drawn on plots 
of pressure versus t emp e r a t u r e . 

FILTER in this program determines whether the wind barbs are 
filtered. 


EXAMPLES 

1. Plot a Skew T chart for stations ACY , IAD, and HAT for the 
latest time in the file. Display the Showalter index and 
the bulk Richardson number. Plot TMPC in color 1 using 
a solid line with width of 5. Plot EWPC in color 2 using 
a dashed line with width of 5. Include theta, theta-e, 
and mixing ratio lines. Plot filtered wind barbs in color 
4 in position 3. 


SNFILE = 
DATTIM = 
AREA 

SNPARM - 
LINE 
PTYPE 
VCOORD = 
STNDEX = 
STNCOL = 
WIND 

WINPOS = 
MARKER = 
BORDER - 


s ound . snd 
last 

@acy ; iad;hat 
tmp c ; dwp c 
1 ; 2/ 1 ; 3/5 
skewt 
pres 

show; br ch 
1 

bm4 

3 

1 
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TITLE 

- 

4 

DEVICE 

* 

qm 

YAXIS 


XAXIS 



FILTER 

- 

ye s 

CLEAR 

- 

ye s 

PANEL 

« 

0 

TEXT 



THTALN 

- 

3/4 

THTELN 

- 

4/4 

MIXRLN 

- 

5/4 


Plot 
Plot 
theta 
no 


relative humidity with a logarithmic pressure axis, 
markers at the sounding data points up to 10 mb. The 
theta-e, and mixing ratio lines are turned off, and 
stability indices are computed. 


PTYPE 

- 

1 og 

STNDEX 

- 

SNPARM 

* 

r e lh 

YAXIS 

- 

1050/10 

THTALN 

- 

0 

THTELN 


0 

MIXRLN 

= 

0 


3. 


Plot equivalent potential temperature on the x axis 
potential temperature on the y axis. 


vs . 


PTYPE - lin 
SNPARM - thte 
VCOORD - thta 
YAXIS - 250/450/50 


ERROR MESSAGES 

[ SNPROF +5] 
[ SNPROF +4 ] 
[SNPROF +3] 
[SNPROF +2] 
[SNPROF +1] 
[ SNPROF - 1 ] 
[ SNPROF - 2 ] 
[SNPROF -3] 
[SNPROF -4] 
[SNPROF -5] 
[SNPROF -6] 
[ SNPROF - 7 ] 

[ SNPROF - 8 ] 

[ SNPROF - 9 ] 
[SNPROF -10] 


Stability indices are specified with color - 0 
Parameter ... was requested with color set to 0. 
Winds cannot be computed. 

Parameter ... is a character type. 

Parameter ... cannot be computed. 

Fatal error initializing TAE. 

Fatal error reading TAE parameters. 

Fatal error initializing GEMPLT. 

The plot type ... is invalid. 

The the x-axis range must be specified in XAXIS. 
The two parms use different temperature units, 
lhe vertical coordinate ... cannot be computed. 
The range along the x or y axis is invalid. 

No valid stations were found. 

SNPARM has not been specified. 
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4.53 PROGRAM SNRAOB 


SNRAOB creates a GEMPAK sounding file from a raobs file created 
by the 604-line ingest program, TRAFIC. 


INPUT PARAMETERS 

RAOBFL *Raw raob data file 

SNOUTF *Output sounding file 


PROGRAM DESCRIPTION 

This program creates a GEMPAK sounding data file from data 
collected from the NWS 604 line and archived at GSFC. 

The program opens the raob file and displays the available 
times and the number of reporting stations. The user must 
enter the number corresponding to the desired time or exi 
the progr am . 

If the sounding file exists, the data will be added to the 
file. Otherwise, a new file is created. If the file exists 
and contains merged upper air data, an error message will be 
written. Mandatory and s i gn i f i c an t - 1 e ve 1 temperature and 
wind data below 100 mb will be stored in the sounding file 
as separate parts. 


EXAMPLE 


1. Create a sounding file TEST . SND from data in the file 
RAOBDAT, which is the logical name for the current 
r e a 1 t ime file. 


RAOBFL - 
SNFILE = 


ERROR MESSAGES 


[ SNRAOB 

+ 1] 

[ SNRAOB 

-1] 

[ SNRAOB 

-2] 

[ SNRAOB 

-3] 

[ SNRAOB 

-4] 

[ SNRAOB 

-5] 

[ SNRAOB 

-6] 


raobda t 
test .snd 


The sounding file ... has been created. 
Fatal error initializing TAE . 

Fatal error reading TAE parameters. 

Error opening raobs file. 

Error opening sounding file. 

Sounding file must contain unmerged data. 
Error reading raobs file. 
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[ SNRAOB - 7 ] 
[ SNRAOB - 8 ] 
[ SNRAOB - 9 ] 


Invalid time input. 

No more stations can be added to the file. 
No more times can be added to the file. 
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4.54 PROGRAM SNRAWD 


SNRAWD decodes upper air data from a file of raw data and adds 
it to a sounding data set. 


INPUT PARAMETERS 

SNEFIL * Sound i ng edit file 
SNFILE * Sound i ng data file 
DATTIM *Da t e / t ime 

ADDSTN *Add station flag 


PROGRAM DESCRIPTION 

SNRAWD decodes raw upper air data and adds it to a sounding 
data set. 

The raw data should be included in a text file. The sounding 
data set must already exist and must be an unmerged data file. 
The individual reports are distinguished by searching for PP 
or TT. 

DATTIM is used to specify the range of data to be included in 
the file. Both the earliest and latest times to be included 
must be specified with the complete date/ time. If the earliest 
time is not included, data from a 24-hour period will be used. 
Since the year and month are not included in the report data, 
the year and month from the last time will be used. Therefore, 
the data may not span months, except that data from the last 
day of the previous month will be decoded if the last time is 
the first day of the month. 


EXAMPLES 

1. Decode the data archived on April 26 in the file APR26.UAIR 
and store it in the file APRIL. SND. 

SNEFIL = apr26.uair 
SNFILE - april.snd 

DATTIM *= 890426/0000-890427/0300 


ERROR MESSAGES 

[SNRAWD +2] Time ... is not within specified range. 

[SNRAWD -1] Fatal error initializing TAE. 
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[ SNRAWD 

-2] 

[ SNRAWD 

-3] 

[SNRAWD 

-4] 

[SNRAWD 

-5] 

[SNRAWD 

-6] 


Fatal error reading TAE parameters. 

An input bulletin is too long. 

Time ... is an invalid last time. 

Time ... is an invalid first time. 

The sounding file is not an unmerged data set. 
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4.55 PROGRAM SNSTNS 


SNSTNS modifies the station information in an upper air file 


INPUT PARAMETERS 

SNFILE * Sound i ng data file 

STNFIL ^Station information file 

ADDSTN *Add station flag 

IDNTYP *STNM or STID 


PROGRAM DESCRIPTION 

SNSTNS updates the station information in a GEMPAK upper air 
file. The station information generally consists of the 
character station identifier, STID, the station number, STNM, 
the latitude, SLAT, the longitude, SLON , the elevation, SELV, 
the stat, STAT and the country, COUN . 

This information, along with the station name, which is not 
used, mu st be stored in a fixed fo rma t in the table file 
specified in STNFIL. The current GEMPAK upper air station 
table for US and Canadian station is GEMTABL : SNSTNS . TBL . 

An upper air file containing stations for the world is 
GEMTABL : SNW3RLD . TBL . 

ADDSTN is a logical parameter which indicates whether stations 
which are in the table file but not already in the upper air 
file will also be added to the upper air file. 

Either STID or STNM may be used to key on the desired station 
identifier by setting the desired value in IDNTYP. 


EXAMPLES 

1. Update upper air file called SOUNDINGS.DAT with station 
information fromMYSTN.DAT adding stations that are not 
already in the file. 


SNFILE 

STNFIL 

ADDSTN 

IDNTYP 


s ound ing s . da t 
my s t n . d a t 
yes 
s t i d 


ERROR MESSAGES 


4-115 


GEMPAK PROGRAMS 


[ SNSTNS +1] 
[ SNSTNS - 1 ] 
[SNSTNS -2] 
[SNSTNS -3] 
[SNSTNS -4] 


WARNING! No stations were updated. 

Fatal error initializing TAE . 

Fatal error reading TAE parameters. 

Invalid input for IDNTYP ; must be STID or STNM. 
STNFIL ... cannot be opened. 
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4.56 PROGRAM SNTSER 


SNTSER draws a time series at a sounding station. 
INPUT PARAMETERS 

SNFILE * Sound i ng data file 

DATTIM *Da t e / t ime 

TAXI S ‘Time 1-time 2-tinc;lbl;gln;tck 

LEVELS ‘Vertical levels 

VCOORD * Vertical coordinate type 

SNPARM * Sound i ng parameter list 

STNDEX ‘Stability indices 

AREA *Data area 

PTYPE *Plot type/h:w ratio/margins 

YAXIS *Ystrt/ystop/yinc/lbl;gln;tck 

BORDER ‘Background c o 1 o r / d a s h/wi d t h 

LINE ‘Line c o 1 o r / d a s h/wi d t h/ 1 a be 1 

MARKER ‘Marker c o 1 o r / t yp e / s i z e /wi d th/hw 
TITLE ‘Title c o 1 o r / 1 i n e / t i t 1 e 

CLEAR ‘Clear screen flag 

PANEL ‘Panel 1 o c a t i on/ c o 1 o r / d a s h/wi d t h 

DEVICE ‘Graphics device 


PROGRAM DESCRIPTION 

SNTSER draws a time series plot for a sounding station. 

Only one parameter may be plotted at a time. This parameter 
may be a level parameter defined in SNPARM or a stability 
index defined in STNDEX. If both SNPARM and STNDEX have 
values, SNPARM wi 1 1 be used. If a level parameter is 
defined in SNPARM, a vertical level must also be set in 
LEVELS . 

The type of y axis can be set in PTYPE. Either LIN or LOG 
is valid. The limits on the y axis can be set in YAXIS. 

The times to plot are specified in DATTIM. Only those 
times specified will be plotted. Lines will be drawn 
connecting the data points, provided that no more than two 
points are missing between segments. 

The time axis is specified in TAXIS using the usual GEMPAK 
date/time conventions, including FIRST and LAST. If 
TAXIS has no increment, a reasonable value is chosen. 

If TAXIS is blank, the limits are taken from the first and 
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last valid times set in DATTIM. 

When the user desires multiple lines on the same graph, TAXIS 
and YAXIS should be explicitly set. Then LINE, MARKER, and 
TITLE can be varied for successive combinations of LEVELS 
VCOORD, AREA, and SNPARM. 

The he i gh t - 1 o -wi d t h ratio of the plot may be specified in 
PTYPE following a /. If no value is entered, a value of 
0.5 will be used. 


EXAMPLES 


Plot the temperature in Celsius at the 500-mb level over 
BWI for the times between 26/12 and 31/00 in the file. 

Use the default line, marker, title and background colors 
Label the Y axis from -40 degrees to 0 degrees in 
increments of 10 degrees. Label the T axis from 26/00 
to 31/12 with internally chosen increments. 


2 . 


SNFILE 

- 

gemda t a : r eal t ime . snd 

DATTIM 

ss 

26/12-31/00 

TAXIS 

- 

26/00-31/12 

LEVELS 


500 

VCOORD 


pres 

SNPARM 

* 

tmp c 

STNDEX 



AREA 

- 

@bwi 

PTYPE 

= 

/ . 5 

YAXIS 

*= 

-40/0/10 

BORDER 

vm 

1 

LINE 

= 

2/1/1 

MARKER 

= 

1 

TITLE 

= 

5 

CLEAR 

=* 

ye s 

PANEL 


0 

DEVICE 

* 

vt 

Now, us 

ing 

the same file, plot a time 

1 i f ted 

index at OKC for all times in 

line with a 

thick, short dash pattern 

seal i ng 

for 

both axes. 

DATTIM 

_ 

all 

TAXIS 

= 


SNPARM 



STNDEX 


lift 

AREA 

= 

©okc 

YAXIS 

ss 



series of 
he file. 


the 

Plot 


the 


4-118 



GEMPAK PROGRAMS 


LINE 


ERROR MESSAGES 

[ SNTSER +1] 
[ SNTSER -1] 
[SNTSER -2] 
[ SNTSER - 3 ] 
[SNTSER -4] 
[SNTSER -S] 
[SNTSER -6] 
[SNTSER -7] 
[ SNTSER - 8 ] 
[SNTSER -9] 
[SNTSER -10] 


2/2/3 


WARNING, more than one parameter specified. 
Fatal error initializing TAE . 

Fatal error reading TAE parameters. 

Fatal error initializing GEMPLT . 

Input for LEVEL is invalid. 

Error defining graph coordinates. 

No points found for plot. 

Session not interactive. 

No parameters entered. 

Pa rm . . . not calculable. 

P a rm ... is a character. 
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APPENDIX A 
GEMPAK PARAMETERS 


This appendix contains a list of the GEMPAK parameters.. 
Algorithms used in computing these parameters are also included. 
The following constants are used in the computations: 


KAPPA - Poisson’s constant 

G = Gravitational constant 
GAMUSD » Standard atmospheric lapse rate 
RDGAS - Gas constant for dry air 
PI - Circumference / diameter 


2/7 

9 . 80616 m/sec/sec 
6.5 K/km 
287.04 J/K/kg 
3.14159265 


References for some of the algorithms: 

Bolton, D. , 1980: The computation of equivalent potential 

t emp e r a t u r e . , Mo n t h 1 y We ather Revi ew , 108, pp 1046-1053. 

Wallace, J .M. , P.V. Hobbs, 1977: Atmospheric Science, Academic 
Press, 467 pp . 

Miller, R.C., 1972: Notes on Severe Storm Forecasting 

Procedures of the Air Force Global Weather Central, 

AWS Tech. Report 200. 


TEMPERATURE PARAMETERS 

TMPC - Temperature in Celsius 
TEMP - T emp erature in Celsius 
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TMPF - Temperature in Fahrenheit 

TMPK - Temperature in Kelvin 

THTA - Potential temperature in Kelvin 

STHA - Surface potential temperature in Kelvin 

STHE - Surface equivalent potential temperature in Kelvin 

TVRK - Virtual temperature in Kelvin 

THTV - Virtual potential temperature in Kelvin 

TVRC - Virtual temperature in Celsius 

TVRF - Virtual temperature in Fahrenheit 

TMAX - Maximum 24 hour temperature in Celsius 

TMIN - Minimum 24 hour temperature in Celsius 

SSTC - Sea surface temperature in Celsius 

TMPC - C TMPF - 32 ) * 5 / 9 

TMPC - TMPK - 273.16 

THTA - TMPK * ( 1000 / PRES ) ** KAPPA 

STHA - TMPK * ( 1000 / PALT ) ** KAPPA 

TVRK - TMPK * ( 1 + ( .001 * MIXR ) / .62197) ) / 

( 1 + ( .001 * MIXR ) ) 

THTV = TVRK * ( 1000 / PRES ) ** KAPPA 

MOISTURE PARAMETERS 

EWPC - Dewpoint in Celsius 
EWPT - Dewpoint in Celsius 
EWPF - Dewpoint in Fahrenheit 
EWPK - Dewpoint in Kelvin 
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DPDC 

DPDF 

DPDK 

MIXR 

MIXS 

SMXR 

SMXS 

RELH 

THTE 

VAPR 

VAPS 

LHVP 


D e wp oint depression in Celsius 
D e wp oint depression in Fahrenheit 
Dewpoint depression in Kelvin 
Mixing ratio in g/kg 
Saturated mixing ratio in g/kg 
Surface mixing ratio in g/kg 
Surface saturated mixing ratio in g/kg 
Relative humidity in percent 

Equivalent potential temperature in Kelvin 
Vapor pressure in millibars 
Saturation vapor pressure in millibars 
Latent heat of vaporization in J/kg 


DPDx - TMPx - DWPx 

MIXR = .62197 * C E / ( PRES - E ) ) * 1000 

E - VAPR * [ 1.001 + ( PRES - 100) / 
900 * .0034 ] 

VAPR = 6.112 * EXP ( ( 17.67 * DVPC ) / 

( DWPC + 243.5 ) ) 

Bolton, MVR (1980) 

MIXS - .62197 * ( Es / ( PRES - Es ) ) * 1000 

Es - VAPS * [ 1.001 + ( PRES - 100 ) / 

900 * .0034 ] 

VAPS - 6.112 * EXP ( ( 17.67 * TMPC ) / 

( TMPC + 243.5 ) ) 

Bolton, MWR (1980) 

SMXR - .62197 * ( E / ( PALT - E ) ) * 1000 

E = VAPR * [ 1.001 + ( PALT - 100 ) / 

900 * .0034 ] 

VAPR = 6.112 * EXP ( ( 17.67 * DWPC ) / 

( DWPC + 243.5 ) ) 

SMXS - .62197 * ( Es / ( PALT - Es ) ) * 1000 

Es - VAPS * [ 1.001 + ( PALT - 100 ) / 

900 * .0034 ] 
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HGHT 

HGTM 

HGTK 

HGTD 

HGFT 

HGFH 

HGFK 

HGML 

DHGT 

MHGT 

ZMSL, 


VAPS = 6.112 * EXP ( ( 17.67 * TMPC ) / 
( TMPC + 243.5 ) ) 

RELH - VAPR / VAPS * 100 


THTE - THTM * EXP [ ( 3.376 / TLCL - .00254 ) * 

C MIXR * ( 1 + .81 * .001 * MIXR ) ) ] 

THTM = potential temperature for moist air 
= TMPK * ( 1000 / PRES ) ** E 
E = 2 . / 7 . * ( 1 - ( .28 * .001 * MIXR ) ) 
TLCL - temperature at the LCL in Kelvin 
Bolton, NfiVR (1980) 


VAPR 


VAPS 


LHVP 


6.112 * EXP [ ( 17.67 
Bolton, MWR (1980) 


DWPC ) / 

( DWPC + 243.5 ) ] 


6.112 * EXP [ ( 17.67 
Bolton, MWR (1980) 


TMPC ) / 

( TMPC + 243.5 ) ] 


( 2.501 - .00237 * TMPC ) * 10E6 

Bolton, XRVR ( 1980) 


HEIGHT PARAMETERS 


- Height in meters 

- Height in me t e r s 

- Height in kilometers 

- Height in dekameters 

- Height in feet (3.28084 * HGHT) 

- Height in hundreds of feet 

- Height in thousands of feet 

- Height in miles (6.2137E-04 * HGHT) 

- Dry hydrostatic height in meters 

Moist hydrostatic height in meters 

Z000, Z900, Z850, Z800 - Estimated height at a pressure 

level 


A- 4 


GEMPAK PARAMETERS 


DHGT and MHGT are computed using the hypsometric 
equation and integrating from the surface pressure 
to the desired level. DHGT (MHGT) is computed 
without (with) the influence of moisture. 

DHGT = HB + ( RDGAS / G ) * LN ( PBOT / PTOP ) * TAVE 

HB = height of the bottom of layer 

PBOT = pressure at bottom of layer 
PTOP = pressure at the top the layer 
TAVE = average t emp erature of the layer 

MHGT = HB + ( RDGAS / G ) * LN ( PBOT / PTOP ) * TVAVE 

HB - height of the bottom of layer 
PBOT = pressure at bottom of layer 
PTOP = pressure at the top the layer 
TVAVE = average virtual t emp erature of the layer 

Zxxx - [ To * ( 1 - ( PRES / ALTM ) ** 

( RDGAS * GAMUSD / G ) ] / GAMUSD 
Zxxx = height at the pressure level 
ZMSL : 1013.25 mb 
Z000 : 1000 mb 

Z900 : 900 mb 
Z850 : 850 mb 
Z800 : 800 mb 

To = sea level temperature in U.S. Std. Atmos. 
- 288 K 


PRESSURE AND ALTIMETER PARAMETERS 


PRES - Station pressure in millibars 

PRES is the actual pressure at a level as reported with 
upp e r air data. 

PMSL - Mean sea level pressure 

PMSL is reported with surface data. 

PALT - Surface pressure in millibars from ALTI 
ALTI - Altimeter setting in inches of mercury 
ALTI is reported with surface data. 

ALTM - Altimeter setting converted to millibars 
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SALT - Abbreviated standard altimeter setting 
AMSL - Abbreviated mean sea level pressure 

PMSL - PRES * EXP ( ( G * SELV ) / ( RDGAS * TVAVE ) ) 
SELV — station elevation 

TVAVE - average virtual t emp erature be twe e n 
station and sea level 

- TVRK + ( DELTV / 2 ) 

DELTV - GAMUSD * SELV / 1000 
Wallace & Hobbs (1977) 

PALT - ALTM * 

( 1 - ( ( SELV / 1000 ) * GAMUSD /To ) ) ** 

( G / ( GAMUSD * RDGAS ) * 1000 ) 

SELV = station elevation in meters 

To - sea level temperature in U.S. Std. Atmos. 

- 288 K 

Wallace & Hobbs (1977) 

ALTM - ALTI * ( 1013.25 / 29.921 ) 

SALT - MOD ( ALTM * 10, 1000 ) 


WINDS 


UWND 

- U-component 

of 

the 

wind 

i n 

me t e r s / s e c ond 

VWND 

- V- component 

of 

the 

wi nd 

i n 

me t e r s / s e c ond 

UKNT 

- U- component 

of 

the 

wi nd 

i n 

knots 

VKNT 

- V* component 

of 

the 

wi nd 

i n 

knots 

DRCT 

- Wi nd direct 

i on 

i n 

degree 

s 



SPED - Wind speed in me t e r s / s e c ond 


SKNT - Wind speed in knots 

PSPD - Packed direction and speed in me t e r s / s e c ond (ddfff) 
PKNT - Packed direction and speed in knots (ddfff) 

GUST - Wind gusts in knots 
GUMS - Wind gusts me t e r s / s e c ond 
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UWND = - SIN ( DRCT ) * SPED 

VWND = - COS ( DRCT ) * SPED 

DRCT = AT AN 2 ( -UWND, -VWND ) * 180 / PI 

SPED = SQRT ( UWND ** 2 + VWND ** 2 ) 

SPED = SKNT / 1 . 9438 

PS PD - JDRCT * 100 + J SPED 

JDRCT = NINT ( DRCT / 5 ) 

J SPED - NINT ( SPED ) 

PKNT = JDRCT * 100 + J SKNT 

JDRCT = NINT ( DRCT / 5 ) 

J SKNT = NINT ( SKNT ) 


LIFTED CONDENSATION LEVEL (LCL) 


TLCL - Temperature in Kelvin at the LCL 
PLCL - Pressure in millibars at the LCL 


TLCL - [1 / ( 1 / ( EWPK - 56 ) + 

LN ( TMPK / EWPK ) / 800 )] + 56 
Bolton, MvVR ( 1980) 

PLCL * PRES * ( TLCL / ( TMPC + 273.16 ) ) ** ( 1 / KAPPA ) 
Poisson’s equation 


STABILITY INDICES 

Note: Default depths given in the definitions below are preceded 

by an ! and may be changed using the in-line depth 
specification, dz in the definitions bel o w defaults to the 
layer thickness in the sounding. 

BRCH - Bulk Richardson number 

BRCH - CAPE / ( 0,5 * U**2 ) 

CAPE = Convective Available Potential Energy 
U - magnitude of shear ( u2 - ul , v2 - vl ) 

ul,vl = average u,v in the 1 owe s t ! 5 0 0 m 
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u2,v2 - average u,v in the lowest 16000 m 

BVFQ - Brunt-Vai sala frequency 

BVFQ - SQRT ( ( G / THTA ) * STAB ) 

BVPD - Brunt-Vai sala period 

BVPD - 2. * PI / BVFQ 

BVSQ - Brun t -Va i s a 1 a frequency squared 

BVSQ - BVFQ ** 2 

CAPE - Convective available potential energy 

CAPE - GRAVTY * SUMP ( DELZ * ( THP - TH ) / TH ) 

SUMP » sum over sounding layers for which 

( THP - TH ) is greater than zero 
THP - potential temperature of a parcel from the 
lowest 1500 m of the atmosphere, raised dry 
ad i aba t i cal ly to the lcl and moist 
adiabatically thereafter 
TH - potential temperature of the environment 

CINS - Convective inhibition 

CINS - - GRAVTY * SUMN ( DELZ * ( THP - TH ) / TH ) 

SUMN - s um over sounding layers for wh i c h 

( THP - TH ) is less than zero. 

THP - potential temperature of a parcel from the 
lowest 1500 m of the atmosphere, raised dry 
adiabatically to the lcl and moist 
adiabatically thereafter 
TH = potential temperature of the environment 

CTOT - Cross Totals index 

CTOT - TD850 - T500 

TD850 - Dewpoint in Celsius at 850 mb 
T500 - Temperature in Celsius at 500 mb 

EQLV - Equilibrium level 

EQLV « level at which parcel undergoing moist adiabatic 
ascent ceases to be buoyant. 

KINX - K index 

KINX - ( T850 - T500 ) + TD850 - ( T700 - TD700 ) 

T850 - Temperature in Celsius at 850 mb 
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T500 = Temperature 
TD850 = Dewpoint in 
T700 = Temperature 
TD700 - Dewpoint in 


in Celsius 

a t 

500 

mb 

Celsius at 

850 

mb 


in Celsius 

a t 

700 

mb 

Celsius at 

700 

mb 



LAPS - Temperature lapse rate in layer 

LAPS - d (TMPK) / dz - d (TMPC) / dz 


LFCV - Level of free convection 

LFCV = level above which a parcel ascending moist 
a d i a b a t i c a 1 1 y is buoyant. 


LIFT - Lifted index 

LIFT = T500 - Tparcel 

T500 = Temperature in Celsius at 500 mb 

Tparcel - 500 mb temperature in Celsius 

of a lifted parcel with the average 
pressure, temperature and dewpoint of 
the layer 1100 mb above the surface 

MLMR - Mean mixed layer MIXR 

Ml MR » Average MIXR in lowest 1500 m 

mLTH - Mean mixed layer TUT A 

MLTH - Average THTA in lowest 1500 m 

RICH - Richardson number in layer 

RICH - BFVQ ** 2 / SHRM ** 2 

SEPA - Isentropic layer pressure thicness 

SEPA - Pressure difference over isentropic layer 15 K deep 

SHCW - S h owa Iter index 

SHOW - T500 - Tparcel 

T500 = Temperature in Celsius at 500 mb 

Tparcel = Temperature in Celsius at 500 mb 
of parcel lifted fr om 850 mb 

SHRD - Wind shear direction in layer 

SHRD - direction of [ du/dz, dv/dz ] 

SHRM - Wind shear magnitude in layer 
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SHRM = magnitude of [ du/dz, dv/dz ] 

STAB - THTA lapse rate in layer 
STAB - d (THTA) / dz 
STAP - THTA change with pressure 
STAP - - d (THTA) / dp 
SWET - SWEAT index 

SWET = 12 * TD850 + 20 * TERM2 + 2 * SKT850 + 

SKT500 + SHEAR 

TD850 = Dewpoint in Celsius at 850 mb 

TERM2 = MAX ( TOTL - 49, 0 ) 

TOTL - Total totals index 

SKT850 = 850 mb wind speed in knots 

SKT500 = 500 mb wind speed in knots 

SHEAR - 125 * SIN ( DIR500 - DIR850 ) + .2 

DIR500 = 500 mb wind direction 

DIR850 = 850 mb wind direction 

If TD850 is negative, then TD850 is set to 0. 
SHEAR is set to 0 if any of the following 
conditions is met: 

wind direction at 850mb is < 130 or > 250 

wind direction at 500mb is < 210 or > 310 

DIR500 - DIR850 < 0 
SPD500 <= 15 or SPD850 <= 15 
Miller (1972) 

TOTL - Total Totals index 

TOTL = ( T850 - T500 ) + ( TD850 - T500 ) 


T850 = Temperature 

i n 

Celsius 

at 

850 

mb 

TD 8 50 = D e wp o i n t in 

Cel 

s i u s at 

850 

mb 


T500 = Temperature 

i n 

Celsius 

a t 

500 

mb 


VTOT - Vertical Totals index 
VTOT = T850 - T500 

T850 - Temperature in Celsius at 850 mb 
T500 = Temperature in Celsius at 500 mb 


CLOUD PARAMETERS 

Cloud coverage may be defined using a cloud code, short code, 


A- 10 



GEMPAK PARAMETERS 


fractional coverage or numeric value. The valid values of these 
parameters ordered from least to greatest cloud coverage, are. 


Cloud Coverage 

Cloud Code 
xCLD 

Short 

missing 


_ 

clear 

CLR 

C 

thin scattered 

-SCT 

-s 

scattered 

SCT 

s 

thin broken 

-BKN 

-B 

broken 

BKN 

B 

thin ove r c a s t 

-OVC 

-0 

overcast 

OVC 

0 

thin obscured 

-X 

-X 

obscured 

X 

X 


Code Fractional Numeric 

xCLO CLCx 

0.000 0 

0.000 1 

0.166 6 

0.333 2 

0.500 7 

0.667 3 

0.833 8 

1.000 4 

1.000 9 

1.000 5 


The following lists the GEMPAK parameter definitions with an 
e x amp le using the s amp 1 e A I KWAY S cloud report: 


22SCT 80-BKN 2500VC 

Note that the character x may be replaced by L, M, or H, 
indicating low, mid or high clouds. Also note that the 
character T indicates the value of the parameter at the level 
o f ma x i mum cloud coverage. 


xCLD - Character cloud coverage code 
Examples: LCLD - SCT 
MCLD = - BKN 
HCLD = OVC 

TCLD - xCLD at maximum cloud coverage 
Example: TCLD ** OVC 

xCLO - Fractional cloud coverage 
Examples: LCLO - 0.333 
MCLO = 0.500 
HCLO - 1 .000 

TCLO - xCLO at maximum cloud coverage 
Example: TCLO - 1.00 

CLCx - Numeric cloud coverage 
Examples: CLCL - 2 
CLCM - 7 
CLCH - 4 

CLCT - CLCx at maximum cloud coverage 
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Example: CLCT - 4 


The next two parameters combine cloud coverage values from the 
three cloud levels. 

CLDS - Combined cloud coverage short code from three levels 
Example: CLDS = S-BO 

CMBC - Combined cloud coverage numeric from three levels 
Example: CMBC = 274 


The next set of parameters are combined cloud height and 
c 1 oud coverage. 

CLHx - Cloud height in hundreds of feet 
Examples: CLHL - 22. 

CLHM - 80. 

CLHH - 250. 

CLDx - Combined cloud height and short code 
Examples: CLDL - 22 S 
CLEM - 80 -B 

CLDH - 2500 

CLDT - CLDx at maximum coverage level 
Example: CLDT - 2500 

CQMx - Numeric combined cloud height and coverage combined as 
CLHx * 10 + CLCx 
Ex amp 1 e s : CQML - 222. 

COMM = 807 . 

CQMH - 2504 

Note: In the case when the sky is partially obscured, the 
value of 10000 is added on to the lowest reporting level. 
For example if AIRWAYS report is: 

-X M5 BKN 19 BKN 
CQML would equal 10053 


COMT - CQMx at maximum coverage level 
Example: COMT - 2504. 


The next two parameters combine the cloud height and coverage 
allowing up to three reports which do not necessarily correspond 
to low, middle and high level clouds. These parameters allow a 
means of storing cloud reports where there may be more than 
one report at a single level. 
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CHCx - Numeric combined cloud height and coverage combined as 
CLHx * 10 + CLCx , where x is the cloud report group 
number from 1 to 3. 

Examples: COM1 = 222. 

COM2 = 807. 

COM3 = 2504 

Note: In the case when the sky is partially obscured, the 
value of 10000 is added on to the first reporting level. 
For ex amp 1 e if A I KWAY S report is: 

-X M5 BKN 19 BKN 
CHC1 would equal 10053 

CHDx - Combined cloud height and short code 
Examples: CHD1 — 22 S 
CHD2 - 80-B 
CHD3 - 2500 


The following set of cloud parameters is the numeric WMO codes 
which are reported by airways data. 

CFRC - Fraction of celestial dome covered by cloud fromWvIO 
Code 2700 

CLFR - Fraction of celestial dome covered by all low and mid 
level clouds fr om VAIO Code 2700 

CTYL - Low-level cloud genera fromWMO Code 0513 

CTYM - Mid-level cloud genera from WO Code 0513 

CTYH - High-level cloud genera from WO Code 0513 

CBAS - Cloud base height from WO Code 1000 


WEATHER CODES 

WCOD - Character weather code 
WNUM - Numeric weather code 

The weather code WEOD may also be accessed as WTHR. The 
weather number consists of 3 parts. A, B, C where 
WNUM = A * 80 * 80 + B * 80 + C. 

Each part corresponds to one of the values: 


0 


(no value) 


1 

R 

(mod rain) 

41 

2 

L 

(mod drizzle) 

42 

3 

S 

(mod snow) 

43 
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4 

A 

(mod hail) 

44 



5 

T 

( thund e r ) 

45 



6 

H 

(haze) 

46 



7 

K 

( smoke ) 

47 



8 

D 

(dust) 

48 



9 

F 

( f og ) 

49 

ZR- 

(it frz rain) 

10 



50 

ZR+ 

(hvy frz rain) 

1 1 



51 

KW- 

(it rain s hwr ) 

12 



52 

RW+ 

(hvy rain shwr) 

13 

R- 

(it rain) 

53 

ZL- 

(it freezing drizzle) 

14 

R+ 

(hvy rain) 

54 

ZL+ 

(hvy freezing drizzle) 

1 5 

ZR 

(mod frz rain) 

55 

S W- 

(it s n ow s hwr ) 

16 

RW 

(mod rain shwr) 

56 

SW+ 

(hvy snow shwr) 

17 

L- 

(it drizzle) 

57 

IP- 

(it ice pellets) 

18 

L+ 

(hvy drizzle) 

58 

IP+ 

( hvy ice pellets) 

19 

ZL 

(frz drizzle) 

59 

SG- 

(it s n ow grains) 

20 

S- 

(it snow) 

60 

SG+ 

(hvy snow grains) 

21 

s+ 

(hvy snow) 

61 

SP- 

(It s n ow pellets) 

22 

sw 

(mod snow shwr) 

62 

SP+ 

(hvy snow pellets) 

23 

IP 

(mod ice pellet) 

63 

IPW 

(mod ice pellet shwr) 

24 

SG 

(mod snow grain) 

64 

IC- 

(it ice crystals) 

25 

SP 

(mod snow pellet) 

65 

IC+ 

(hvy ice crystals) 

26 

A- 

(It hail) 

66 

TRW 

(mod thunder shwr) 

27 

A+ 

(hvy hail) 

67 



28 

T- 

(it t hund e r ) 

68 



29 

T+ 

(hvy thunder) 

69 



30 

IF 

(ice fog) 

70 



31 

GF 

(ground fog) 

71 



32 

BS 

(blowing snow) 

72 



33 

BD 

(blowing dust) 

73 



34 

BY 

(blowing spray) 

74 



35 

BN 

(blowing sand) 

75 

IPW- 

(it ice pellet s hwr ) 

3 6 

IC 

(mod ice crystals) 

76 

IPW+ 

(hvy ice pellet shwr) 

37 

IN 

(ice needles) 

77 

TRW- 

(.It rain thunder shwr) 

38 

AP 

( sma 11 hail) 

78 

TRW+ 

(hvy rain thunder shwr) 

39 

KH 

( smoke , haze) 

79 



40 






The 

following correspond to 

a single character code: 

-1 

TORNA 

(tornado) 

-3 

WYTER (water 80001 ") 

-2 

FUNNE 

( f unn e 1 cloud) 





WTHR - Character weather codes for W604 data 
W604 - Numeric weather codes used for 604 data 


0 - 22 

1 - T 23 

2 — R 24 

3 - S 25 

4 - F 26 


L+ 

44 - 

R+S 

SP 

45 - 

RS- 

SG 

46 - 

RS + 

BY 

47 - 

ZL- 

BN 

48 - 

ZL+ 


A- 14 



GEMPAK PARAMETERS 


WIMO 

VMMO 


5 = 

H 

27 = BD 

49 

= 

TSW 

6 - 

K 

28 = IF 

50 

= 

TSW- 

7 - 

L 

29 = RS 

51 

- 

TSW+ 

8 - 

R+ 

30 = TRW 

52 

s* 

TRW- 

9 - 

R- 

31 - TR- 

53 

= 

TRW+ 

10 - 

S + 

32 = TR+ 

54 

55 

TRWA 

11 - 

s- 

33 = TS- 

55 

- 

R-S- 

12 - 

RW 

34 = TS+ 

56 

= 

R-S + 

13 - 

sw 

35 - ZR- 

57 

- 

R+S- 

14 = 

TR 

36 = ZR+ 

58 

= 

R+S + 

15 - 

TS 

37 = SW- 

59 

= 

TRW+A 

16 - 

ZL 

38 - SW+ 

60 

* 

TRW- A 

17 - 

ZR 

39 - SG- 

61 

= 

TORNA 

18 - 

IP 

40 = SG+ 

62 

2= 

FUNNE 

19 - 

GF 

4 1 - RW- 

63 

= 

WATER 

20 = 

BS 

42 = RW+ 




21 - 

L- 

43 - R-S 




Character WMO 

we ather code 




Numeric VAIO weather code 




The 

transf o rma t i o n is: 




0 


34 = BD+ 

67 


ZR 

1 

* 

35 = BD+ 

68 

= 

R-S- 

2 


36 - BS 

69 


RS 

3 

* 

37 = BS+ 

70 

* 

S- 

4 

- K 

38 - BS 

71 

* 

S- 

5 

- H 

39 - BS+ 

72 

* 

S 

6 

- D 

40 = 

73 

* 

S 

7 

- BD 

41= F 

74 

= 

S + 

8 

- BD 

42 = F 

75 

= 

sh- 

9 

- BD 

43 = F 

76 


in 

10 

- F- 

44 = F 

77 


SG 

11 

- GF 

45 = F 

78 


S- 

12 

- GF 

46 = F 

79 

* 

IP 

13 


47 = F 

80 

* 

RW- 

14 


48 = F 

81 

= 

RW 

15 

OB 

49 = F 

82 

= 

RW+ 

16 

at 

50 = L- 

83 

= 

RW-SW- 

17 

- T 

51 = L- 

84 

= 

RWSW 

18 

- Q 

52 = L- 

85 

= 

S W- 

19 

= FUNNE 

53 = L 

86 

= 

SW 

20 

a se 

54 = L+ 

87 

= 

SP- 

21 

= 

55 = L+ 

88 


SP 

22 

* 

56 = ZL- 

89 

= 

A- 

23 

s 

57 = ZL 

90 

= 

A 

24 

B 

58 = R-L- 

91 


R- 

25 

s 

59 = RL 

92 

= 

R 

26 

» 

60 = R- 

93 

= 

RS 

27 


61 = R- 

94 

= 

R+S + 
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28 

- 


62 = R 

95 = TRW- 

29 

« 


63 = R 

96 = TRW- A 

30 

- 

BD 

6 4 - R+ 

97 - TRW+ 

31 


BD 

65 - R+ 

9 8 - TD 

32 

33 


BD 

BD+ 

66 = ZR- 

99 - TRW+A 


WASH - Asheville surface data weather code 
RWSH - Numeric Asheville surface data 


The data are stored using 7-digits as 1234.567* The 
column digits have the following translation: 


Co 1 umn 1 : 

Co 1 umn 2 : 

Column 3: 

Column 4: 

0 - None 

0 = None 

0 = None 

0 


None 

1 - R- 

1 = None 

1 - S- 

1 

=* 

SW- 

2 - R 

2 = None 

2 - S 

2 


SW 

3 = R+ 

3 = None 

3 = S + 

3 


SW+ 

4 = RW- 

4 - L- 

4 - SP- 

4 


None 

5 = RW 

5 = L 

5 = SP 

5 


None 

6 = RW+ 

6 = L+ 

6 = SP+ 

6 

d 

None 

7 - ZR- 

7 = ZL- 

7 = none 

7 


SG- 

8 = ZR 

8 = ZL 

8 - IC 

8 

= 

SG 

9 = ZR+ 

9 = ZL+ 


9 

= 

SG+ 

Co 1 umn 5 : 

Co 1 umn 6 : 

Co 1 umn 7 : 




0 = None 

0 = None 

0 = None 




1 = IP- 

1 = F 

1 = K 




2 = IP 

2 = IF 

2 = H 




3 = IP + 

3 = GF 

3 = KH 




4 - None 

4 - BD 

4 - D 




5 = A 

5 = BN 

5 = BS 




6 - None 


6 = T 




7 - None 


7 = T+ 




8 = AP 


8 = TOR 

9 = Q 





PWTH - Character past weather WvlO code 
PVWM - Numeric past weather WvlO code 

The past weather WvlO numeric codes are: 

0 * Cloud covering less than 1/2 sky 

1 * Cloud covering more than 1/2 during part of 

period and less than 1/2 sky during part 

2 - Cloud covering more than 1/2 sky 

3 — Sandstorm, dust storm or blowing snow 

4 * Fog, ice fog, thick haze or thick smoke 
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5 = Drizzle 

6 = Ra i n 

7 = Snow, rain and snow mixed or ice pellets 

8 = Showers 

9 - Thunderstorm with or without precipitation 
The translation is: 

0 - 5 = L 

1 » 6 = R 

2 = 7 = S 

3 = BD 8 = RW 

4 - F 9 = T 


STATION PARAMETERS 

STID - Character station identifier 
STNM - Station number 

STNM is the 5-digit WMO identifier for upper air data; 6 
digits for surface data, usually the "WMO identifier with a 
zero appended . 

SLAT - Station latitude in degrees 

SLON - Station longitude in degrees; West longitude is negative 

SELV - Station elevation in meters 

RANG - Range in kilometers (specialized use) 

AZIM - Azimuth in kilometers (specialized use) 

LATI - Latitude in degrees from r ang e / a z imu t h 

LONG - Longitude in degrees from range/azimuth 

LATI and LONG are calculated from the RANG and AZIM using 
equations developed for AO I P S / RADPAK . 

DELT - Delta time in seconds (specialized use) 


MISCELLANEOUS PARAMETERS 
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VS BY 

VSBK 

P06 I 

P24 I 

P06M 

P24M 

SNOW 

HAIL 

HLSZ 

WHGT 

WPER 

DDEN - 

PSYM - 


- Visibility in statute miles 

- Visibility in kilometers 

Precipitation over last 6 hours in hundreths of inches 
Precipitation over last 24 hours in hundreths of inches 
Precipitation over last 6 hours in mill ime t e r s 
Precipitation over last 24 hours in millimeters 

- Snow depth in inches 

■ Ha i 1 flag 

• Hail size in centimeters 

■ Wave height in meters 
Wa ve period in seconds 
Density of dry air in kg/(m**3) 

Montgomery stream function in m* *2/ ( 100* s * *2 ) 
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APPENDIX B 

GRID DIAGNOSTIC FUNCTIONS 


The following describes the computation of GEMPAK grid diagnostic 
funct i ons . 


Each grid in a grid file is identified by a parameter name, 
time, level, and vertical coordinate. A scalar grid is a 
single grid, while a vector grid is composed of two grids 
containing the u and v components. 

The parameter name is used to retrieve a grid from the file, with 
a few exceptions: Certain special parameters will be computed 

from other data in the grid file if the parameter name itself is 
not found in the grid file. These special scalar parameters are 


TMPK 

THTA* 

MIXR* 

TVRK 

DRCT 


TMPC 

STHA 

MIXS 

TVRC 

SPED 


TMPF 

THTE* 

SMXR 


XVAL 


DWPC 

STHE 

SMXS 


YVAL 


EWPF 

RELH* 





EWPK 






where * indicates names 

wh i c h also ma y 

be used 

as operators. 

Other 

special scalar 

par ame ter 

names denote 

cons tant 

value grids: 



DTR 

E 

GRAVTY 

KAPPA 

PI 

RTD 

nnn 


Conversion factor for degrees to radians = PI / 180 

Base of natural logarithms = 2.71828182 
Gravitational constant = 9.80616 (note spelling) 

Gas c ons t an t / spe c i f i c heat = 2/7 
3.14159265 

Conversion factor for radians to degrees «= 180 / PI 

Any number (i.e., 2, -10.2, ... ) 


Another class of special parameter names provides information at 
grid points depending on the navigation of the grid file: 


CORL Coriolis force 
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LATR 

LONR 

XVAL 

YVAL 


CORL - 2. * OMEGA * SIN ( LATR ) 

Latitude in radians 
Longitude in radians 

Value of the x coordinate in graph coordinates 
Value of the y coordinate in graph coordinates 


Finally, scalar grids may be identified by their location within 
the grid file. The grid number must be prefixed with the symbol 
#. Note that grids may be renumbered as grids are added to or 
deleted fr om t h e f i 1 e . 


Vector grids are two separate grids containing the u and v 
components. Special vector parameter names may be used to identify 
the following vectors: 

OBS Observed wind 

GEO* Geostrophic wind 

AGE* Ageos trophic wind 

I SAL* Isallobaric wind 

THRM* Thermal wind 


where * indicates names that also may be used as operators 
that all of these wind vectors will have 
meters per second. The observed wind 
VWND in the grid file if the 
as UGRD and VGRD if the 


Note 

u and v components in 
mu st be stored as UWND and 
c omp onents are north relative and 
components are grid relative. 


GIFVFT 1 !In 1 rvrnDn Ver i iCal coordinate ma Y b e specified in GDATTIM, 
GLEVEL and GVCORD. However, any of these values may be 

overridden by in line parameters appended to an operand in the 

form of timealevel%ivcord. In line parameters are only allowed 

tor operands, since they modify parameters for individual grids. 

Grid operators may be nested, allowing a complicated diagnostic 
function to be computed. One limitation is that layer and 
time range operators expect to work on operands read directly 
from the grid file or computed from special names. 

01 1 ° wing llst of diagnostic operators, scalar operands are 
named Si and vector operands are Vi. Lower case u and v refer to 
the grid relative components of a vector. All meteorological 
output grids are in MKS units, except as noted. Operators using 
un ctions are described in the GEMPAK PARAMETER APPENDIX " 

niii^t lar f nd V K Ct u r differentia l operators are valid in any map 
projection for which the map scale factors can be computed. At 

• thlS f ° r the s 1 e r e 0 g r a ph i c , cylindrical and conic 

projections available in GEMPAK . In the definitions below, only 
he cartesian form of the operators is shown. The general 

given 1111631 C00rdinate forms involving the scale factors are not 
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The operators which are designated for use in polar coordinates are 
specific to that coordinate system. 


SCALAR OUTPUT GRID 


Algebraic and trignometric functions 
(angles are expressed in radians) 


ABS 

Absolute value 

ACOS 

Arc cosine function 

AS IN 

Arc sine function 

ATAN 

Arc tangent function 

ATN2 

Arc tangent function 

COS 

Cosine function 

EXP 

Exponential to real 

EXP I 

Exponential to integer 

LN 

Natural logarithm 

LOG 

Base 10 logarithm 

SIN 

Sine f unc t i on 

SQRT 

Square root 

TAN 

Tangent function 

ADD 

Add i t i on 

AVG 

Average 

MUL 

Mu ltiplication 

QUO 

Division 

SUB 

Sub t r ac t i on 


ABS (S) 

ACOS (S) 

AS IN (S) 

ATAN ( S ) 

ATN2 (SI, S2 ) - ATAN ( SI / S2 ) 
rn <1 ( ^ 1 

EXP (SI , S2 ) - SI ** S2 

EXP (SI, S2 ) - SI ** NINT ( S2 ) 

LN (S) « LOG (S) 

LOG (S) = LOGIO (S) 

SIN (S) 

SQRT (S) 

TAN (S) 

ADD (SI, S2 ) - SI + S2 

AVG (SI, S2) = ( SI + S2 ) / 2 

MUL (SI, S2 ) - SI * S2 

QUO (SI, S2 ) - SI / S2 

SUB (SI, S2 ) « SI - S2 


ADV Advection 

ADV ( S, V ) - 


u * DDX (S) + v * DDY (S) ) 


AVOR Absolute vorticity 

AVOR ( V ) - VOR ( V ) + CORL 


DDEN Density of dry air 

DDEN ( PRES, TMPC ) = PR_DDEN ( PRES, TMPC ) 

DDR Partial derivative with respect to R 

DDR ( S ) is computed using centered finite differences, with 
backward or forward differences at the boundary. 
Polar coordinates are assumed, and ( R, THETA ) maps 
into ( X, Y ) . 


DDT 


Time derivative 

DDT ( S ) = ( S (timel) - S (time 2) ) 
where the time difference 


/ ( timel - t ime 2 
is in seconds. 


) 


DDX Partial derivative with respect to X 

DDX ( S ) is computed using centered finite differences, with 
backward or forward differences at the boundary. 
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DDY 

DEF 

DIRN 

DIRR 

DIV 

DOT 

DTH 

FCNT 

FRNT 


JCBN 

KNTS 

LAP 


Partial derivative with respect to Y 

DDX ( S ) is computed using centered finite differences, with 
backward or forward differences at the boundary. 

Total deformation 

DEF ( V ) - ( STR (V) ** 2 + SHR (V) ** 2 ) ** .5 

North relative direction of a vector 
DIRN ( V ) = PR W J)RCT ( UN (V), VN (V) ) 


Grid relative direction of a 
DIRR ( V ) - PR_DRCT ( u, v ) 


vector 


Dive rgence 

DIV C V ) - DDX ( u ) + DDY C v ) 
Vector dot product 

DOT ( VI. V2 ) - Ul * U2 + vl * V2 


rvruV^ c | erivative with respect to THETA 
TH C S ) is computed using centered finite differences, with 
backward or forward differences at the boundary, 
ro 1 a r ( c^o rd^i na t e s are assumed, and ( R, THETA ) 


maps 


FCNT^ l (\ ? *1 6 center of a Polar coordinate grid 

( S ) can be computed only for lat/lon grids which have 
been mapped to polar (R, THETA) coordinates and 
or which the center lat/lon have been stored 
wi th each grid. 


Fr ontogenesis 

FRNT ( THTA, V ) - CONV * 1/2 * MAG ( GRAD ( THTA ) ) * 

( DEF * COS (2 * BETA) - DIV ) 

CONV = unit conversion factor 
- 1.08E4 * 1.E5 

BETA = AS IN C -DDX (THTA) * COS (PS I) / 
MAG ( GRAD (THTA) ) ) 

PS I - 1/2 ATAN ( SHR / STR ) 


Jacobian determinant 

JCBN ( SI, S2 ) - DDX (SI) * DDY (S2) - DDY (Si) 


Convert meters / second to 
KNTS ( S ) - PR_MSKN (S) 

- S * 1 .9438 


kno t s 


* DDX ( S2 ) 


Laplacian operator 

LAP ( S ) - DIV ( GRAD (S) ) 
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LAV 

LDF 

MAG 

MASS 

MDIV 

MIXR 

MSDV 

NORM 

PLAT 

PLON 

POLF 

PVOR 


Layer average 
LAV ( S ) = ( 


(2 levels) 

S ( 1 eve 1 1 ) + S ( 1 eve 1 2 ) ) / 


2 . 


Layer difference (2 levels) 

LDF ( S ) - S (levell) - S (leve!2) 


Magnitude of a vector 

MAG ( V ) = PR^SPED ( u, v ) 


Mass per unit volume in a layer 

MASS - 100 * LDF (PRES) / ( GRAVTY * (levell - leve 2) 

The 100 converts mb to Pascals. Levell and eve 
are also converted to Pascals when VCOORD=PRES . The 
volume is expressed in units of m * m * (units o t e 
vertical coordinate). This is an operand. 


Laye r - ave r ag e mass 

MDIV ( V ) = DIV ( 


divergence 

[ MASS * LAV (u), MASS * LAV (v) 


) 


Mi xing ratio 
MIXR ( IXVPC, PRES 


PR_MIXR ( DWPC, PRES ) 

The units are g/g internally, but g/kg 
on output. 


Layer-average mass- scalar flux 
MSDV ( S, V ) - DIV ( t MASS * 

MASS * 


d i ve rgence 

LAV ( S ) * LAV ( u ) , 

LAV (S) * LAV (v) ] ) 


Note: MASS is computed using the in-line 
parameter values for V rather 


for S. 


than 


those 


scalar vector component normal to a crosssection 
NORM ( V ) = DOT ( V, unit normal vector ) 

If the starting point for the cross section is 
on the left, the unit normal vector points 
into the cross section plane. 


Latitude at each point in polar coordinates 

PLAT ( s ) Note: only the header, which contains the center 

latitude and longitude, is used. 

Longitude at each point in polar coordinates 

PLON ( S ) Note: only the header, which contains the center 

latitude and longitude, is used. 

Coriolis force at each point in polar coordinates 

POLF ( S ) Note: only the header, which contains the center 

latitude and longitude, is used. 


Potential 
PVOR ( V ) 


vor t i ci ty in a layer 
= - AVOR ( VLAV (V) ) * LDF ( THTA ) / 
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RELH 

RICH 

ROSS 

SAVG 

SAVS 

SDIV 

SHR 

SM5S 

STAB 

STR 

TANG 


( 100 * LDF ( PRES ) ) 

The 100 converts millibars to Pascals. 

Units are Kelvins / seconds / Pascals (note that 
GRAVTY is not included). PVOR wo rks on a layer 
in PRES or THTA coordinates. 


Relative humi d i t y 

RELH ( TEMP. DWPT ) = PR RELH ( TEMP. DWPT ) 


Richardson stability numb er in a layer 

RICH ( V ) = GRAVTY * DZ * LDF (THTA) / 

( LAV (THTA) * MAG ( VLDF (V) ) ** 2 ) 

Note that DZ = change in height across the layer. 
RICH can be evaluated in PRES, THTA or HGHT 
vertical coordinate. 

Rossby number 

ROSS ( VI, V2 ) = MAG ( INAD ( VI , V2 ) ) / 

( CORL * MAG ( VI ) ) 

Average over whole grid 

SAVG ( S ) = average of all non-missing grid point values 

Average over subset grid 

SAVS ( S ) = average of all non-missing grid point values in 
the subset area 


Flux divergence of a scalar 

SDIV ( S, V ) - S * DIV ( V ) + DOT ( V, GRAD ( S ) ) 


Shear deformation 

SHR ( V ) = DDX ( v ) + DDY ( u ) 


Smo oth scalar grid using a 5-point 
SM5S ( S ) = .5 * S (i.j) + .125 * 


smo o the r 
( S (i+l.j) 
S (i-l.j) 


S (i , j+1) + 
S (i.j-1) ) 


Thermodynamic s tabi 1 i ty wi thin a 
STAB ( TMPC ) = LDF ( TMPC ) / DZ 


layer (lapse rate) 


DZ - change in height across the layer 
= - ( RDGAS / GRAVTY ) * LAV (THTA) * 

( LAV (PRES) / 1000 ) ** KAPPA * 
LDF (PRES) / LAV (PRES) in THTA 
coordinates 

The units are degrees / kilometer. 


Stretching deformation 

STR ( V ) - DDX ( u ) - DDY ( v ) 


scalar vector component tangential to a croi 
TANG ( V ) . DOT ( V, uni, tangent vector ) 


s e c t i on 
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TAV 

TDF 

TIIT A 

THTE 

UN 

UR 

VN 

VOR 

VR 

XAV 


XSUM 


YAV 


If the starting point for the cross section is 
on the left, the unit tangent vector points to 

the right. 


Time average (2 times) 

TAV (S) = ( S (timel) + S (time 2) ) / 2. 

Time difference (2 times) 

TDF (S) = S (timel) - S (time 2) 


Potential temperature 

TUT A ( TMPC, PRES ) = PR_THTA ( 

Equivalent potential temperature 
THTE (PRES, TMPC, DWPC) = PR_THTE 


TMPC, PRES ) 

(PRES, TMPC, DWPC) 


North relative u component 
UN ( V ) = zonal wind component 


Grid relative u c omp o n e n t 
UR ( V ) - u 

North relative v component 

VN ( V ) = meridional wind component 


Vor t i c i ty 

VOR ( V ) - DDX ( v ) - DDY ( u ) 


Grid relative v component 
VR ( V ) = v 


Av e r a g e 
XAV (S) 


long a grid r o w 

( S (XI) + S (X2 ) + ... + s (KXD) ) / KNT 
KXD - number of points in row 
KNT - number of non-missing points in row 
XAV for a row is stored at every point in 
In polar coord. XAV is the average along a 


that row. 
radial . 


) 


Sum along a grid row . 

XSUM (S) =• ( S (XI) + S (X2 ) + ... + S (KXD) 

KXD = number of points in row 

XSUM for a row is stored at every point in that row. 
For polar coord., XSUM is the sum along a radial. 


A v e r a g e 
YAV (S) 


value along a grid col umn 

- ( S (Yl) + S (Y2) + . . . + S (KYD) ) / KNT 
KYD = number of points in column 
KNT - number of non-missing points in column 
YAV for a col umn is stored at every point in that 
col umn . For polar coordinates, YAV is the average 
around a circle. If the theta coordinate starts at 
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0 degrees and ends at 360 degrees, the first radial 
is not used in computing the average. 

YSUM Sum along a grid column 

YSUM (S) = ( S (Yl) + S (Y2) + ... + s (KYD) ) 

KYD = number of points in column 

YSUM for a column is stored at every point in that 
column. For polar coordinates, YSUM is the sum 
around a circle. If the theta coordinate starts at 
0 degrees and ends at 360 degrees, the first radial 
is not used in c omp uting the s um . 


VECTOR OUTPUT GRID 


AGE Ageos trophic wind 

AGE ( S ) = [ U COBS) - u (GEO( S ) ) , V (OBS) - v (GEO(S)) ] 

C I RC Circulation (for cross sections) 

CIRC ( V, S ) = [ TANG (V), s ] 


GEO 


Geos t rophi c wind 

GEO ( S ) - [ - DDY (S) 

CORL ] 


const / CORL, DDX (S) 


const 

S 

vert coord 

GRAVTY 

ZMSL 

none 

GRAVTY 

HGHT 

PRES 

1 

PSYM 

THTA 

100/RO 

PRES 

HGHT 

RO = PRJDDEN 

( PRES, TMPC 


GRAD Gradient of a scalar 

GRAD ( S ) = [ DDX C S ) , DDY ( S ) ] 

I NAD Inertial advective wind 

INAD ( VI , V2 ) - [ DOT ( VI, GRAD (u2) ), 

DOT ( VI , GRAD ( v2 ) ) ’ ] 

I SAL Isallobaric wind 

ISAL ( S ) - [ - DDT ( v (GEO(S)) ) / CORL 

DDT ( u (GEO( S ) ) ) / CORL] 

KNTV Convert meters / second to knots 

KNTV ( V ) = [ PR^MSKN (u), PRJtfSKN (v) ] 


LTRN 


!'TDM r r aVerag x d transport of a scalar 
LTRN ( S , V ) = [ MAS S * LAV ( S ) * LAV ( u ) , 

MASS * LAV (S) * LAV (v) ' ] 


const / 
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NORMV 

ROT 

SMUL 

SM5V 

TANGV 

THRM 

VADD 

VAVG 

VAVS 

VECN 

VECR 

VLAV 

VLDF 

VMUL 


Note: MASS is computed using the in-line parameter values 
for V rather than those for S. 

Vector component normal to a cross section. 

NORMV ( V ) = NORM ( V ) * unit normal vector 


Coordinate rotation 
ROT ( angle, V ) = 


u * COS (angle) + v * SIN (angle), 
u * SIN (angle) + v * COS (angle) J 


Multiply a scalar with each component of a vector 
SMUL (S,V)=[S*u, S*v] 


Smo oth vector grid using a 5-point 
SM5V ( V ) = .5 * V (i,j) + .125 * 


smoo the r 
( V ( i+1 , j ) 
V ( i - 1 , j ) 


V ( i , j+1) + 

V ( i , j - 1 ) ) 


Vector component tangential to a cross section. 
TANGV ( V ) - TANG ( V ) * unit tangent vector 


Thermal wind 

THRM ( S ) - [ u (GEO(S)) (levell) 

v (GEO(S)) (levell) 


u (GEO(S)) (level2) , 
v (GEO(S)) ( 1 eve 1 2 ) ] 


Add the components of two vectors 
VADD ( VI , V2 ) = [ u 1 +u2 , vl+v2 ] 

Av erage over wh olegrid 

VAVG ( V ) - average of all non-missing grid point values 


Av erage over subset grid 

VAVS ( V ) = average of all non-missing grid point values in 
the subset area 

Create a vector grid from two north relative scalar components 
VECN ( SI , S2 ) = [ SI , S2 ] 

Create a vector grid from two grid relative scalar components 
VECR ( SI , S2 ) - [ SI , S2 ] 


Layer average for a vector 

VLAV ( V ) = [ ( u (levell) + u (level2) ) / 2 . , 

( v (levell) + v (level2) ) / 2. J 

Layer difference for a vector 

VLDF ( V ) = [ u (levell) - u (levell), 

v (levell) - v ( 1 eve 1 2 ) ] 

Multiply the components of two vectors 
VMUL ( VI, V2 ) = [ ul *u2 , vl*v2 ] 
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VQUO 

VSUB 


Divide the c omp onents of 
VQUO ( VI, V2 ) = [ ul/u 

Subtract the components 
VSUB ( VI, V2 ) = [ ul-u 


two vectors 
2, vl/v2 ] 

of two vectors 
2 , vl - v2 ] 
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APPENDIX C 

LINE TYPES AND WEATHER SYMBOLS 


This appendix shows the line types and weather symbols which have 

been added to GEMPAK. 

There are ten distinct line type: 

1 - solid 

2 - sho rt dashed 

3 - medium dashed 

4 - long dash short dash 

5 - long dash 

6 - long dash three short dashes 
7-1 ong dash dot 

8 - long dash three dots 

9 - medium dash dot 

10-dotted 

These patterns can be expanded or compressed by prefixing the 
single digit with a number from 1 to 9. A prefix of 1 
compresses the pattern, 2 is the default and prefixes 3 -- 9 
expand the basic pattern. For example, 32 expands line type 2 
while 12 c omp resses the s ame pattern. 

These patterns are shown on the next page. On the following 
page, the we a t h e r s ymb ols are sh own . 
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APPENDIX D 

CHANGES FROM GEMPAK4 TO GEMPAK5 


This appendix describes some of the changes in GEMPAK5 from 
GEMPAK4 . Many changes have been made in the programs and input 
parameters, and several new programs have been added. 

TIME SERIES: 

The following time series programs have been added: 

SFGRAM - surface time series 
SNTSER - sounding time series 
GDTSER - grid time series 

HODOGRAPH : 

The new program, SNHODO, draws hodographs. 

DECODER PROGRAMS: 

SFRAWD and SNRAWD decode surface and upper air data from an 
archived text file. These programs may require modification 
to use with data archived in formats different than that used at 
Goddard Space Flight Center. 

GDCROSS : 

GDCROSS is a new program which draws a vertical cross section 
between two points in a grid field. Both scalar and vector 
quantities, including circulation, can be displayed. 

GENERALIZED GRID DIAGNOSTICS: 

The grid diagnostic functions have been generalized so that 
differential operators are valid for grids in any projection. 

GRID DIAGNOSTIC FUNCTIONS: 

Several new grid diagnostic functions have been added. New 


functions include: 





JCBN - 

Jacobian 

d e t e rmi nan t 


SM5S - 

5-point 

smoother 

for 

scalars 

SM5V - 

5-point 

smoother 

for 

vectors 
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Internally, all vectors are stored as grid relative c omp o n e n t s . 
Operators to specify grid relative or north relative vector 
components have also been defined. 

GRID PACKING: 

Grids stored in grid files can be packed in several ways. These 
include GRIB, the standard gridded binary format, and NMC, the 
standard format used for NMC history tapes. 

GRID PROFILE: 

GDPROF has been extensively modified to produce plots similar to 
SNPROF . Winds can be plotted in the right margin. Potential 
temperature, equivalent potential temperature, and mixing ratio 
background lines can be added to the plot. A listing of the 
quantity being profiled can also be obtained. 

GRID POINT LOCATION: 

Grid point locations for GPOINT in GDPROF and GDTSER and for 
CXSTNS in GDCROSS may be specified using surface or upper air 
character or integer station identifiers, as well as lat/lon 
coordinates or grid x and y coordinates. 

LINE DASHING PATTERNS: 

New line dashing patterns have been defined which include dotted 
patterns. The capability to scale these patterns has also been 
added . 

CONTOUR ALGORITHM: 

The contouring algorithm has been rewritten to eliminate spurious 
loops and overdrawn lines. In addition, the contours no longer 
overdraw the labels. 

WEATHER SYMBOLS : 

Weather symbols are plotted in SFGRAM and SFMAP when WSYM is 
spec i f i ed . 

LATITUDE /LONGITUDE LINES: 

La t i tud e / 1 ong i t ude grid lines can be drawn in all the programs that 
plot data on a map. s 

COLORS : 

A specific color may now be assigned to a color number by entering 
=NAME after the color number. For example, COLORS-1 = r ed ; 2-b 1 u e , will 
set color number 1 to red and color number 2 to blue. If the last 
character in COLORS is a ?, then the current color names will be 
listed. 

SYMBOL ATTRIBUTES: 

All symbols, including text, arrows, barbs, markers, and weather 
symbols, have an additional attribute which specifies the line 
width. 
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CHANGES FROM GEMPAK4 TO GEMPAK5 


STREAMLINES: 

The algorithm for drawing streamlines has been rewritten. The 
a r r ow heads on the stre am lines will no longer change size as the 
plot is zoomed. Also, GD STREAM correctly passes grid relative 
winds to the streamline package. 

GRAPHICS UTILITY PROGRAMS: 

The GP programs have been modified and reorganized for GEMPAK5 . A 
new program, GPMAP , can be used to define a coordinate system and 
add a map and 1 a t i t ud e / 1 ong i t ud e grid lines. 

MAP FILES: 

A US county map has been added. Maps containing only the Northern 
and Southern hemispheres have been created. With a stereographic 
map projection, these maps will not plot beyond the equator. 

SNCROSS : 

SNCROSS, the sounding cross section program now draws a cross section 
for any me t e r o r o 1 og i c a 1 parmameter that can be computed. In addition, 
a time section is drawn if a single station and multiple times are 
entered . 

COND I T I ONAL PARAMETERS : 

Parameters input in SFPARM and SNPARM may include conditions and 
arithmetic operations. Conditions include the operators <, >, =. 

Arithmetic operations include + , -, *, /, followed by a real number. 
For example, TMPF>100 selects those stations where the temperature 
is greater than or equal to 100 degrees. TMPF*10 multiplies the 
temperature by 10. This replaces the variable, SCALE, in SFMAP 
and SNMAP . 

CENTERED PARAMETERS: 

In SFMAP and SNMAP, a parameter is plotted centered at each station 
if it is the last parameter in the list and an @ is appended. 

SURFACE DECODER: 

The RA library, which decodes surface airways data, is included. 
NON-TAE INTERFACE: 

The following new features have been added to the user interface 
used when running programs outside the TAE : 

- The user will be prompted for parameters before the 
program is executed. Previously, the program was 
executed with current global values first. 

- A DISPLAY function has been added which displays 
the level 1 help along with values for the program 
par ame t e r s . 

- SAVE and RESTORE functions have been added. 

- If $RESPOND is NO, the first tutor will still be 
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displayed. 

- IF $MAPFIL is changed, the new map file will be 
activated immediately. In the TAE and in earlier 
versions of the non-TAE, the value of SMAPFIL was 
not changed until after the program was started 
again. 

MULTIPLE FILE INPUT: 

The objective analysis programs can read data from two files. This 
allows ship data to be included in an objective analysis. 

Multiple grid files can also be used in the grid diagnostic programs. 
This allows easy comparisons between two different models, although 
they must have the same navigations. 

STATION REPORT TIME: 

The station report time, STIM, is now a station parameter. The time 
will be printed as hours and minutes in the form HHMM. 


BACKGROUND PROFILE LINES: 

A minimum, maximum, and increment can be specified for the 
background lines on a profile in THTALN , THTELN , and MIXRLN . 

UNMERGED DATA: 

Da t a wh ich have not been me rged, i.e., data stored at mand atory and 
significant levels separately, can be listed in the separate parts 
in SNLIST, read from a text file in SNEDIT, and transferred to a new 
file in SNMOD. The choice between merged and unmerged data is 
specified in the par ame ter, MRGDAT . 


LAYER QUANTITIES: 

Several layer quantities can be computed in the sounding programs 
i nc 1 ud i ng : 


RICH 

BVFQ 

BVSQ 

BVPD 

LAPS 

STAB 

STAP 

SHRD 

SHRM 

SEPA 


Richardson number 
B r un t - Va i s a 1 a frequency 
B r un t - Va i s a 1 a frequency squared 
B r un t -Va i s a 1 a period 
T emp erature lapse rate 
Potential temperature lapse rate 
Potential temperature change with pressure 
Wind shear direction 
Wind shear magnitude 

Pressure difference in isentropic layer 


STATION DIAGNOSTIC QUANTITIES: 

Several stability indicies can be computed at sounding stations. 
These new quantities, which can be specified in STNDEX , include: 


BRCH - Bulk Richardson number 

CAPE - Convective Available Potential Energy 
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